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Fronts that start from a local perturbation and propagate into a linearly unstable state come in 
two classes: pulled fronts and pushed fronts. The term "pulled front" expresses that these fronts are 
"pulled along" by the spreading of linear perturbations about the unstable state. Accordingly, their 
asymptotic speed v* equals the spreading speed of perturbations whose dynamics is governed by 
the equations linearized about the unstable state. The central result of this paper is the analysis of 
the convergence of asymptotically uniformly traveling pulled fronts towards v* . We show that when 
such fronts evolve from "sufficiently steep" initial conditions, which initially decay faster than e~ x x 
for x — > oo, they have a universal relaxation behavior as time t — > oo: The velocity of a pulled front 
always relaxes algebraically like v(t) = v* - 3/(2A*i) + (30F/2) DA* /(DA* 2 i) 3/2 + 0(l/t 2 ). The 
parameters v* , A* , and D are determined through a saddle point analysis from the equation of motion 
linearized about the unstable invaded state. This front velocity is independent of the precise value 
of the front amplitude which one tracks to measure the front position. The interior of the front is 

essentially slaved to the leading edge, and develops universally as <f>(x, t) =*„ (t) (x-fdt' »(0) + 

0(l/t 2 ), where & v (x — vt) is a uniformly translating front solution with velocity v < v*. Our result, 
which can be viewed as a general center manifold result for pulled front propagation, is derived in 
detail for the well known nonlinear diffusion equation of type dt4> = d 2 4> + (j) — <j) Z , where the invaded 
unstable state is 4> = 0. Even for this simple case, the subdominant t~ 3 ^ 2 term extends an earlier 
result of Bramson. Our analysis is then generalized to more general (sets of) partial differential 
equations with higher spatial or temporal derivatives, to p.d.e.'s with memory kernels, and also to 
difference equations such as those that occur in numerical finite difference codes. Our universal 
result for pulled fronts thus implies independence (i) of the level curve which is used to track the 
front position, (ii) of the precise nonlinearities, (Hi) of the precise form of the linear operators in the 
dynamical equation, and (iv) of the precise initial conditions, as long as they are sufficiently steep. 
The only remainders of the explicit form of the dynamical equation are the nonlinear solutions 
and the three saddle point parameters v* , A*, and D. As our simulations confirm all our analytical 
predictions in every detail, it can be concluded that we have a complete analytical understanding of 
the propagation mechanism and relaxation behavior of pulled fronts, if they are uniformly translating 
for t — > oo. An immediate consequence of the slow algebraic relaxation is that the standard moving 
boundary approximation breaks down for weakly curved pulled fronts in two or three dimensions. 
In addition to our main result for pulled fronts, we also discuss the propagation and convergence of 
fronts emerging from initial conditions which are not steep, as well as of pushed fronts. The latter 
relax exponentially fast to their asymptotic speed. 
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I. INTRODUCTION 

A. Outline of the problem 

In this paper we address the rate of convergence or 
"relaxation" of the velocity and profile of a front that 
propagates into an unstable state. The particular fronts 
we analyze separate two nonequilibrium homogeneous 
states, one of which is stable and one of which is un- 
stable, and are such that the asymptotic front solution 
is a uniformly translating one. We assume that the un- 
stable state is initially completely unperturbed in a large 
part of space, and that thermal and other noise are neg- 
ligible. Examples of such situations arise in one form or 
another in physics [fj]-f29[ , chemistry |29|-|35|] , and biology 
p0| , [32|]36| . If the unstable state domain is not perturbed 
by imperfect initial conditions or thermal noise, it can 
only disappear through invasion by the stable state do- 
main. We analyze the propagation of fronts formed in 
this process, in particular the temporal convergence to- 
wards an asymptotic front shape and velocity, and show 
that it is characterized by a universal power law behav- 
ior in the so-called pulled regime. We concentrate on 
planar fronts, which thus can be represented in one spa- 
tial dimension. However, our results for these and for the 
dynamical mechanism also have important implications 
p7| for the derivation of moving boundary approxima- 
tions f3§| , |39| for weakly curved fronts in higher dimen- 
sions, as well as for the evaluation of the effects of noise 
on fronts jffl-Q , especially the effect of multiplicative 
noise 

The problem of front propagation into an unstable 
state has a long history, which dates back |49| to the pio- 
neering work by Kolmogoroff, Petrovsky and Piscounoff 
(= KPP) |(| and by Fisher @ on the nonlinear diffu- 
sion equation 



= %<t> + f(<j>) , 



(1.1) 



where f(4>) is such that it has a homogeneous stable 
state <fi = 1 and a homogeneous unstable state <fi = 0. 
The early work on this equation p0| , pl| was motivated 
by the biological problem of gene spreading in a popu- 
lation. Since this work, the nonlinear diffusion equation 
( |l . l[ ) , in particular the one with a simple nonlinearity of 
the type 

/ = /kpp(0) = <t> - k , k > 1, e.g., k = 2 or 3 , (1.2) 

has become a standard problem in the mathematical lit- 
erature [|0|j3|,|^,||-§|. For the F-KPP equation de- 
fined by ( |1 . 1[ ) and Ql.2| ) , there exist dynamically stable 
uniformly translating front solutions 4>(x, t) = <fr v (x — vt) 
for every velocity v > v* — 1\J /'(0), and hence every one 
of these solutions is a possible attractor of the dynam- 
ics for long times t. The resulting dynamical behavior or 
"velocity selection" depends on the initial conditions and 
has been investigated by a variety of methods KG p3,p4[ ; 



essentially all its relevant properties have been derived 
rigorously |^2| . E.g., following the lines of KPP [|50|| , 
Aronson and Weinberger proved rigorously |Q , that ev- 
ery initial condition, that decays spatially at least as fast 
as e~ A x (A* = v* /2) into the unstable state for x — » oo, 
approaches for large times the front solution (x — v*t) 
with the smallest possible velocity v* . Most of the rig- 
orous mathematical methods can, however, not be ex- 
tended to higher order equations 

In physics, the interest in front propagation into un- 
stable states initially arose from a different angle. Since 
the late fifties, the growth and advection of linear per- 
turbations about a homogeneous unstable state has been 
analyzed through an asymptotic large time expansion of 
the Green's function of the linear equations |58|-|60| . Only 
ten to fifteen years ago did it become fully clear in the 
physics community |pl|-[7l[|, that there was actually an 
empirical but deep connection between the rigorous re- 
sults for the second order equations and some aspects of 
the more general and exact but nonrigorous results for 
the growth of linear perturbations. This has given rise to 
a number of reformulations and intuitive scenarios aimed 
at understanding the general fro nt p ropagation problem 
into unstable states |^ , |63] , |65| , |67| -[7lll . 

Although our results bear on many of these ap- 
proaches, our aim is not to introduce another intuitive 
or speculative scenario. Rather, we will introduce what 
we believe to be the first systematic analysis of the rate of 
convergence or "relaxation" of the front velocity and pro- 
file in the so-called "linear marginal stability" (6^^^] or 
"pulled" fn| , f70| , f7l| regime. In this regime the asymptotic 
front velocity is simply the linear spreading speed deter- 
mined by the Green's function of the linearized equations. 
Quite surprisingly, our analysis even yields a number of 
new and exact results for the celebrated nonlinear diffu- 



sion equation (1.1), but it applies equally well to (sets 
of) higher order partial differential equations that admit 
uniformly translating fronts, to difference equations, or 
to integro-differential equations. We will discuss such 
equations in general, and then illustrate our results on 
the example equations from Table I. 

For all such equations, our results have a remarkable 
degree of simplicity and universality: Pulled fronts al- 
ways converge in time with universal power laws and 
prefactors that are independent of the precise form of 
the equations and independent of the precise initial con- 
ditions as long as they obey a certain steepness crite- 
rion. To be precise, for equations such that the dynam- 
ically selected asymptotic front is a uniformly translat- 
ing pulled front, and for so-called sufficiently steep initial 
conditions defined such that lim^oo </>(x, 0) e Xx = for 
some A > A* , we derive that the asymptotic velocity con- 
vergence is given by the universal law 



v(t) 
X = — 



X 



2X*t 



1 - 



(X*) 2 Dt 



(1.3) 
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The velocity v*, the inverse length A* and the diffu- 
sion constant D are in general obtained from a saddle 
point expansion J60[ for the equation of motion linearized 
about the unstable state. In a frame moving with veloc- 
ity v* the quickest growing mode k* is identified by the 
complex saddle point equation dk [oj(k) — v*k]\ k=k , = 
where u(k) is the dispersion relation of a Fourier mode 
In the more usual decomposition into real 



-itjt+ikx 



functions this implies that p8 63 6 



91m 



91m k 



dim 



k' 



dRe k 



= 



(1.4) 



The speed of the frame is asymptotically the same as the 
speed of the front if 



Im Lu(k*) 
Im k* 



(1.5) 



For the uniformly translating fronts that we will analyze 
here, we have 

Im k* = A* > , Re k* = , Re uo(k*) = , (1.6) 

and a real positive diffusion coefficient D 



D 



id 2 w 



2dk 2 



d 2 lm uj 



2d\ 2 



(1.7) 



While the velocity of a front is converging to its asymp- 
totic value, so is the profile shape. Note that v(t) dLq ) 
does not depend on the "height" </> = h, which is being 
tracked. In fact, if we define the velocity of the fixed 

amplitude <f> = h through Cj> (x + f dT 

v <t>( T ) > t ) = K 



then up to order 1/t 2 the velocity v^t) 



pendent of the "height" 4> = h. Moreover, it is determined 
solely by properties of the equation linearized about the 
unstable state, as Eqs. (|F4|)-([L7]) show. In this sense, 
we can indeed speak of pulling of the front by the leading 
edge of the front. 

The above expression for v(t) contains all the univer- 
sal terms, since the next 1/t 2 term in the long time ex- 
pansion does depend on initial conditions. The above 
analytic results for the universal velocity convergence as 
well as related ones for the relaxation of the front profile 
which arc summarized in Table II and discussed in more 
detail below, are fully confirmed by extremely precise nu- 
merical simulations. Taken together, this study therefore 
yields the understanding of the pulled front mechanism 
that so many authors P,p3 65 -67,p$-f72j have sought for. 



In this paper, first the asymptotic long time behavior 
is worked out in det ail and to high orders for the F-KPP 
equation (1.1), (1-2) in two matched asymptotic expan- 
sions in \j\Jt. Once we will have laid out the struc- 
ture of this expansion, it is clear that essentially the 
same matched expansions can be applied to other more 
complicated types of equations, provided that they ad- 
mit a family of uniformly translating front solutions in 



the neighborhood of the asymptotic "pulled" velocity v* . 
Moreover, the two lowest order equations in the 1/yt 
expansion in the so-called leading edge region together 
with a boundary condition suffice to calculate the uni- 
versal convergence. The structure of these equations is 
virtually independent of the precise form of the dynami- 
cal equation. For more general equations, 



The EFK eq.: 

d t 4> = -jd^ + d 2 ^ + /(</>) with < 7 < 1/12 , 

The streamer eqs.: 

d t a = Ddla + d x {aE) + a\E\e- x l\ E \ , 
d t E = -Dd x a - oE , 

A difference equation from kinetic theory: 
d t C t (t) = -C z + Cl 1 , 

A second order extension F-KPP eq.: 

T 2 d 2 <p + d t( f> = d 2 x( f> + 0-(!> 3 , 

An equation with a memory kernel: 

d t <f>(x, t) = d 2 0(x, t) + /* dt' K(t - t 1 ) <f>(x, t') 
-0 fc OM) , 

Finite difference versions of the F-KPP eq. like: 

Uj(t + At) - Uj(t) _ u j+ x(t) - 2uj(t) + Uj-i(t) 



At 



(Ax) 2 
+ u 3 {t)-u k 3 (t) . 



■X is inde- Table I: Summary of the equations studied in detail in 



Section V F as examples of the general validity of our re- 
sults for higher order equations, coupled equations, dif- 
ference equations, and equations with a kernel. All these 
equations have pulled front solutions whose asymptotic 
speed relaxes according to (1.3). 



we hence limit the discussion to the motivation and anal- 
ysis of these two equations. Although we will give some 
discussion of the assumptions that underly the expansion 
(like the one that there is a nearby family of moving front 
solutions), a full analysis of these as well as of the ex- 
tension to nonuniformly translating fronts, such as those 
arising in the EFK equation of Table I for 7 > 1/12, 
in the Swift-Hohenberg equation |f73| , or in the com- 
plex Ginzburg-Landau equation |6^] , will be left to future 
publications [[74[-f76|. 

For equations fll.l| ), (1.1 ) we sim ply have v* = 2, 
A* = D = 1. The first term in (|1.3 |) then reduces 
to a wellknown result of Bramson [[77| , who rigorously 
proved that the convergence to the asymptotic velocity 
v* is v(t) = v* — 3/(2X*t) uniformly, i.e., independent of 
the amplitude <f> whose position one tracks. The factor 
3/2 in this expression has often been considered puzzling, 
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since the linear diffusion equation with localized initial 
conditions yields v(t) = v* - l/{2X*t) + ... . In 0, it 
was argued that the factor 3/2 in this result applies more 
generally to higher order equations as well, but a system- 
atic analysis or an argument for why the convergence is 
uniform, was missing. Apart from this and a recent red- 
erivation j72| of Bramson's result along lines similar in 
spirit to ours 1 and a few papers similar in spirit to that 
of Bramson [|^J7^J8^], we are not aware of systematic 
calculations of the velocity and profile relaxation. Even 
for the convergence of the velocity in the celebrated non- 
linear difusion equation, our 1/i 3 / 2 term appears to be 
new. 

From a different perspective, Powell et al. |39| also 
considered the convergence properties of pulled fronts. 
These authors studied the shapes of the front profiles in 
the nonlinear diffusion equation and argued that they re- 
lax along the the family of unstable uniformly translating 
front solutions. Although they realized that the velocity 
relaxation was algebraic and from below, they did not 
seem to realize that the dominant — 3/(2t) velocity cor- 
rection was known from earlier work | J77] , p7[ . As we shall 
see below when we will discuss the shape relaxation of 
fronts, our derivation is the first analytic derivation and 
confirmation of the picture of Powell et al., and identifies 
the connection with the velocity relaxation. 

Our results are not only of interest in their own right, 
but they also have important implications. Since the 
asymptotic convergence towards the attractor $* is al- 
gebraic in time, the attractor alone might not give suf- 
ficient information about the front after long but finite 
times, since algebraic convergence has no characteristic 
time scale. In particular, there is no time beyond which 
convergence can be neglected. Such slow convergence 
means that in many cases, experimentally as well as the- 
oretically, one observes transients and not the asymptotic 
behavior. In fact, in the very first explicit experimental 
test of front propagation into unstable states in a pat- 
tern forming system ||, viz. Taylor-Couette flow, the 
initial discrepancy between theory and experiment was 
later shown to be related to the existence of slow tran- 



sients 
oretica 



1( 



The slow convergence is important for the- 
studies as well: it is a common experience (see, 
e.g., |Jl2|,^6 81 1) that when studying front propagation in 
the "pulled" regime numerically, the measured transient 
front velocity is often below v*. This is so even though 
the asymptotic front speed can never be below v* , be- 
cause no slower attractor of the dynamics exists. This 



observation finds a natural explanation in our finding 
that the rate of convergence is always power law slow, 
and that the front speed is always approached from be- 
low. 

A second important implication of the absence of an 
intrinsic time scale of the front convergence is the follow- 
ing. When we consider the propagation of such fronts 
in more than one dimension in which there is a coupling 
to another slow field (as, e.g., in the phase field mod- 
els p9| , p2| , ^3| ) , the front dynamics does not adiabatically 
decouple from the dynamics of the other field and from 
the evolution of the curvature and shape of the front it- 
self. This implies that the standard moving boundary 
approximation |38],[39 84 1 (which actually rests on the as- 
sumption that the convergence on the "inner scale" is 
exponential) can not be made. Though this is intuitively 
quite obvious from the power law behavior of the front 
convergence process, the connection between the conver- 
gence and the breakdown of a moving boundary approx- 
imation also emerges at a technical level: the divergence 
of the solvability integrals that emerge when deriving a 
moving boundary approximation turns out to be related 
to the continuity of the stability spectrum of pulled fronts 
p7[ . The break-down of the solvability analysis for per- 
turbations of the asymptotic front in the pulled regime 
also has consequences for the evaluation of multiplicative 
noise in such equations |37|| 



B. Pushed versus pulled fronts, selection and 
convergence 

Let us return to the well understood nonlinear diffu- 



sion equation (1.1) and discuss to which nonlinearities 
f(<p) our prediction of algebraic convergence applies and 
why. If /'(0) < 0, the invaded state 4> = is linearly sta- 
ble, and the construction of a uniformly translating front 
4>(x,t) — $> v (x — vt) poses a nonlinear eigenvalue prob- 
lem. The solution with the largest eigenvalue v is the 
unique stable and dynamically relevant solution (unique 
up to a translation, of course). As is well known and dis- 
cussed in Section ||, any initial front that separates the 
(meta)stable state <j) = at x — > oo from another stable 
state at x — > — oo will converge exponentially in time to 
this unique attractor However, whenever /'(0) > 0, 
4> = is unstable, there is not a unique asymptotic at- 
tractor <!>.„, but a continuous spectrum of nonlinear eigen- 



x The main focus of the work by Brunet and Derrida j72| is 
actually the correction to the asymptotic velocity if the func- 
tion f(4») has a cutoff h such that fh{<t>) = for cj> < h. The 
method the authors use to derive this, is actually closely re- 
lated to the one they use to rederive Bramson's result, and to 
our approach. See in this connection also the recent paper by 
Kessler et al. Q. 
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values v which constitute the velocities of possible attrac- 
tors $„. The existence of a continuum of attractors of 
the dynamics poses a so-called selection problem: From 
which initial conditions will the front dynamically ap- 
proach which attractor? The attractor with the smallest 
velocity plays a special role, as its basin of attraction are 
all "sufficiently steep" initial conditions, as defined in 
Section ||. It therefore will be refered to as the selected 
front solution. 

When we concentrate on these "sufficiently steep" ini- 
tial conditions and analyze the dependence on the non- 
linearity / in (1.1), the transition from exponential to 
algebraic convergence does not coincide with the tran- 
sition from stability to instability of the invaded state 
(j) = 0, but with the transition between two different 
mechanisms of front propagation into unstable states. In- 
deed, it is known (see also Sect. |l]), that for /'(0) > 0, 
there are two different mechanisms for how the selected 
front & se i and its speed v se i are determined. Either & se i 
is found by constructing a so-called strongly heteroclinic 
orbit for $„ from the full nonlinear equation. This case 
is also known as case II [ |63| or nonlinear marginal sta- 
bility (6|@, or as pushing ||l|,^o|j7l) . Or, the selected 
velocity v se i is determined by linearizing about the un- 
stable state <fi — 0, which case is known as case I or linear 
marginal stability, or as pulling. We henceforth will use 
the terms "pushing" and "pulling" for the two different 
propagation mechanisms of a selected front evolving from 
steep initial conditions, since they very literally express 
the different dynamical mechanisms. 

In a pushed front just like in a front propagating into 
a (meta)stable state, the dynamics is essentially deter- 
mined in the nonlinear "interior part" of the front, where 
(f> varies from close to (j) = to close to the stable 
state. The construction of the selected front as a strongly 
heteroclinic orbit in the pushed case continuously ex- 
tends into the construction of the heteroclinic orbit of 
the unique attractor if the invaded state is (meta)stable 
(/'(0) < 0). For both pushed fronts and fronts propa- 
gating into linearly stable states, the spectrum of linear 
perturbations is bounded away from zero, so that con- 
vergence towards the asymptotic front is exponential in 
time. 

In a pulled front, the dynamics is quite different: As 
we shall see, it is determined essentially in the region 
linearized about the unstable state. We call this region 
the leading edge of the front. Eq. ([hi]) is appropriate for 
analyzing the front interior. We will see in Section [ID, 
that a stability analysis performed in this representation 
is not able to capture the convergence of a steep initial 



condition towards a pulled front. Rather the substitution 

f/- = ^e^ , i = x-v*t , (1.8) 

which we s hall term the leading edge representation, 
transforms ([hi]) into 



9tif> 
f 



(1.9) 



o 



(v 



/(V 



2 e -A*« 



/'(0) 



This equation will turn out to be appropriate for analyz- 
ing a leading edge dominated dynamics. Note that / is at 
least of order tp 2 with an exponentially small coefficient 
as £ — ► oo . For large £, the dynamics is purely diffusive. 
If the nonlinearity obeys f(<fi) — f'{0)4> < for all <fr > — 
which is known as a sufficient criterion for pulling — the 
nonlinearity / is always negative. Then / purely damps 
the dynamics in the region of smaller £. The dynamics 
evolving under (l.£) is equivalent to simply linearizing 
(1.1) about the unstable state in the large £ region - 
there is only one subtle but important ingredient from 
the requirement that the dynamics in the linear region 
crosses over smoothly to the nonlinear front behavior at 
smaller £, that actually enters our leading edge analysis 
in the form of a bou ndary condition. In the leading edge 
representation (|l.9| ), this is brought out by the presence 
of the sink-type term / which is nonzero in a localized 
region behind the leading edge. With this small caveat 2 , 
we can conclude that the leading edge of the front "pulls 
the rest of the front along" , which is precisely the mecha- 
nism that gives rise to the universal algebraic convergence 
behavior. In a pushed front, in contrast, the nonlinear- 
ity "pushes the leading edge forward" and convergence is 
exponential. 

To illustrate this discussion by a concrete example, we 
note that when the function f(cj>) in the nonlinear diffu- 
sion equation is of the form 



LU + l 



L2n+1 



n>0 



(1.10) 



we can rely on known analytic solutions for $ v . In this 
case, the state <\> = is (meta)stable for e < 0. For 
< e < (n + l)/n 2 , the selected front is pushed, and for 



e > (?i + l)/n' 2 , it is pulled (see Section || and Appendix 

S)- . ' . 

At this point, a brief explanation of our use of the word 
"metastable" may be appropriate. For systems with a 
Lyapunov function, the word metastable is often used in 
physics to denote a linearly stable state, which does not 
correspond to the absolute minimum of the Lyapunov 



2 Note though, that this subtle point is quite important — 
as we shall see, the saddle point or pinch point analysis gives 
precisely the wrong prefactor for the leading 1/t convergence 
term because this boundary condition is not satisfied. 
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function or "free energy" . A domain wall or front be- 
tween the absolutely stable and a metastable state then 
moves into the metastable domain; one may therefore 
loosely call a linearly stable state "metastable", if it is 
invaded by another "more stable" state through the mo- 
tion of a domain wall or front. 

The understanding of the two different dynamical 
mechanisms of pushing and pulling in the nonlinear diffu- 
sion equation ( |l.l[ ) lays the basis for the analysis of equa- 
tions like those listed in Table I. The essential step to- 
wards a generalization of the leading edge representation 
(1.9) is done by a saddle point analysis, that identifies 
which Fourier modes of linear perturbations of the un- 
stable state will dominate the long time dynamics. This 
analysis yields the parameters v* , A*, the diffusion con- 
stant D and possible higher order terms required for the 
leading edge representation. 



C. Sketch of method and results on front relaxation 
in the pulled regime 



Bramson's method |77| to calculate algebraic conver- 
gence is specifically adapted to equations of type (1.1). 
It is based on a representation of the diffusion equation 
by Brownian processes, which are evaluated probabilisti- 
cally. Instead, we construct the asymptotic convergence 
trajectory towards a known asymptotic state by solving 
the differential equations in a systematic asymptotic ex- 
pansion which, though nonrigorous, extends immediately 



to higher order equations. Our approach leads to exact 
results since the expansion parameter are inverse powers 
of the time t, so these terms become arbitrarily small in 
the asymptotic regime. 

The idea of the method is that in a pulled front, the 
speed is essentially set in the leading edge, where lin- 
earization of the equation of motion about the unstable 
state is justified. This leading edge has to be connected 
to what we will refer to as the interior part of the front, 
defined to be the region where we have to work with the 
full nonlinear equation. For the interior, we use the fact 
that for large times the shape of the converging front 
will resemble the asymptotic front, and thus can be ex- 
panded about it. We also explicitly make use of the fact 
that the initial state <fi(x, 0) for large x is steeper than 
the asymptotic front profile <&* = in the leading edge, 
i.e., 4>(x, 0) e x x — > as x —* oo. The structure of the 
problem then dictates the expansion in \j\ft. 



The structure of the expansion in 1/s/t is the only real 
input of the analysis; its selfconsistency becomes clear a 
posteriori and it can be motivated from the earlier work 
on the long time expansion of the Green's function of the 
linearized equations. Equivalently, the self-consistency 
emerges from the observation that the equation govern- 
i ng t he convergence towards the asymptotic front profile 
(1.9) reduces essentially to a diffusion equation in the 
leading edge of the front. The derivation of the exact re- 
sults summarized in Table II is essentially based on this 
ansatz. 
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FIG. 1. Illustration of the fact, that even though the shape of a front profile is quite close to $*, the position of a front is 
shifted logarithmically in time relative to the uniformly translatin g pr ofile $*(£)• 

Solid lines: evolution of some initial condition <f>(x, 0) of the form (fif3) under dt(f> = di<f> + (j> — <f> A at times t — 5, 10, 15. Dotted 
lines: evolution of <f>(x,t) = $*(£), £ = x— 2t, at times t — 5, 10, 15. $* is placed such, that the amplitude 0=1/2 coincides 
with that of <f>(x,i) at time t = 5. The logarithmic temporal shift is indicated by the fat line. 



The shape convergence is also obtained explicitly from 
our analysis. The crucial input for the analysis is the 
right frame and structure to linearize about. At first 
sight, a natural guess would be that for large times, the 
actual shape of the front <j>(x, t) should be linearizable 
about the shape of the asymptotic front <j>*( x — v*t). 



However, the algebraic velocity convergence (1.3) im- 
plies, that if a converging front profile <f> is close to the 



asymptotic uniformly translating front profile $* (x— v*t) 
at some time t$, the distance between the actual pro- 
file and <£* will diverge at large times t as X(t) = 
-(3/2A*)ln(i/*o) + • ■ ■• This result which is illustrated 
in Fig. 1 implies that if we want to linearize (f> about $* 
at all times, we have to move $* along with the non- 



asymptotic velocity v(t) (1.5) of the converging front. A 
crucial step for the analysis is thus to linearize about 



8 



$*(£x) in a coordinate system 

(fr,t) , &-=£-X(t) = a;-i;*t-X(f) , (1.11) 

moving with the conv ergin g front. If we expand cj> about 
^*(£x) with !;x from (1.11) and then resum, we find that 
the interior shape of the front is given by 

<j>(x,t) =® v (t)(tx+x )+o(J^) (1.12) 



for £x *C V^Dt. xq expresses the translational degree 
of freedom of the front. The uniformly translating front 
$„(£) is a solution of the ordinary differential equation 
for the uniformly translating profile cj)(x, t) = $„(x — vt) 
but with v replaced by the instantaneous value v(t) o f th e 
velocity. E.g., for the nonlinear diffusion equation (1.1), 
$„(£) is the solution of 



vd^ v (0 = £>f<MO + /(MO) 



(1.13) 



Eq. ( 1.12] ) also confirms that to leading order the interior 
is slaved to the slow dyna mics of the leading edge. The 
transient profiles & v (t) i n ( 1 -12| ) propagate with velocity 
v(t) smaller than v* according to ( |1.3| ). 

For the special case of Eq. flLl| ) it is well known (see 
also Section ||), that when constructing a front & v start- 
ing from = 1 at £ — > — oo, it eventually will become 
negative for finite £ whenever v < v* , and that globally 
such fronts either do not exist or are dynamically un- 
stable, depending on the properties of / for negative </>. 
However, only the positive part of &v(t) from £ — > — oo 
up to £ -C s/t plays a role as a transient. That the con- 
vergence trajectory is approximately given by $> v m, was 



already observed numerically in equations of type (1.1) 
by Powell et al. [[39]]. Our analytical derivation of this 
result actually holds for a larger class of equations, but 
at the same time we find that it only holds up to a cor- 
rection term of order l/t 2 . This non-universal correction 
is always non-vanishing. 

For £x 3> V^Dt, the transient crosses over to 



(1.14) 



The analytical expre ssion for the unive rsal c orrect ion o f 
order \jlft in fll.l4j ) is given by Eqs. (|3 65|) and (|3 67|) 
for th e non line ar di ffusion equation and is generalized by 
Eqs. ( 5.39 ) or (|5.69|) , while the correction of order l/t will 
depend on initial conditions, and is thus non-universal. 

A crucial insight implemented above is that the front 
consists of different dynamical regions which have to be 
matched to each other. The situation is sket ched in Fig. 
2. For a pulled front, the Gaussian region ( 1.14 ) of the 
leading edge e ssent ially determines the velocity while the 
front interior ( 1.12 ) is slaved to leading order. The Gaus- 
sian region might be preceeded by a region of "steep- 
ness" A being conserved in time which for sufficiently 
steep initial conditions A > A* has no dynamical im- 
portance (where the steepness A is defined in Eq. ( |2.6| ) 
below). Likewise, for flat initial conditions the dynam- 
ics is dominated by the conserved A region, while pushed 
dynamics is dominated by the front interior. In both of 
these cases, the intermediate Gaussia n re gion is absent. 
For the nonlinear diffusion equation (1.1), the different 
cases are discussed in Section [n] and summarized in Table 
IV. 



Our results (E3) - ( 1.14 ) are universal in four ways: 

• They are independent of which "height" or level curve 
is being tracked to define the front velocity. 

• The predicted convergence behavior is independent of 
the precise initial conditions, provided they decay quicker 
than e~ A ' x ' far in the unstable regime. 
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• The leading edge behavior (1.3) and ( 1.14 ) is indepen- 
dent of the precise nonlinearities. For Eq. (1.1), the con- 
stants v* , A* and D depend on /'(0) only. For the more 
general equations, these constants are completely deter- 
mined by the saddle point expansion in the equation lin- 
earized about the unstable state. 

• If we analyze general equations like those listed in Table 
I, our prediction for the interior part of the front ( 1.12j ) 
stays unchanged, as long as the front speed stays deter- 
mined by the linearization about the unstable state, i.e., 
the front stays pulled, and as long as the state behind the 
front stays homogeneous. The effect of the nonlinearities 
just gets absorbed in appropriate functions $.„. 

The results summarized in this Subsection are the most 
central new results of this paper. They are summarized, 
for easy reference in Table II. 





intprinr 


leading edge: 






Gaussian region conserved A, region 




<i> ~ * v(t) (^xKs 


^J> ~ % e** e^ /4Dt <> ~ e XiK - vm 



X 







FIG. 2. Sketch of a pulled front with the different dynamical regions: interior = nonlinear region, leading edge = region 
linearized about the unstable state. Depending on the initial conditions, the leading edge might still consist of two different 
regions: a Gaussian region and a region of conserved steepness A. For A > A* (defining "sufficiently steep" initial conditions), 
the (intermediate) asymptotic Gaussian region determines the velocity relaxation. 



Height independent velocity: 
v(t) = v* +X 



3 <, /=Z= +0 I 



2X*t V V {^*) 2 Dt ) \t 2 
where the saddle point analysis of the 
linearized equation yields v* , A*, D, 



cf. Table V or Eqs. (1.4)-(|l.7|) 



=> Use the coordinate: £x — x ~ v*t — X(t) 
Front for £x <C V <ADt (front interior) : 

(x,t) = *„ (t )(60 + o(^) , 



**(£x)+.*»fch(£x) + 



1 

where $ V (Q solves <j>(x,t) = & v (x — vt) 
and r] ah (£) = 8<$> v (£)/5v\ v , . 



Front for £x ^ V 'ADt (leading edge) : 
,t) = ctxe-W* e-&/(«*) + _ . 



Table II: The central results on the universal algebraic relaxation towards uniformly translating pulled 
fronts, see also Fig. 2. These results apply to steep initial conditions in the nonlinear diffusion equation in 
the pulled regime (Case IV of Table IV, see Sections III and IV) and to more general equations (see Section 



D. Organization of the paper 



Before embarking on our explicit calculation of the ve- 
locity and shape convergence in the pulled regime, we 
review in Section || rather well-known results on the 
multiplicity, stability and convergence of pushed fronts 
in the nonlinear diffusion equation, and discuss how far 
these results can be extended to pulled fronts or fronts 
emerging from "flat" initial conditions. Since the conver- 
gence towards pulled fronts cannot be derived by linear 
stability analysis, we set the stage for Section III by in- 
troducing the leading edge transformation. In the central 
Section III the detailed analysis of pulle d fr ont relaxation 
in the the nonlinear diffusion equation (1.1 ) is given. The 
detailed numerical simulations that fully confirm our an- 
alytical predictions are presented in Section IV. In this 



spatial discretization and system size arising in the nu- 
merical solution of pulled front propagation. In Section 
^ we extend our analysis to more general equations, dis- 
cuss the example equations listed in Table I, and present 
numerical results, again in excellent agreement with our 
analytical predictions. Here the picture of a new center 
manifold theorem for pulled front propagation emerges. 
We then close the main body of the paper with a sum- 



mary and outlook in Section VI 



Section we also pay attention to the specific problems of 



Since this is a long paper with a large number of de- 
tailed results of various types, and since we have made an 
attempt to make our results accessible for readers from 
different fields, we introduce Table III as a "helpdesk" for 
the reader who wants to focus on a particular aspect of 
the front propagation problem only, or who wants to get 
only an idea of the essential ingredients of our approach 
and the main results. 



10 



We finally note that a brief sketch of our results can be 
found in [B5| and the lecture notes [pal . Later extensions 



of the present work can be found in 



THE READER'S HELPDESK 



If- 



Th 



en our a 



dvice 



you want to focus right away on the relax- 
ation calculation without analyzing what a 
stability calculation tells and does not tell 
about the relaxation of pulled fronts 



if you know what is meant with the "pulled" 
velocity v* you can start with S ectio n 111 
immediately; if not, read Section JIF4 flF 



and possibly Section V C 1 first 



you are already familiar with previous ideas 
concerning front selection in the physics lit- 
erature, but want to get an idea of our 
change of emphasis and of the new detailed 
results of this paper 



to skim Section II A for notation and a sum- 



mary of the most important background ma- 
terial, to then check the appropriate Sections 
of on p oints which are not clear from Sec- 
tion IIA, and then to proceed to Section III 



you (mainly) want to read about the con- 
nection between stability, selection and 
relaxation 



to read Section |l| with Table IV and for 
the generalization Sections VB and V C with 
their appendices 



you only want to get an idea of the concep- 
tual basis of the algebraic convergence 



to re ad Secti on III A and possibly Sec- 
tions |VC] - |VE| for the arguments concerning 
higher order partial differential equations or 
other types of equations 



you are unfamiliar with the concept of pulled 
velocity v* for higher order equations and 
want to know how it is determined 



[LTa| an< 
ions|VD 



to read Sections III A^ and V C_l J_and possi- 



bly parts of Sections |V D and VE ) 



you just want to sec the numerical support 
for the algebraic relaxation prediction from 
Tables II and V, or want to read about the 
numerical intricacies of studying the pulled 
front convergence 



read Section IV on the nonlinear diffusion 



equation and Section VF for higher order 
and coupled equations 



you just want a toolkit for when to apply the 
predictions from Tables II and V 



to read Section VI Ej 



Table III: A guide through the paper for the efficient reader who wants to read about specific results only, 
or who already has some background knowledge on the problem of front propagation into unstable states. 
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II. STABILITY, SELECTION AND 
CONVERGENCE IN THE NONLINEAR 
DIFFUSION EQUATION 

In this Section, we provide the necessary back- 
ground information on fronts propagating into un- 
stable states by reviewing a number of results on 
the multiplicity and stability of uniformly translat- 
ing front solutions of the nonlinear diffusion equation 

llllj3|^||-||^||||3S0ii&ii- We ais ° 

summarize to what extent the linear stability analysis of 
these uniformly translating fronts allows us to solve the 
selection problem, i.e., to determine the basins of attrac- 
tion of these solutions in the space of initial conditions 
and for different nonlinearities /, and to what extent it 
allows us to answer the related question of the conver- 
gence rate and mechanism. It will turn out that the 
linear stability analysis fails to explain how pulled fronts 
emerging from sufficiently steep initial conditions relax to 
their asymptotic speed and profile. This sets the stage 
for a different approach to pulled fronts by introducing 
the leading edge representation. 



A. Statement of problem and essential concepts 

In Sections [n] - |v|, we analyze the nonlinear diffusion 
equation 



dt<t>{x,t) = e£<i> + f(4>) 



(2.1) 



where f((f>) is assumed to be continuous and differen- 
tiable. For studying front propagation into unstable 
states, it is convenient to take 



/(0) = = /(l) , /'(0) = 1 
/((/>) > for all < (j) < 1 



(2.2) 



so that in the interval [0, 1] /(</>) has one unstable s tate 
at 4> = and only one stable state at cj) = 1. Eq. (2.2) 
implies that /'(l) < 0. Note, that we have specified the 
behavior of f(4>) only on the interval < < 1. This 
is all we need since it can be shown by comparison argu- 
ments |52],[)6| that an initial state with < 0(x, 0) < 1 for 
all x co ns erve s this property in time under the dynamics 
of 0, (O 



In passing we note that for a nonlinearity like (1.10) 
a general equation of the form 



d T tp = Dd^v + F e (<p) , F e (0) = = F £ (^ S ) , 

F £ '(0) = e , v? s >0, (2.3) 

results. It allows e to take either sign. For e < 0, the 
state 4> = is linearly stable, for e > 0, it is unstable. 
Fronts propagating into metastable states (e < 0) will 
someti mes also be discussed briefly for co mpa rison. If 
e > 0, (|2.3|) transforms to the normal form (2.1) as 



Fe(<P) 



(2.4) 



t = er , x = y/e/D y , 4> = tp/tfs , f{4>) = 



Hence velocities transform as dx/dt = [dy/dr] /y/De. 

The front propagation problem can now be stated as 
follows. Consider some initial condition < <fi(x,0) < 1 
with 



lim cj)(x, 0) = 



(x, 0) > for some x 



(2.5) 



t hat evolves under the equation of motion (2^2) with 
( |2.2| ) into a front propagating to the right. Which time- 
independent profile and which velocity will this front ap- 
proach asymptotically as time t — > oo, if any? How quick 
will the convergence to this asymptotic front be? Can we 
identify the mechanisms that generate such dynamical 
behavior? Can we rephrase it in such terms that we can 
generalize results to equations other than (2.1)? These 
questions essentially concern the nature of the front se- 
lection mechanism. 

As is well known, the answers to these questions de- 
pend on more specific properties of the initial condition 
as well as of the nonlinearity f(4>). For the nonlinear dif- 
fusion equation, the answer to the selection problem is 
known in full rigor, but we will only review here those 
concepts which are important in a more general context 
and which play a role in the subsequent relaxation anal- 
ysis. We now briefly outline these main concepts and 
results and explain them in more detail in the rest of 
Section II. 

Existence of a family of front solutions. For front 
propagation into unstable states the selection problem 
is different and more intricate than for bistable fronts 
(fronts between two linearly stable states), since when 
one solves the o.d.e. for the uniformly translating pro- 
file <p(x — vt) — $„(?) one finds that there is a family 
of fronts solutions parametrized by the continuous vari- 
able v that are possible attractors of the dynamics. This 
is in contrast to the situation for bistable fronts where 
the selected velocity v is obtained simply as a nonlinear 
eigenvalue problem. 

Steepness of a front. Most of our discussion focusses 
more than earlier work on the central and unifying role 
of the steepness A of the leading edge of a front, defined 
as the asymptotic exponential decay rate: 



4>{x,t) 



x — >oo — \ 7 

~ e 



A = — lim 



d In (, 
dx 



(2.6) 



When (f>(x, t) decays faster than exponentially as x — > oo, 
this implies A = oo. 

Pulled and pushed fronts. The family of uniformly 
translating and dynamically stable fronts <5> v can be 
uniquely parametrized either by the velocity v or by the 
spatial decay rate or steepness A. The difference be- 
tween pushed and pulled solutions is especially clear if we 
characterize them by A. A given nonlinearity / defines 
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two particular steepnesses: A se ; which characterizes the 
pushed and pulled front solutions and X s teep which char- 
acterizes the basin of attraction of these so-called selected 
fronts. The front solution with A = A se ; > 1 defines the 
pushed front, while the pulled one has A = A se ; = A* = 1. 
The continuous family of dynamically stable front solu- 
tions that exists in addition to these selected fronts, is 
parametrized by A < X s teep < 1- The nature and con- 
struction of the fronts is discussed in more detail in Sec- 
tion [IB, together with a simple property of pulled fronts 



which will play an important role in our later relaxation 
analysis, namely the fact that the asymptotic large time 
profile of a pulled front is as <£v (£) ~ £ e~ A 5 for £ ^> 1. 

We will characterize also an initial condition by its 
steepness A and call it a sufficiently steep initial condi- 
tion, if <j>(x, t — 0) decays to zero exponentially faster 
than e~ Asl 



for some A steep < 1, 



sufficiently steep: <f>(x, 0) < e~ Ax , (2.7) 

for some A > A steep , 

otherwise we call it flat: 

flat: 0(x,O)™e- Aa; , X< \ steep , (2.8) 



How A 

Section [II D|. We will see that always < X s t ee p < 1 for 



s teep is determined by /(</>), will be discussed in 



Eq. (2.1), and in particular that for pulled fronts 



pulled fronts: A 



steep 



= A* = 1 



(2.9) 



while for pushed fronts X s teep < 1- The criterion (2.7) 



for steepness includes all initial conditions with bounded 
support or, e.g., the initial condition (j)(x,0) = 9(—x) 
with 9 the step function. 

Note that the intermediate case <p(x, 0) ~ x~ v e x is 
neither sufficiently s teep nor flat, according to our def- 
we shall recover Bramson's 



initions. In Section III 
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observation that such special initial conditions also lead 
to a \jt relaxation of the velocity profile, but with a in- 
dependent prefactor for v < 2. 



Conservation of steepness. In Section II E we discuss 



what we term conservation of steepness: if an initial con- 
dition is characterized by a steepness A, then at any finite 
time the steepness of <fi(x,t) is the same as that of the 
initial condition 4>(x, t — 0). (Note that the limits t — > oo 
and x — > oo do not commute.) 

The linear stability analysis of front solutions can be 
performed in detail for the nonli near diffusion equation. 
As summarized in Section II C, pushed fronts have a 



gapped spectrum, while pulled fronts have a gapless spec- 
trum within their natural Hilbert space. In the selection 
analysis, we in general also need perturbations from out- 
side this Hilbert space. 

Stability and selection. In Section II Dj , we discuss the 
connection between the stability of front solutions and 
the selection mechanism; this connection, which under- 
lies much of the marginal stability scenario paEa,B7[ , 



hinges on the fact that the conservation of steepness al- 
lows one to relate the steepness of the initial condition 
to the steepness of the late stage evolution of the front 
that can be decomposed into an asymptotic front profile 
plus a linear perturbation. The spectral decomposition of 
this perturbation is largely determined by the steepness 
of the initial and the asymptotic state. 

Basins of attraction and rate of convergence are also 
discussed in Section II D. Flat initial conditions ( |2.8j ) 
approach a front characterized by their initial A. Suffi- 
ciently steep initial conditions fl2.7fl in the pushed regime 
(X se i > 1) evolve at late times into a pushed front cor- 
rected by linear perturbations that can be represented by 
eigenfunctions of the stability operator, whose spectrum 
has a gap. Hence the convergence of a pushed front is ex- 
ponential in time. In contrast, the rate of convergence of 
pulled fronts (X se i = 1) can not simply be obtained from 
the spectrum, as it is gapless, and generic perturbations 
are not spanned by the natural eigenfunctions. 

Leading edge and interior dominated dynamics. Both 
the stability analysis and our relaxation analysis bring 
out the importance of distinguishing leading edge domi- 
nated from interior dominated dynamics. The most ob- 
vious form of leading edge dominated d yna mics results 
from flat exponential initial conditions ( |2.8| ) with finite 
steepness A. In this case, the asymptotic front speed is 
just the speed 



v(X)=X+- 



(2.10) 



with which the exponential tail e Xx propagates accord- 
ing to the linear dynamical equation 



d t( p = d 2 x (j> + + o{(j) 2 ) 



(2.11) 



This equation is obtained by linearizing about the unsta- 
ble state (f> = 0, and is appropriate in the leading edge 
region. The more important leading edge dominated dy- 
namics occu rs, however, for sufficiently steep initial con- 
ditions (2.7) converging to a pulled front. As already 
mentioned, for pulled fronts the asymptotic front speed 
is just the linear spreading velocity v* determined in the 
leading edge where the dynamics is essentially governed 
by the linearized evolution equation. This type of leading 
edge dominated pulled dynamics occurs when the nonlin- 
earities in /(</>) are mostly saturating so that they slow 
down the growth. In passing we note that we rederive in 
Appendix the well-known sufficient criterion for pulling 
in the nonlinear diffusion equation, viz. 



f'(0) = max 

O<0<1 



(2.12) 



with the help of a transformation that we call the lead- 
ing edge transformation |85| , [76| ]. This form of a proof is 
generalizable to some other equations |7J| . Pulled fronts 
are actually at the margin of leadi ng e dge domination: 
although the linearized equation ( [2.1l| ) is sufficient to 
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determine v se i = v* = 2, we will see in Section [II that 
the convergence towards this velocity is governed by a 
nontrivial interplay of the dynamics in the leading edge 
and the "slaved" interior. 

Leading edge dominated dynamics contrasts with in- 
terior dominated dynamics, which occurs when the non- 
linear function /(0) is such that steep initial conditions 
give rise to pushed fronts. For interior dominated or 
pushed dynamics, v se i is associated with the existence 
of a strongly heteroclinic orbit in the phase space as- 
sociated with $ v (£) (Section |l B| ). This means that 
the whole nonlinearity f(4>) is needed for constructing 
v se i, not only the linearization /'(0) about the unst able 
state. The linear stability analysis of Section II C im- 
plies that pushed fronts converge e xpon entially in time 
to their asymptotic speed (Section [IE). This type of 
dynamics extends smoothly towards fronts p rop agating 
into metastable states, i.e., towards e < in (2.3). 

While in this Section we consider the nonlinear 
diffusion equation (2.1), ( |2.2| ) only, the straightfor- 
ward extension to generalized p.d.e.'s of the form 
F((j), d x 4>, d%4>, d t (j>) = can be found in Appendix |B]. 

In the following subsections the above assertions are 
further substantiated. Readers familiar with most of the 
concepts and results listed above can proceed to Section 
III. 



B. Uniformly translating fronts: 
candidates for attractors and transients 



In this Section, we recall some well known properties 
@^,||j6||6|j7l|,||@ of uniformly translating front 
solutions of the nonlinear diffusion equation (2.1), (2.2). 
We transform to a coordinate system moving with uni- 
form velocity v. (x,t) — >■ (£,£), £ = x — vt, so that the 
temporal derivative transforms as dt\ x — dt\^ — vd^ 1 t . 
For a front <p(x — vt) — $ v (£) translating uniformly with 
velocity v, the time derivative vanishes in the comoving 
frame 9t|^$„ = 0, and so $«(£) obeys the ordinary dif- 
ferential equation 



e%$ v + vd c $ v + /(*<,) = o 



(2.13) 



In view of the initial condition ( J2.5D , we will throughout 
this paper focus on the right-moving front and hence we 
impose the boundary conditions 



$„(£)->! fore 
$„(£)-► fore 



-co 



(2.14) 



Close to the stable state = 1, the differential equa- 
tion can be linearized about 0=1 and solved explicitely. 
The general local solution is a linear combination of 
e~ A± £ with 



According to ( |2.2| ), /'(l) is negative. Thus for any real v, 
X + is positive and A_ is negative. With the convention 
( |2.14 ), only the negative root is acceptable. So 



$„(£) =l±e 



o(e 



for e — * — oo 



(2.16) 



A+ = 



■±y/v*-4f'(l) 



(2.15) 



The free integration constant multiplying e~ A ~^ here has 
been decomposed into a sign ± and a free parameter eo 
accounting for translation invariance. Apart from trans- 
lation invariance, there are two solution for <E>„ close to 
0=1 distinguished by ±. 

A global view of the nature and multiplicity of solu- 
tions can be obtained with a well known simple particle- 
in-a-potential analogy. This analogy has of course been 
exp loited q uite often in various types of approaches 
p3,100 "l02| , and only works for the nonlinear diffusion 
equation, not for equations with higher spatial deriva- 
tives; for these, we have to rely on a construction of so- 
lutio ns as trajectories in phase space as sketched around 
Eq. fl2~24[ ). 

The particle-in-a-potential analogy is based on the 
identification of equation (2.13) with the equation of mo- 
tion of a classical particle with friction in a potential. 
One identifies <&„ with a spatial coordinate, £ with time, 
v with a friction coefficient, and / with the negative 
force, / = —force = d l f ) V{4)) derived from the poten- 
tial V{4>) — f d</)' f {<!>')■ The potential has a maximum 
at = 1 and a minimum at = 0. The construction 
of is equivalent to the motion of a classical particle 
with "friction" v in this potential, where at "time" — oo 
the particle is at rest at the maximum of V. Obviously 
for any positive "friction" v > 0, the particle will never 
reach the minimum at = 0, if it takes off from the 
maximum at = 1 towards > 1. It will always reach 
= if it takes off towards < 1. Thus for every v > 0, 
there is a unique uniformly tran slatin g front (unique up 
to a translation), that starts as ( 2.14 ) and reaches = 
monoto nical ly. Close to = 1 it is given by the — branch 
in Eq. (§t|). 

Let us be more specific on how = is approached. If 
the "friction" v is sufficiently large, the motion of the par- 
ticle will be overdamped when it first approaches = 0, 
it will reach = only for "time" £ — > oo, and form 
a monotonic front over the whole £ axis. This behavior 
continues down to a critical value of the "friction" v c . It 
defines the critical velocity v c as the smallest velocity at 
which monotonically reach es - ► at £ — * oo. 

(As we will discuss in Sect. [I C , a uniformly translating 
front is dynamically stable if and only if v > v c .) If 
v < v c , the particle will reach = at a finite "time" 
e and cross it. What then happens, depends on /(0) for 
negative arguments. If /'(0) = 1 for both positive and 
nega tive argu ments as in the case of the nonlinearities 
(1.2) or ( l.lOD , the particle might oscillate a finite or an 
infinite number of times through = and reach = 
asymptotically for £ — > oo as 
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where 



A v e^ A -« + B v e~ A +« 
K + /3) e- A *« 



for u > 2 , 
for v — v* - 
for \v\ < 2 



<M0 



(2.22) 



(2.17) 



A±(«)=Ao(«)±/i(») («>2), A (w) = - (aU«), 

(2-18) 

(2.19) 

A*=A (v*) = \±(v*) = 1 (u = v*=2). (2.20) 



The solution ( 2.17 ) of the equation linearized about = 
contains two free parameters for every v. These parame- 
ters are determined by the unique approach of the front 
<5> v from <f> — 1 and will in general both be non- vanishing. 
The special value v* = 2 is determined by linearization 
about the unstable state. As can be seen from (2.17), 
it is a lower bound on the critical velocity v c . At this 
value of the velocity the two roots A + and A_ coincide. 
As a result, the asymptotic profile is not the sum of two 
exponentials, but an exponential times a first order poly- 
nomial in £. 

Depending on the nonlincarity /, the critical v c can be 
determined by two different mechanisms that turn out to 
distinguish pushed (v c > v*) or pulled (v c = v*) fronts. 
Suppose first that upon lowering v the front solutions $„ 
remain monotonic till v = v* . In this case v c = v* is de- 
termined by the equation linearized about the unstable 
state, and we will see, that sufficiently steep initial condi- 
tions (2/7) evolve into pulled fronts. A second possibility 
is the following. At very large v, the front solution is cer- 
tainly monotonic, since in the particle-on-the-hill analogy 
the particle slowly creeps to the minim um of the poten- 
tial for large "friction" v. Hence A v in ( 2.17 ) is positive 
for large v. Now, depending on the nonlinearities, it may 
happen upon lowering v that at some velocity v = v\ 
A u t = 0. The front is nonmonotonic for v < as A v 



will be negative for v < v> . Hence in this case 



= «t 

and pushed fronts result. The pushed velocity ir thus 
emerges from the global analysis of the whole nonlinear 
front, and not only from linearization about the unstable 
state. 

For uniformly translating pulled fronts we will use the 
short hand notation $* = For large £ they are 

asymptotically 



since in general the coefficient a in (2.17) is nonzero. This 
particular form will in Section [II turn out to have impor- 
tant consequences for the convergence of pulled fronts: it 
determines the prefactor of the 1 ft relaxation term. 

For fronts with velocity v > v* , the smaller A will dom- 
inate the large £ asymptotics, so generically 



However, for a front solution with velocity v<, we have 
A,jt = 0, and so 



(2.23) 



An alternative formulation that can be generalized to 
higher order equations is the following. A construction of 
front solutions of Eq. (2.13) is equivalent to a construc- 
tion of trajectories in a phase space = d^ v ) in 
which the flow is given by 



dt 



-v9 v - /(*„) 



(2.24) 



Front solutions correspond to trajectories between the 
fixed points = (1,0) and (0,0). These are thus 

heteroclinic orbits in phase space. Out of the (1,0) fixed 
point come two trajectories i n opp osite directions along 
one eigenvector according to ( 2Tq ) . When we follow the 
direction for which decreases for increasi ng £, its be- 
havior near the (0, 0) fixed point is given by (2.17). Now, 
since t he flo w depends continuously on v, so will A v and 
B v in (2.17). For large v A v is positive, and from the con- 
struction of the flow in phase space one sees that A v may 
change sign on lowering v. The largest v with A v = 
determines the change from monotonic to non-monotonic 
fronts. At this v = v\ the trajectory flows into the sta- 
ble (0, 0) fixed point along the most strongly contracting 
eigen direction — this is precisely what is expressed in 
(IH). For this reason, the solution is referred to 
by Powell et al. |6{| as a strongly heteroclinic orbit. In 
IpSf , this solution was referred to as "the nonlinear front 
solution" . 

In summary, the main results of the preceding analysis 
are: 

• For every v > v c , there is a uniformly translating front 

with velocity v, which monotonically connects <fi = 1 
at £ — > — oo to (j> = at £ — > oo. All with these 
properties are uniquely determined by v up to transla- 
tion invariance. 

• For every < v < v c , there is a unique front solution 

that translates uniformly with velocity v, and that 
monotonically connects cf> = 1 at £ — > — oo to cf> = at 
some finite £ = £. 

• Depending on the nonlinearities, the change from 
monotonic to nonmonotonic behavior can either occur 
at the velocity v* , with v* = 2 for (2.1) and (2.2), or at a 



(2 21) larger velocity v^: v c — max[v* ,v<]. If exists, it is the 



largest velocity at which there is a strongly heteroclinic 
orbit. 

The results for invasion into either metastable (/'(0) < 
0) or unstable states (/'(0) > 0) and for v c = > v* 
and v c — v* are summarized in v(X) plots in Fig. 3, 
which show the multiplicity of stable uniformly trans- 
lating fronts parametrized by either v or A. 
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The results of this subsection play a role in the subse- 
quent analysis: 

(i) There are important connections p3| between the 
properties of the uniformly translating front solu tions 
and the stability of these fronts (see Section II C). In 
particular front solutions with velocity v > v c are dy- 
namically stable and possible attractors of the long time 
dynamics. Fronts with velocity v < v c either do not exist 
or are unstable. 

(ii) The results for front selection can be easily for- 
mulated in terms of th e prop erties of these uniformly 
translating solutions Q,^5 67 1: for sufficiently steep ini- 
tial conditions the dynamically selected velocity coincides 
with v c : v se i — v c . If v se i = d*, we speak of the pushed 
regime, while if v se i — v* we speak of pulled fronts. 



(Hi) We will see in Section [II, that the positive mono- 
tonic part of the front solutions $> v (0 with velocity 
v < v* plays a role in the convergence behavior in the 
interior region of pulled fronts. Note, however, t hat w hile 
a solution of the o.d.e. has according to (2.17) an 

oscillatory leading edge for large £, that causes the dy- 
namic instability of these solutions, the relaxing front 
is approximated by & v (t) only in the interior front re- 
gion and crosses over to a different functional form in 
the leading edge. This behavior is in agreement with the 
conservation of positivity of the solution in a nonlinear 
diffusion equation, if the initial condition was positive. 
All arguments essentially also apply to higher order equa- 
tions, though then the positivity and monotonicity prop- 
erties of the solutions loose their distinguished role. 




FIG. 3. Steepness A (2.£) versus velocity v(X) = A + /'(0)/A with solid line for real A and dotted line for real part of complex 
A. is the pushed velocity derived from global analysis, v* the linear spreading velocity. A fat line or point on the axes denotes 
the possible attractors <& v (x — vt) of the dynamics, parametrized either by velocity v or by steepness A. a) The case /'(0) < 
corresponding to front propagation into a (meta)stable state. In this case, there is a unique attractor with velocity v ae i = 
and steepness A se i = A+(t)'). b) and c) The case /'(0) > corresponding to front propagation into an unstable state. In this 



case there is a continuum of attractors parametrized by v > v c . b) The pushed regime: u se ; 



> v*. The steepness 



X se i = A+(v T ) of the steepest attractor is isolated just as in case a. There is a continuous family of fronts parametrized by 
< A < \ a teep = A_(tr). c) The pulled regime: v se i = v c — v* . The steepness A* = X 3 teep = A se ; of the steepest attractor is 
at the margin of the A-continuum of attractors. 



C. Linear stability analysis of moving front solutions 

To study the linear stability of a uniformly translating 
front <!>„, we linearize about it in the frame £ — x — vt 
moving with the constant velocity v, by writing 



(2.25) 



Inserting (2.25) into ( J2.l[ ), we find to linear order the 
equation of motion for ?y(£, t) 

d t r) = £ v r) + 0{r) 2 ) 

with the linear operator 



(2.26) 



(2.27) 



C v is not self-adjoint, so left and right cigenfunctions 
will differ. The trouble is caused by the linear derivative 



v dc . It can be removed by the following transformation 

IIS 



e<' 2 r, 



e</ 2 C v e"</ 2 . 
H. v is the linear Schrodinger operator 



(2.28) 
(2.29) 



-dl + V{i) , V{t) = £ -/'($.„ (0) , (2.30) 



and the equation of motion (2.26) transforms to 

- d t ip = H V *P + O (V e-"« /2 ) . (2.31) 

This Schrodinger problem is, of course, well known to 
physicists Q (see also p8| , p0 94,5^]), as long as ip lies 
in the natural Hilbert space of TL V . 
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However, the transformation ( |2.28 ), (2.29) increases 
the weight of the leading edge (£ — > oo) by a factor e 1 ^ 2 , 
while it enhances convergence at £ — > — oo. Therefore, 
only perturbations with 



lim 



with \q(v) 



(2.32) 



are spanned by the eigenfunctions within the conven- 
tional Hilber t space of 7i v . For the selection analysis in 
Section IIP below, this function space in general is not 
sufficient. As is discussed in detail in Appendix D, one 
can constru ct eig enmodes of C v outside the Hilbert space 
defined by ( 2.32Q . With this extension of the function 
space, scalar products of arbitrary eigenfunctions might 
be divergent, so one looses the efficient tool of projection 
onto eigenfunctions by taking inner products. Neverthe- 
less, in most cases generic perturbations still can be de- 
composed into these eigenfunctions, except in the case of 
pulled fronts: the linear perturba tion 77 of a sufficiently 
steep front <f> ~ e _A ^ with A > 1 ( |2.7|) abou t t he as ymp- 
totic pulled front - « + /3) (T% ( fT?] ), (^2l|) will 
decay asymptotically as 



n 



-K + /3) e 



(2.33) 



Since there is only one zero mode of translation with a 
slightly different asymptotic behavior 



£ v r]o = 



Vo 



-(of + j8- a) e~« , (2.34) 



the asymptotics of 77 cannot fully be decomposed into 
eigenfunctions. Here the double root structure of the 
leading edge with a ^ plays a crucial role, as it later 
will do again. 

The most important conclusions from the present dis- 
cussions and the detailed Appendix [5] are: 

1) Non-monotonic fronts are intrinsically unstable, and 
generically will not be approached by any initial condi- 
tion. 

2) Monotonic fronts propagating with velocity v are 
stable against perturbations steeper than e~ x -^ v ^. 

3) Perturbations 77 about pushed fronts , that decay 
more ra pidly than A-('ir) have a gapped spectrum 
see Eq. (D14). The same holds for perturbations about 
fronts with a velocity v > v\ if their steepness is 
larger than A_(u). 

4) The spectrum of pulled fronts is gapless and cannot 
be decomposed into eigenfunctions of 7i v even outside 
the conventional Hilbert space. 

Before closing this subsection, we note that although 
the particle-on-a-hill analogy for or the mapping onto 
the Schrodinger equation for r\ a are insightful and very 
efficient ways to arrive at our results for existence and sta- 
bility of uniformly translating front solutions, the analy- 
sis by no means relies on these. In fact, much of the phase 
space analysis can easily be generalized to higher order 
equations as those shown in Table I. E.g., in the stability 



analysis of non-monotonic fronts, the discrete set of solu- 
tions with A v = plays a particular role. For equations 
like the EFK equation from Table I, monotonicity ceases 
to be a criterium, but conditions like A v = defining 
so-called strongly hetcroclinic solutions continue to play 
a central role in the stability analysis, as is discussed in 
Appendix |] to Sect. [v|. 



D. Consequences of the stability analysis for 
selection and rate of convergence; marginal stability 

Suppose now, that we start with an initi al c ondition 
4>(x,0) in the nonlinear diffusion equation (2.1) with a 
given nonlinearity /(</>), and then study the ensuing dy- 
namics. What will the linear stability analysis tell us 
about the asymptotic (t — > 00) state and the rate of con- 
vergence? It turns out that the issue of selection is more 
closely related to that of stability than one might expect 
at first sight. The reason is the co nserv ation of steepness 
discussed in more detail in Section II E below: If initially 
at t = the steepness A defined in (2.t) is nonzero (finite 



or infinite ), then at any finite time t < 00 the steepness 
is conserved: 



<t>(x,0) 



4>{x,t) 



for all t < 



(2.35) 



Note that the limits x — + 00 and t — + 00 do not commute. 
We characterize the initial condition by its steepness Xinit 
defined by 



(x,t = 0) 



(2.36) 



As a consequence of ( 2.35 ), we can use Xinit to character- 
ize not only the initial conditions but also the profile at 
any later time < t < 00, when the front velocity might 
be already close to its asymptot ic val ue. 

The conservation of steepness ( 2.35 ) entails that a front 
characterized by an initial steepness Xi n u, will be char- 
acterized by the same steepness after any finite time, so 
also at a late stage when the velocity and shape of a front 
are close to their asymptotic limits. At such a late stage, 
the front <fi can be decomposed into a possible attractor 

( x — v t) of the dynamics plus a linear perturbation r] as 
in ( |2.25| ). Characterize the attractor that we investi- 
gate by its steepness X aS ympt ■ The resulting perturbation 
r](x,t) = (f>(x,t) — & v (x — vt) then will have steepness 

(2.37) 



A r; = min [A 



init 1 X aS y m pt\ 



Whether the perturbation 77 will grow or decay, that 
means, whether with a particular velocity v is the at- 
tractor of the evolution of <f> or not, is determined by the 
decomposition of the perturbation 77 into eigenmodes of 
the linear operator. Whether this spectrum has growing 
eigenmodes, depends on the operator and the functions 
space defined by the steepness A^. With the tools of the 
stability analysis from Appendix the selection ques- 
tion can therefore be rephrased purely in terms of A^, 
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A. 



asympt 



and the two steepnesses A steep and A se ; charac- 



terizing the nonlinearity /: for pushed fronts 



A 



steep 



A_(« f ) , A seJ =A+(» + ) , v sel =J , (2.38) 



and for pulled fronts 



A 



— A se ; — A* 



V se l 



(2.39) 



in the notation of ( |2.18D - (|2.20|) . 

A detailed discussion of the question to what extent 
one can understand the selection and rate of conver- 
gence of fronts following this line of analysis, is given 
in Appendix || and summarized in Table IV below. The 
two most important conclusions for our purposes concern 
fronts evolving from sufficiently steep initial conditions: 

(i) The gapped spectrum of the conventional Hilbert 
space for pushed fronts implies that the relaxation to- 
wards pushed front solutions is exponential in time. 

(ii) Even after extending linear stability analysis be- 
yond the Hilbert space, it is not possible to derive the 
rate of convergence of pulled fronts from the stability spec- 
trum, since it is gapless, and generic perturbations can- 
not be decomposed into eigenmodes of the linear stability 
operator even in an enlarged functions space. 

We finally note that the usual marginal stability view- 
point is to characterize the family of stable front solu- 
tions <&„ by the velocity v; from this perspective, the 
front velocity v se i selected by the sufficiently steep ini- 
tial conditions is at the edge of a continuous spectrum 
of stable solutions with v > v se i ■ In this sense, both the 
pushed and the pulled attractors are marginally stable 
p3| , |65|j67| . The picture changes, however, when the at- 
tractors are not characterized by the velocity v, but by 
their asymptotic steepness A, see Fig. 3. The pulled front 
then still is at the margin of a continuous spectrum, while 
the pushed front is isolated just like the bistable front. 



E. The dynamics of the leading edge of a front 

In this Section, we reconsider the dynamics in the lead- 
ing edge in more detail, first to d emon strate the conser- 
vation of steepness expressed by (2.35), second to clarify 
the dynamics that ensues from flat initial conditions, and 
third to lay the basis for the quantitative analysis of the 



relaxation of pulled fronts in Section HI. 



1. Equation linearized about <j) = 

When we analyze the leading edge region of the front, 
where \cf>\ <C 1, we to lowest order can neglect o((j> 2 ) in 
( [2.11 ) and analyze 



d t cf> = dl 



(2.40) 



We first explore the predictions of this equation, before 
explori ng the corrections due to the nonlinearity / in 



Section HE 2 



Eq. (2.40) is a linear equation, so the superposition of 
solutions again is a solution. A generic solution is, e.g., 
an exponential e~ Xx . It will conse rve s hape and propa- 
gate with velocity v(X) = A + 1/A (|2.10[ ): 



<j)(x,t) 



Mm 



(2.41) 



The minimum of v(X) is given by v* — v(X* = 1) = 2. 

Consider now a superposition of two exponentials 
ci e~ XlX + C2 e~ X2 ' x . Without loss of generality, 
we can assume the maximum velocity to be v max = 
max[w(Ai), f(A2)] = In the coordinate system 

£i = x — v(\\)t the temporal evolution then becomes 



(x,t) = ci e~ AlCl + c 2 e 



o = \ 2 (»(Ai) -«(A a )) > 



(2.42) 
(2.43) 



Clearly, the contribution of A2 decays on the time scale 
l/<7, and so for large times ^> 1/cr, the velocity of a so- 
called level curve of (j) = const. > in an a;, t diagram 
will approach v(Xi) and the profile will converge to e~ Al ^ 
(sec for a similar type of analysis) . The steepness of 
the leading edge at £ — > 00, on the other hand, will be 
given by A m in = min[Ai, A2] for all times t < 00. 

This simple example already backs up much of our 
discussion of perturbations outside the Hilbert space in 
Appendix |], that apply to the Cases II and III in Table 
IV: 

1) The limits £ — > 00 and t — > 00 in general do not com- 
mute. 

2) The steepness A = minj[Aj] is a conserved quantity 
at x — > 00 and t < 00. As the explicit example of [10E] 
shows, for equations for which one can derive a compari- 
son theorem, the conservation of steepness can easily be 
derived rigorously. 

3) The velocity of a constant amplitude (f> = const. > 
will be governed by the quickest mode present v — 
maxj[z;(Aj)] at large times t > 1. 

Let us now analyze initial conditions steeper than any 
exponential. Q uite generally, an initial condition <f>(x, 0) 
evolves under ( 2.40 ) as 3 



(x,t) 



dy cf)(y,0) 



(x-y) 2 -4t 2 /(4t) 



V^t 



(2.44) 



3 Eliminate the linear growth term in (2.4C) by the transfor- 
mation 4> = e*<^, solve the diffusion equation dt<j> = d^.(f>, and 
transform back. 
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Assume for simplicity, that the initial condition 4>(y, 0) 
is strongly peaked about y = 0, so that for large times, 
we can neglect the spatial extent of the region where 
4>(y, 0) =/= initially. Upon introducing the coordinate 
£ = x — 2t we get 



e -?-? 2 /(4t) 

(x, t) oc ■= for t >• 1 



(2.45) 



This general expression leads to three important obser- 
vations: 

1) The steepness of the leading edge characterized by 
A = oo at £ — > oo indeed is conserved for all finite times 
t < oo. 

2) At finite amplitudes <j> — const. > and large times 
t, the steepness of the front propagating towards £ — > oo 
approa ches A* = 1 and the velocity approaches v* = 2. 

3) Eq. (2.45) furthermore implies, that a steep initial con- 
dition like <f>(y, 0) approaches the asymptotic velocity v* 
as 



v{t)i 



0, 



1 ( 1 
2 -2t +0 [¥ 



(2.46) 



where we defined the position £h{t) of the amplitude h 
in the comoving frame £ = x — 2t as 0(Oi(t),t) = h. 
Eq. (2.46) is then obtained simply by solving \n<fi — 
-~%fj±t - (lni)/2 = const 
This algebraic convergence is consistent with the gap- 
less spectrum of linear perturbations, and as such it iden- 
tifies the missing link in the analysis of the relaxation of 
pulled fronts. However, Bramson's work j77j shows, that 
the qualitative prediction of convergence as 1/t is right, 
but the coefficient of 1/t is wrong. In fact, the mathe- 
matical literature |54j has established ( 2.46j ) as an upper 
bound for the velocity of a pulled front in a nonlinear dif- 
fusion equation. The algebraic convergence clearly comes 
from the 1 / \ft prefactor characteristic of the fundamen- 
tal Gaussian solution of the diffusion equation ( 2.45 ) - 
this qualitative mechanism will be found to be right in 
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We fini sh ou r discussion of solutions of the linearized 
equation ( [2.40 ) with another ill ustra tive example. After 
the discussion of the solution ( [2.41 ) one might be wor- 
ried about initial conditions with A> 1. Such an initial 
condition is steep according to our definition, so it should 
approach the velocity v* . But according to ( 2.41 ), it ap- 
proaches the larger velocity w(A). However, even in the 
framework of the linearized equation, this paradox can 
be resolved: An initial condition e~ Xx on the whole real 
axis is, of course, unphysical, and we in fact only want 
this behavior at x 3> 1, where <j> is small. Let us there- 
fore truncate the exponential for small x by writing, e.g., 
4>(x, 0) = 6(x) e~ Xx , with 9 the step function. Insertion 



into (2.44) yields the evolution 



-v(\)t] 



1 + erf 



:-2\t 
V4t 



where erf x = 27r~ 1 / 2 J* dt e~ l is the errorfunction. For 
t ^> 1 the crossover region where x w 2Xt separates two 
different asymptotic types of behavior: 



( e -\[x-v(\)t\ 

</>(x,t) W < e -(x-2t)-(x-2t) 2 /4t 

[ V4irt \(l-x/(2\t)) 



for x > 2Xt 

for x < 2Xt ^■ 48 ' 1 



In the region of x 3> 2Xt we find our previous solution 
fl2.4l| ) with conserved leading edge steepness and veloc- 
ity v(X), whi le in the region of x <C 2Xt we essentially 
recover Q2.45 ), with £ = x — 2t. 

Considering the three different velocities — v(X) for 
the region of conserved A, v* = 2 for the "Gaussian" re- 
gion behind, and 2A for the crossover region between the 
two asymptotes — the distinction between flat and steep 
initial conditions now comes about quite naturally: 

a) For flat initial conditions, we have A < 1, and an or- 
dening of velocities as 2A < v* < v(X). The crossover 
region then moves slower than both asymptotic regions, 

so for large times the region of finite <f> will be dominated 

by e -A[a:-«(A)i]_ 

b) For steep initial conditions, we have A > 1, and the 
velocities order as v* < v(X) < 2X. The crossover region 
then will move quicker than both asymptotic regions, and 
the region of finite <f> will be dominated by e~^~^ ^ 4t /\/t, 
where £ = x — 2t. 

We finally note that the above results can also be rein- 
terpreted in terms of the intuitive picture advocated in 
H0: The group velocity v gr (X) — dv(X)/dX of a near 
expon ential profile in the leading edge is, according to 
(ETO), negative for A < 1 and positive for A > 1. In this 
way of thinking, the region with steepness A in the case 
considered above expands when A < 1 since the crossover 
region moves back in the comoving frame [case (a)] , and 
it moves out of sight towards £ — > oo for A > 1 [case (&)], 
since the crossover region moves faster than the local co- 
moving frame. 



2. Leading edge representation of the full equation 

Just as the linear stability analysis of the front was 
insufficient to cover the full dyn amical behavior of the 
nonlinear diffusion equation (2.1) and in particular the 
dynamics of the leading edge, so is the linearized equa- 
tion (fOcj). In Section || we will see that only through 
joining these complementary approaches, we can gain a 
quantitative understanding of the convergence of steep 
initial conditions towards a pulled front <&*. 

The shortcomings of the linearized equation ( 2.40[ ) be- 
come quite clear by confronting it with what we will call 
the leading edge representation of the full equation (JO): 



cW = df</> + /(</>,£) 



where we transformed with 



(2.47) 



ip = <f) e 



t = x-v*t 



(2.49) 



(2.50) 
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The parameters are A* = v* /2 = y/ f'(0). This transfor- 
mation eliminates the terms of order ip and d^ip from the 
linear part of the equation. The nonlinearity is 



MO=e yi (/(V> 
= O (ip 2 e 



- /(0) - /'(0) V e 



s 2 



(2.51) 



This transformation is quite comparable to the transfor- 
mation of a linear perturbation r\ into the Schrodinger 
picture as in (2.27), (2.28). However, we here transform 
the full nonlinear equation, and not only the linearization 
about some asymptotic solution. 

For, e.g., /(</>) = (j) — <fi 3 we have / = — tp 3 e~ 2 ^. When 



we neglect / in (2.4E), th e equ ation is equivalent to the 
linearization about <f> = (2.4C). The linearization is cor- 
rect for £ 3> 1, but the presence of the crossover towards 
a different behavior for smaller £ has important conse- 
quences for the solutions of the full nonlinear diffusion 
equation. 

In particular, for the leading edge of a pulled front 
$* ~ (a£ + (3) (2.17), we generical ly fin d a^O and 
accordingly the leading edge behavior (2.21). This lead- 



ing edge b ehav ior will play a central role in Section III 



In Section [I B we derived a ^ from the uniqueness of 
the trajectory in phase space, i.e., from the construction 
of the whole front from <f> = up to <j> = 1. We now 
will give a different argument for a^O from the analysis 
of (2.49), that does not rely on constructing the whole 
solution up to <j) — 1. 

The front propagates uniformly with velocity v* = 
2, so in the frame £ = x — 2t it is stationary. , 3>* = $*e^ 
then solves 

df#* + /(**, £) =0 . (2.52) 
The boundary conditions ( |2.14 ) for $* imply for \I>*: 



a£ + (3 for £ — > oo 
0(e x * 6 ) for£^-oo 



(2.53) 



The solution = a£ + (3 for £ — > oo can directly be 
derived from ( |2.52 ) and the condition, that $* vanishes 



at £ — > oo. Now integrate (2.52) over the real £ axis, and 
find 



(2.54) 



The integral on the right hand side is well-defined, since / 
vanishes exponentially, both for £ — * — oo and for £ — > oo. 
Clearly, a nonlinearity / ^ generically implies a ^ 
and hence the leading edge behavior ( 2.21[ ) for a pulled 
$* front. Only for particular nonlinearities /, we occa- 
sionally find a = (see (3.67) and Appendix]^). Having 
a = is obviously only possible if / has terms of op- 
posite sign, so that its spatial average vanishes. Fo r the 
nonlinearity of form / = <\> — <fi k with k > 1 (1.2), we 
find a always, and in this case the term / act s like a 
localized sink term in the diffusion equation ( 2.52Q for ijj. 
This interpretation is especially useful for the discussion 
of the non-uniformly translating fronts [|74],f7Gl . 



F. Concluding remarks — interior and edge 
dominated dynamics 

Table IV summarizes the results of this Section for 
various nonlinearities and initial conditions. For pushed 
fronts, the stability analysis gives essentially all the in- 
gredients to determine the rate of convergence for pushed 
fronts. For pulled fronts emerging from sufficiently steep 
initial conditions one finds however (see Case IV of Ap- 
pendix that the linear stability analysis is not the 
appropriate tool: the spectrum is gapless and the rate of 
convergence can not be determined from the spectrum. 

The crucial insight for the further analysis is that a re- 
laxing front can be decomposed into different dynamical 
regions. Linear stability analysis is the appropriate tool 
for the interior dominated dynamics of a pushed front. 
Pulled fronts and fronts evolving from "flat" initial con- 
ditions are leading edge dominated. This calls for differ- 
ent methods of analysis. The relaxation of pulled fronts 
emerging from sufficiently steep initial conditions will be 
addressed in the next section. 
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INITIAL CONDITION 

00,0) ~ e~ Xin " x 
as x — > oo 


7 

mctastable (/'(0) < 0): 
w se ; = ut > 


\T ONLINE A RITY f((f>) 

unstable (/'(0) > 0), 
pushed regime: 

W se ; = W t > 

= 2v//'(0) > 


unstable (/' (0) > 0), 
pulled regime: 

Vsel = V* 

v* = 2V7'(0) > 


steep: 

Aimt > Ao(w S ei) 

(including A — oo) 


Case I- 
pushed dynamics, 

A — ► A+(«+) 
u(t) = + ©(e-' 74 ) 


Case I- 

V.. CIO v.y _L ■ 

pushed dynamics, 

A — ► A+(«+) 
u(t) = ^ + 0(e- CT *) 


Case IV- 
pulled dynamics, 

A — ► A* 
«(t) = v* + 0(l/t) 


steep: 


Case II: 
pushed dynamics, 

v — ► v^ 
generically: 
v(t) = «t + 0(e-' Tt ) 


Case II: 
pushed dynamics, 

A ► A_i_(Vh 

t> — > 
generically: 
= u f + 0(e~ CT *) 


not applicable, since 
A+Cu* 1 ) — Xn(v*) — A* 


flat: 

< Ai„it < A_(f ;je i) 


not applicable, since 
A_(«t) < 


Case III: 
leading edge 
dominated dynamics, 
A — » A < A_(«t) 
w — ► w(A) > 
generically: 
v(t) = «(A) + 0{e- at ) 


Case III: 
leading edge 
dominated dynamics, 
A — ► A < A* 
v — ► U (A) > v* 
generically: 
v(t) = v(X) + 0(e- at ) 



Table IV: Table of initial conditions and nonlinearities, resulti ng in relaxation cases I - IV from Appendix 
[5|. Fronts at all times t are characterized by their steepness A (2J3) in the leading edge, and an arrow — > 
indicates the evaluation of the quantity for t — > oo. The nonlinearity only ente rs through the existence of 
a strongly heteroclinic orbit Q v (x — vt) wit h v se i = v' > v* (see Section II B ) or its non-existence (then 
v se i — v*). Vsei determines X± l o(v se i) as in (2.18), which in turn classifies the initial conditions. Pushed or 
pulled dynamics are special cases of interior or leading edge dominated dynamics for steep initial conditions. 
Cases I - III are treated in Appendix ||] with stability analysis methods and generically show exponential 
relaxation. Case IV is not amenable to stability analysis methods. It shows algebraic relaxation and is 
treated from Section III on. 
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III. UNIVERSAL PULLED CONVERGENCE OF 
STEEP INITIAL CONDITIONS IN THE 
NONLINEAR DIFFUSION EQUATION 

In the present Section, we will combine our under- 
standing of the dynamics of the leading edge and of the 
interior of a front into one consistent analytical frame, 
that allows us to calculate the long time convergence of 
steep initial conditions towards a pulled front — as we 
discussed in the previous section, the relaxation can in 
this case not be obtained from the linear stability anal- 
ysis of the asymptotic solution. The different dynamical 
regions of such a front are sketched in Fig. 2. We match 
an expansion in the interior, that resembles features of 
the linear stability analysis, to an expansion of the lead- 
ing edge. Both expansions are asymptotic expansions 
in \j\fi. This approach allows us to derive the power 
law convergence of the front velocity and the front pro- 
file towards <&*. This convergence is universal in leading 
and subleading order and we calculate all universal con- 
vergence terms analytically. For clarity we present the 
detailed calculation for the nonlinear diffusion equation 
in this Section first, and then discuss the generalization 
m Section [\]. 



A. Observations which motivate our approach 

1. Asymptotic steepness of leading edge determines rate of 
convergence 

Our calculation of the spreading of t he le ading edge 
under the linearized equation in Section II E gave qual- 
itatively the right results, but failed to reproduce the 
quantitative results for the nonlinear equa tion: Insert- 
ing sufficiently ste ep in itial conditions (2.7) into the lin- 
earized equation (2.11), we found that the asymptotic 
shape (2.45) approaches times a Gaussian for t — > oo 
and £ 3> 1 and that this implies for the asymptotic con- 
vergence that v(t) Un = 2 - l/(2t) + . . . Q2.461 ). For the 
nonlinear equation we know that the asymptotic front 
profile $* behaves as $* ~ £ for £ » 1 (2.21) and 
Bramson has derived with probabilistic methods, that 
v(t) = 2 — 3/(2i) + . . . independent of the height at which 
the velocity is measured j77j . 

How the exact result of Bramson jnj comes out nat- 
urally and generally is brought out quite clearly by 
rephrasing an argument of |67]] as follows (see also 



Let us work in the leading edge representation ( 2.49 ), 
(2.50), and let us from here on use the co-moving vari- 
able £ specifically for the frame moving with the pulled 
velocity v* = 2, 



v t 



2t 



(3.1) 



The funda ment al similarity solution of the diffusion-type 
equation ( 2.49D for the leading edge variable ip — 4> eX ^ 



in the region where the nonlinearity can be neglected, is 
of course the Gaussian 



(3.2) 



It reproduces our solution ( 2.45 ) for cf>. But also any 
derivative of the Gaussian ip n = d^ipo solves (2.4E) for 
£ 1. The ipn/ipo are simply Hermite polynomials 



[ 103, 104 1 . In particular, the dipole solution 



-r/(4t) 



t 3/2 



(3.3) 



also solves the diffusion equation ( [2.4S| ) for £ ^> 1 and 
has the p rop er asymptotics <I> oc £ for t — > oo. Trans- 
forming (3.3) back to 0, we find 



<f>(x,t) cx (x-2t) e -[*-«+(3/2)M e -(x-2t)V(4t) _ (3 4) 

If we now trace the position 2t + Xf l (t) of the point where 
4> reaches the amplitude h in the ori gina l x frame, we find 
by solving (j>(2t + X h , t) = h from (O) for X h (t) < y/U 



v(t) 



X h 



3 
2i 



(3.5) 



in agreement with Bramson's result. This indicates that 
for large times t> 1 and far in the leading edge £ 1, 
the conver ging front is approximately given by ( |3.3| ) if 
a ^ in ( 2.17 ) — remember that a ^ implies tha t 
$*(£) ~ £e 4 for large £. We will see indeed that fOj ) 
does emerge as the dominant term in a systematic asymp- 
totic expansion in the leading edge region. For reasons 
explained below, it is, however, more convenient to for- 
mulate this expansion in a slightly different frame. 



2. Interior follows leading edge: uniform convergence 

The above argument shows that the leading — 3/(2t) 
velocity correction is due to the diffusion-type dynam- 
ics in the leading edge. Why would the convergence to 
the asymptotic profile be uniform, i.e., be independent 
of the level curve whose position is tracked (the level 
curve is the curve in an x-t plot that indicates the val- 
ues where 4> reaches a particular level or "height")? The 
answer to this question is intuitively quite simple. As 
Xh — —3/(2t), Xh — —3/2 Int. If we compare the posi- 
tion Xh ± of a height hi in the leading edge (hi <§; 1) with 
a position Xh 2 of a height in the interior (hi = 0(1)) 
where the dynamics of <f> is described by the nonlinear 
equation, we will have Xh 2 — Xh t — W(hi,h2), where 
W is the width of the front between these two heights. 
Clearly, if W approaches a finite value for long times, we 
need to have also Xh 2 — 3/2 In t in dominant order as 
t — > oo, and hence also X} l2 = — 3/(2i) + • • •. But an 
equation of motion like ( |1.1| ) has front solutions whose 
width is finite, so we expect indeed that W — 0(1) for 
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large times. Our analysis will confirm this expectation. 
In other words, the leading order velocity correction as 
— 3 / (2t) is set by the dynamics of the leading edge, and 
because of the finite asymptotic width of the front, the 
convergence is uniform, i.e., independent of h. 



3. Choose proper frame and subtraction for the interior 

The above observations have another important conse- 
quence. After the front has evolved for some time, we will 
find it selfconsistent to assume, that its shape will resem- 
ble the asymptotic shape . If we want to understand 
the interior part of the front, it might at first sight seem 
appropriate to linearize the converging front </> about the 
asymptotic front However, the profile $* propagates 
uniformly with velocity 2, while as we saw above, the 
transient profile <f> propagates with velocity v* — 3/(2t). 
Thus, if the interior regions of the <f>- and the <I>*-fronts 
are at about the same part of space at time to, their dis- 
tance will diverge as (3/2) ln(t/to) as t grows! This was 
already illustrated in Fig. 1. Hence, linearization of <p 
about the asymptotic profile 4>* during the whole time 
evolution requires to move <I>* along with the velocity 
2 — 3/(2i) + ... of and not with its proper velocity 2. 
Our expansion is therefore based on writing <fi as 



where 



u = e - x(t) 



2t - X(t) . 



(3.6) 



(3.7) 



This Ansatz anticipates that we need to shift the profile 
$* an appropriate distance X{t) oc lnt, and that with a 
proper choice of X(t), r\ becomes a small and decaying 
perturbation. 



4- Choose proper expansions and match 
interior 



ie to 



We will need two different expansions for the leading 
edge and for the interior. The expansions have to be 
chosen such that they can be matched in overlapping in- 
tervals through resummation of the e xpa nsions. 

Since we use the coordinate system ( |3.7|) in the interior, 
we also should use it in the leading edge. The leading 1/t 
contribution from the leading edge suggests to expand rj 
in the interior as rji(£ t x)/t + ■ ■ ., and we shall see indeed 
that such a form emerges automatically from the ansatz 



( |3.6[ ) . The appropriate variable for £x ^> Vi in the lead- 
ing edge, on the other hand, is the similarity variable of 
the diffusion equation 



f 2 
22L 

At 



(3.8) 



as suggested by (3^) - (3^). Expressing £x by z and 
t introduces a dependence on We find, that it 

is actually consistent to expand the interior in powers 
of \j\ft (instead of 1/t) times functions of £x, and the 
leading edge also in powers of l/y/t times functions of z. 

The structure of these expansions is essentially our 
only input. Given this structure, the leading and sublead- 
ing order universal terms of the expansions are uniquely 
determined. 



B. Expansion in the interior region 

We first analyze the interior part of the front where 
4> varies from close to to close to 1. We work in the 
comoving frame £x = x — v*t — X(t) of (|3.7|), where X 
will have to be determined. We expand </> about 4>*(£x) 
as in (|]^). Because of translation invariance, we have 
the freedom to fix the position of $* and the zero of the 
coordinate system by imposing 



mt) = \ and $*(0) = i 



V(0,t)=0. (3.9) 



For <!>*, one has <&*(— oo) = 1, and we also assume that 
4> approaches 1 for —> — oc. This results in the second 
condition on r) 



fx 



lim v (Z x ,t) = 



(3.10) 



We insert <f> into the equatio n of motion (1.1), trans- 
form x to the coordinate ^x (3.7) and find for rj the 
equation 4 



d t rj = <9f 77 + v*d i n + Xd^* + /($* 



V) -/(**) • 

(3.11) 



Once 77 is small enough because time has evolved suf- 
ficiently long, /($* + 77) can be expanded in rj and we 
find 



d t r) = C*r/ + X + X d^rj 



/«($*) 



v 2 + o( v 3 ) , 

(3.12) 



4 Throughout this paper, we shall suppress the index X 
on partial derivatives with respect to £x for notational con- 
venience. Since d^ x \ t = d^\ t , this does not lead to any 
ambiguities. 
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where 



£* = dl + v*d e + /'($* (U)) 



(3.13) 



is the lineariz ation oper ator ( 2/27) for v = u*. 

In Sections [I E 1 and [II A 1 we have argued, th at one 
expects X(t) — 0(t _1 ). Asymptotic balancing in ( 3.12j ) 
then requires, that the leading order term of r] is of the 
same order r\ = 0(t ). We therefore try to expand as 
V = r li{£x)/t + ■ ■ ■■ We have argued, that connecting the 
interior expansion to the leading edge expansion requires 
an ordening in powers of 1/ y/t . So we choose the ansatz 



X 



c l | C 3/2 C 2 
t i3/ 2 t 2 + ' 



(3.14) 
(3.15) 



Substitution of the above expansions into ( 3.15 ) and or- 
dering in powers of yields a hierarchy of o.d.e.'s of 
second order: 

£*Vi = -ci<9 e $* , (3.16) 
C* n = -c^* , (3.17) 

C*m = -Vi - cid 6Vl - c 2 a 5 $* - /"(S>?/2 , (3.18) 
-/"(^*)'7i»?| etc., generally: (3.19) 

n-2 



n-2 



V 



2? ~ E C ¥ 9 ?^ 



m=2 



fc=2 



(3.20) 



It is important to realize that we do not need to drop 
nonlinear terms, but that the expansion of /($* + rf) 
in powers of rj is also ordered in powers of l/yt. So the 
higher order terms n n find their natural place as inhomo- 
geneities in the equations for r]i for i > 2. The hierarchy 
of o.d.e.'s is such that the differential equation for r\i con- 
tains inhomogeneities that depend only on r/j with j < i. 
The equations therefore can be solved successively. Each 
r}i solves a second order differential equation, and the two 
cons tants of int egration are fixed by the two conditions 
@ and (|3~T0|) 5 . 

Note also, that the time dependent collec tive c oordi - 
nate X(t) in £ x = £-X(t) only enters Eqs. ( |3.16D -(^2C|) 
in the form of the constants c„/ 2 , which at this point 
are still undetermined, and that the functions ry„/ 2 obey 
o.d.e.'s. 

Let us now compare 77 = <f> — $* to the variations of 
the profile shape with velocity v, 



5 = $ 



v*+X 



- $* = X T] sh 



X 2 



(2) 



(3.21) 



where f? s h = 8& v /5v\ v * is a "shape mode", which gives 
the change in the profile under a change in v. By consid- 
ering variations of v in the o.d.e. for the profile <fr v , we 

(2) 

find that ?/ s h and 77V obey 

C* Vsh + d^* =0 , 
C* V H ] +2d^ sh + /"(**) feh) 2 = . 



(3.22) 
(3.23) 



Upon comparing ( |3~22|) - ( |3~2^ ) with (§T|) - 
we can identify 



V3/2 = c 3 / 2 ?? s h , 

C? (21 

V2 = C 2 T] sh + —rfj + Cip , 



(3.24) 
(3.25) 

(3.26) 



with p a correction term, that solves the equation 

C*p + r] sh =0. (3.27) 

In these differential equations, 77 s h , % h and p obey the 
conditions 

W0) = 0, 7^(0) = 0, P(0) = 0, (3.28) 
!fch(-cx>) = , ?7s ( h (-°°) = . P(~«>) = . ( 3 - 29 ) 



cf. (|3_9|) and fl3.10[ ). 

p is the first nonvanishing term that indicates the dif- 
ference between the transient profile 4>(x, t) and the uni- 
formly translating front solution with the instantaneous 
velocity 



v{t) = v* + X . 
as resummation of (f> yields 



Cl 



(3.30) 



(3.31) 



This equation confirms that up to order 1/t 2 , the profile 

shape is given by the solution $„(t) of the o.d.e. with the 

instantaneous velocity v(t). 

Some remarks on these results are in place: 

1) We see, that the dynamics in the front interior is 

slaved to the evolution of v(t) i mposed by the leading 

edge, as we anticipated in Section [II A 3 . 



5 Had we introduced an 771/2, we would have found the equa- 
tion £*??i/2 = with the unique solution ni/2 = 0. 
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2) The fact that the profile <f> is up to order t 2 given are 



by 3>v(t){£x) can be traced back to the fact that since 
v(t) varies as t -1 , the time derivative term dt4> in the 
dynamical equation generates terms of order t~ 2 . This 



is w hy th e first two equations in t he hi erarchy, (3.16) 
and ( 3.17 ), coincide with the o.d.e. ( 3.22 ) for the shape 
mode. 

3) Based on numerical data, Powell et al. |69) have con- 
jectured, that <f> converges along the trajectory in func- 
tion space formed by the <&t,'s with v < v* . We here have 
derived this result analytically, and identify the velocity 
v of the transients <fr v with the actual instantaneous ve- 
locity v — v* + X of the front. We find a non- vanishing 
correction of order 1/i 2 to <j> as & v » +x - 

4) The transients $ v have always v < v* at late times, 
since we will find that c\ — —3/2, in accord with the 
discussi on o f Section III A 1 Note that as discussed in 
Section [II B|, such are positive from £x — > — oo up to 



C* = df + 2% + 1 , <F = (a£ x + 0) e~ u . (3.32) 
We remove the exponential through the transformation 

C* = e~^ x 9| e^ x , r]2± = e~^ x , <f> = e~^ x V • 

(3.33) 



The differential equations determining the ip& are explic- 
itly 



a finite value of £x only. For the transient ( p.3l|) wc need 
only the positive part of <lv The transient ( 3.31 ) crosses 
over to a different functional form, before becomes 
negative. 

5 ) Th ere is a non-universal contribution of order l/t 2 
to (|3.31 ). It is non-universal, because it depends on ini- 
t ial c onditions: The structure of our expansion ( [3.14 ), 
( |3.15 ) is an asymptotic expansion about t — > oo, that 
does not fix t — 0. We thus can expand in l/(t— t ) = 
l/t + to/t 2 + 0(l/t 3 ) just as well as in l/t. This all ows u s 



Ci« x +7) ) 7 = /?-", 
C|(a^x +7) , 

[-1 + ci(l - d £ )]tpi + c 2 {aix + 7) , 

[-j +Cl (l_^)],/, s+C 3(l-^)Vl + 

+C5(a£ x +7) 

etc., generally: 
[-V + ci(l- Wa=a 



(3.34) 



n-2 

E« 

rn—3 



.(1 - +C^(aCx +7) 



where we have omitted exponentially small corrections 
of order e~^ x in the inhom ogen eities on th e r.h.s. of the 
equations. The conditions ( |3.9| ) and ( 3.1 0| ) on r/ do not 
influence the solution in the leading edge. 



to add an arbitrary multiple of Vsh /t 2 to <\> in Eq. (plh. The equations (]3J3J) are easily solved. For $ 



The order l/t term in (3.31) is thus always non-zero, 
because the functions p(£x) in (|3.32j ) is non- vanishing 



and not a multiple of T] s h,(£x)j but its precise value will 
depend on initial conditions. 



we find in the region £x 3* 1 



OO / 

2-^1 in/2 
n=2 



(3.35) 



6) The expansion is an asymptotic expansion [106|. 
Thus, when we w ill have determined the coefficients c\ 
and C3/2 in ( |3.14 ) later, these are the exact prefactors if 
we expand the velocity and shape in inverse powers of 
t in the limit t — > 00. However, the expansion will not 
have a finite radius of convergence in 1 / \ft. 



C. Interior shape expanded towards the leading edge 



We now will se e, tha t for £x > 0(y/i) the structure 
of our expansion ( 3.1 5| ) breaks down. We then have to 
resum the terms and use a different expansion 6 . 

Let us calculate the contributions rji from (|3.16|) - 
(3.20) explicitly in the leading edge r egion , where £ x ^* 1 
and </»,$* < 1. In this region C* (|3T3|) and $* fl2.17|) 
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Obviously, for £^ > y/i, the expansion is not properly or- 
dered in powers of 1/y/i anymore, since, e.g., £jf/t even- 
tually will become larger than £x ■ A quick inspection of 



Actually, the interior expansion also breaks down for £x — > 
—00. There too, a different expansion can be used, and this 
expansion can be matched to the one we introduced for the 
interior region. We will not discuss this further here, as it is 
of no further consequence. 
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( |3 .35 ) shows that we can continue to work in an l/y/i 

expansion if we use the v ariable z = ^ ( |3.8| ) instead of 
£x- The expression ( 3.35| ) can be identified with 



, Cl ( Cl -l)( Cl -2)(4z) 7 / 2 , 



7! 

+t „ ^ +a cg -° )(te) + "°' 4 : )3/ ' + ^') N 



2! 



3! 



+ 0(1/Vt 



(3.36) 



This resummed expansion anticipates the crossover to 
the expansion in z and 1/y/t below for £x>^ 3* 1, 
z = £i:/4i = 0(1), which will fix the coefficients C\, 
etc. Note that for c\ = —3/2, terms of order y/i sum up 
to ay/ Azt e~ z = ^xe~^/ 4 ' , which is, in dominant order, 



the behavior already anticipated in Section Iff A 1 



Instead of resumming the interior expansion explicitly, 
it is much more transparent to write an expansion di- 
rectly in terms of powers of 1/y/t and the similarity vari- 
able z of the diffusion equation. This approach, which 
amounts to a matching procedure, is the subject of the 
next subsection. 



D. Analysis of the leading edge 

We now take up the analysis of the leading edge region 
£x > 0(y/t) in the case that the initial conditions are 
sufficiently steep, so that for tp = <f> e^ x 



lim ip(£x,t) < e 



S > 



(3.37) 



Note that according to the discussion of Section II E , this 



condition holds at any finite time t < oo if it is obeyed 
initially at t = 0. 



We have already argued in Sections [IE and III A 1 



that the asymptotic profile of the leading edge might be 
expected to be somewhat like a Gaussian in £x and t 
times a Hermite polynomial. Also the resummation of 
the interior front solution suggests such a form for large 
£x- We now investigate this expansion more systemati- 
cally, and will show that it actually takes the form of a 
Gaussian times a generalization of Hermite poly nomials, 
namely confluent hypergeometric functions 104] 



In passing, we stress that the arguments from III A 1 



can be compared directly to our calculation here only to 
lowest order, because we now work with the coordinate 
z = (x — 2t — A(i)) 2 /(4i), while we presented our earlier 
intuitive arguments in the coordinate z* — (x — 2t) 2 / (At). 
Of course, one can also set up a systematic expansion in 
the latter coordinate z* , but this requires the introduc- 
tion of logarithmic terms for a proper matching to the 
interior part of the front. Working throughout in the 



shifted frames £x = x — 2t — X(t) or z avoids this alto- 
gether. 

In the coordinates £x and t, the equation of motion 
for tp in the leading edge region i s (Re call that the earlier 
leading edge representation in ( [2.49] ) was in the frame 
i = x-2t) 

d t /ip = d\i\) + A(<9 e - 1)V> + o(e~ 5x ) . (3.38) 

The differential operators transform under change of 
coordinates to z = /(4t) and t as 



d t \ g = d t \ 



z » I 
t '* 



d i \ t = W- d 



(3.39) 



Motivated by the form ( p. 45 ) and the discussion of Sec- 
tion 



IIIA1, we extract the Gaussian e £^/( 4t ) 



from ij) by writing: 

W,t) = e" G(z,t) > * = (3-40) 
This extraction also allows us to make contact later with 



funct ions tabulated in [104]. The dynamical equation 
( |3.38 ) is equivalent to the equation for G: 



zdl 
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ci 



G 



(Xt- Cl )+xVt^{l-d z ) 



G 



(3.41) 



The equation is organized such, that the differential op- 
erators of order t° are on the l.h.s. of the equation, while 
the r.h.s. has the operators of order t~ x / 2 and smaller. 

In analogy to our earlier expansion (3.15), we now 
make an ansatz for G in powers of X/yi times functions 
of z. A glimpse at the form of the i nterior shape ex- 
panded towards the leading edge ( 3.36 ) tells us, that the 
expansion should start with the order yft. We write 



G(z, t) = Vig^_ (z) + g (z) + ^£ + 



(3.42) 



Insertion of this ansatz into (3.41) again results in a hi- 
erarchy of ordinary differential equations, that can be 
solved successively: 



zd 2 



zdl 



zdl 



z\d z 



ci 



ci 



9=1 =0 



9o 



cs + ci\/z(l - d z ) 



1 



z 

2 



C2 + cs\/z(l — d z ) 



(3.43) 
(3.44) 

(3.45) 



C3 + Ciyfz(l - d z ) 



9a 
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etc. The general solution of the homogeneous equations 
with two constants of integration ka and i» can be found 

j 1 ° 2 2 

in [104], they are confluent hypergeometric functions. 
These special solutions of the inhomogeneous can also 
generally be expressed in terms of double integrals over 
known functions, as is discussed in Appendix |]. Below 
we will, however, just guess the series expansion of the 
special function g^ p we need. We write the general solu- 
tion as 



- i 2 — n 3 
i f ^Af( Cl + — -,2 



(3.46) 



where the funct ions M(a,6, z) can be expressed by the 
Kummer series [104] 



, ,, , . „ a z a(a + 1) z 2 



(a) n z n 

(b) n n\ 



+ 



with (o)„ = f[(a + k - 1) = r( °,t" } • I :! - 47 ] 



fe=i 



r(a) 



Just as in the integration of the interior shape in Sec- 
tion IIIB. there are two constants of integration to be 
determined in every solution g™ . In addition, however, 
the Cj are not just parameters of the equations as in Sec- 
tion III B, but they now have to be determined also. The 



conditions we use to dete rmine these three constants per 
equation, are now (3.36) and (3.37) in analogy to the 
two conditions (3.9) and (3.10) for the r)&: (i) The solu- 



tion 5™ has to agree with the e xpansion of the interior 
towards the leading edge ( 3.36|) fo r z <C 1. Then the 
coefficients of z° and z 1 !" 1 in ( 3.36] ) determine the con- 
stants of integration fca and ia. (ii) The transients at 
any finite time have to be sufficiently steep in the sense 
that they obey ( 3.37 ) at an y finite time t. Because of the 
form the expansions ( |3.41 ) and (3.43), we require that 
each term g in the expansion diverges for z > 1 at most 
as a power law of z, not exponentially as e z . In addition, 
tp = e z G should not diverge as t — > oo, but approach 
a time independent limit. This gives another condition 
on the constants of integration, that can be obeyed only 
for a particular choice of Cn+3. With these choices of the 

const ants , the small z expansion ofV> = e^Gfrom ( |3~40| ) 
and ( 3.42 ) becomes identical with the interior shape ex- 
panded towards the leading edge (3.36) 



We will solve the first two equations (|3.43| ) and ( |3.44| ) 
explicitly, since they determine the universal terms of 
the velocity correction X. In particular, the solution for 



9=1 ( |3.43| ) will connect to our qualitative discussion of 



the leading \jt veloci ty convergence term (3.5) in Sec 



tion 



III A 1 



Eq. ( |3.44 ) will give the universal subleading 
term 7 of order l/i _3/2 . 

Let us now start with the solution of the homogeneous 
leading order equation (3.43), where g sp 1 (z) = 0. The 



constants of integration are fixed by (3.36) as k^i = 
and l=± = 2a. Therefore q-i (z) is after matching to the 
interior 



g-i(z) = 2aVzM[ci + 



3 3 

2' 2' 



(3.48) 



In order to analyze, how c\ is determined by the match- 
ing and the requirement that all transients are exponen- 
tially steeper for £x — ► oo than the asymptotic profile, 
we first rec all t he large z behavior of Kummer functions 
M(a, b, z) [104]: For positive b each term of the series 
fl3.47p is finite. For a not zero nor a negative integer, the 
series is infinite. For a zero or a negative integer a = —n, 
the series is finite, since all terms from order z n+1 on 
contain the factor (a + n) = 0, and for b = 1/2 or 3/2, 
these finite polynomials arc Hermitc polynomials. The 
large z asymptotics of M(a, b, z) for positive b is 



i-b 



M(a, b, z) 



T(b) 

(a) H * |a| 
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- a G A/a 



(3.49) 



where A fp de notes zero plus the positive integers. If one 
inserts (3.4E) into ( |3.48| ) ones find for £_ x > \/i 



(e x /t) Cl - Cl - 3/2 £ .A/a 

(£lAr ci ~ 3/2 e~«i/( 4t ) -ci - 3/2 g Mo 

(3.50) 



For — ci — 3/2 not a positive integer, we see from ( 3.50| ) 
that 4>(£x,t) does not converg e exp onentially fast to zero, 
in violation of the condition ( 3.37 ) . Accordingly, f or so - 
called sufficiently steep initial conditions that obey (3.37) 
we conclude that c\ -I- 3/2 has to be zero or a negative 
integer. Possible solutions are 



ci = =£ , gl^{z) = 2ajz~ (f-f) , 
etc., with the g=i given by Hermite polynomials. 



(3.51) 



7 DonJt confuse the expansion in 1/ s/i of the vel ocity in ( 3.5 ) 
or (3.14) with the denominators in ( [j.2j ) and (3.3). These 
powers of l/\/t in the ^-representation are absorbe d in to the 
X(t) of ^ in the ^-representation, as is sketched in (3.4). 
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There are two ways to argue, why generically ci = 
—3/2 is the appropriate solution, (a) If the initial condi- 
tion is always non- negative, e.g., because is a density, 
the transient may not have nodes, so c\ — —3/2 is the 
only possible solution, (b) If one can create a front with 
nodes whose leading edge after some evolution is the su- 
perposition of the solutions in (3.51), the solution with 
ci = —3/2 propagates quickest, so the other contribu- 
tions will be convected to the back, a nd th e c\ = —3/2 
solution will dominate at large times 1 107 ] . This argu- 
ment coincides with the argument from Section |II C| , that 
fronts with nodes generically are not attractors for the 
l ong time dynamics for the nonlinear diffusion equation 
(1.1). A similar reasoning for the leading edge r egion 
can be developed from the arguments in Section VI F . 
Furthermore we have checked various initial conditions 
with nodes numerically and we have found, that either 
the node gets stuck behind the evolving front or moves 
away to £ — > oo with velocity larger than v* , leaving in 
both cases a leading edge of the front behind that devel- 
ops with ci = —3/2. We thus find for initial conditions 



(3.37) steeper than $* generically 



ci 



~ 3 f \ 



2a 



(3.52) 



This sol ution is identical with the order yfi of i/ie z with 
ip from ( [3.36] ). For <j) we find in the region £x 3> 1 lin- 
earizable about the unstable state in leading order 

= atx e- c *-&/<**> (l + 0(1/60 + 0{1/Vt)\ (3.53) 



6c = x - v*t + - hxt+ 0(1/ y/i) 



(3.54) 



consistent with the arguments from Section III A 1 



Integration of go now give s the sublc adin g universal 
ter ms, w hich a rc (1/ y/i) in ( 3.53| ) and ( |3.54 ). Insertion 
of J3 52|) into 03.44) results in 



z\d z + l 



(3.55) 

We now can follow Appendix |] for the general solution 
of the inhomogeneous equation, or we rather can guess a 
special solution of the inhomogeneous equation by noting 
that the function 



(l)„-2 Z* 



(3.56) 



is proportio nal to a truncated Kummer series 
M(-l,\,z) ( |3.47|) and solves 



zd 2 



z)d z + l 



F N (z) 



y N-l 



ih) N -r (N-l)' 



(3.57) 



The special solution of the inhomogeneous equation 
03.55 ) is then easily seen to be 



/A,"!:) =2n [^ + 2cs^-z- 3 -F 2 (z) 



(3.58) 



Upon comparing ( |3.56 ) to (3.47) and ( |3.49 ), one finds 



°-aJ^ z- z ' 2 e z 



(3.59) 



The general solution (3.46) of ( 3.55j) is thus 



-1 3 

9o(z) = g s p (z) + k a (l-2z) + l ^M [— ,-,z 



(3.61) 



- I - , - ^ ) : -'• 2 < 



where we have used ( p.49| ) and ( 3.5£ ) for the la rge z 
asy mpto tics. Compare now the small z expansion (3.6C) 
to ( 3.36| ). One obviously has to identify 



3a 



+ k = (3 , 4ac3 + l = 



(3.62) 



If qp w ould decay asymptotically as z 3 / 2 e z for large 
z (3.61), the subleading contribution of order 1/y/i in <fi 
( |3.53 ) would not decay like a Gaussian e 

-e x /m as the 

leading order term does, but it would decay algebraically 
like £y 3 (4t ) 3 / 2 . This would destroy the ordering of our 



expansion (3.42) and lead to a divergence of ip for t — > oo. 
Th us the coefficient of the leading order term z~ 3 ' 2 e z in 
go (3.61) has to vanish: 



a y/ir - 



lo 
4 







(3.63) 



Eqs. ( 3.6 2| ) and (3.63) fix all constants ko, lo and ca. The 
velocity correction of order 1 /t 3 ^ 2 is 



3y/^ 
C3 = — — , 
2 2 

and the analytic solution for go (z) is 

F2(z) 



(3.64) 



go(z) = f3 (1 -2z) + 3a [z- 



(3.65) 



with a and (3 the coefficients of the asymptotic leading 
edge shape $*(£) = (a£ + /3) for $* < 1. Note, that 
the subleading term /3 contributes only the rather trivial 
(1 — 2z) term, while the coefficient of the leading a con- 
tains all nontrivial terms. The result (3.64) and ( 3. 65] ) 
reproduces the order t° in ( 3.36| ) identically. 
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We summarize the results obtained from the analysis 
of the leading edge: The appropriate coordinate system is 
£x = x — v*t — X(t) 1 and the universal velocity correction 
is given by 



X 



3 
2t 



O 



(3.66) 



The shape in the leading edge, where ^ 1, is given in 
terms of the variables £x and t by 



<t>{tx,t)=e-tx-tyW G ^k^ (3.67) 

= e -tx-&/m ^ x + go ^ + (3 . 68) 



i (e 



x 



with ffp( z) fr om 

Eqs. (3.65) - ( |3.67| ) is the second part of our final re- 
sult, valid in the leading edge of the front, where 4> <C 1. 




It complements our earlier result ( |3.31 ), va lid in the in- 
terior of the front, with functions n v from (3.22) and p v 
from (pOTh. 



E. Summarizing remarks 

Let us end this section by putting these analytical re- 
sults into perspective: 

1) The requirement that the leading edge remains 
steeper than the asymptotic profile <j>* at any finite time 
together with the requirement that it converges to 
as t — ► oo determines the velocity convergence constants 
c « . These constants arc thus determined in the leading 

2 ..... . 

edge by the initial conditions. They are just parameters 
in the equations for the interior (3.16) - (3.20). 



2) The leading order velocity correction c\ reproduces 
Bramson's result [J77j, which he derived through solv- 
ing the (nonlinear) diffusion equation with probabilistic 
methods. The universal subdominant 1/t 3 / 2 is new. 

3) According to our discussion in connection with the 
interior expansion, gi and C2 should be termed non- 
universal, because the change from l/\ft to l/\/t — to 
in the asympto tic ex pansion about t — > oo changes these 



terms. As for ( |3.3lD , we conclude that at least parts of 
these terms depend on initial conditions and are therefore 
non-universal. 

4) We stress once more that the full expansion is only 
asymptotic in 1/Vt, but that the prefactors of the 1/t 
and 1/t 3 / 2 terms are exact. 

5) The leading edge expansion is an intermediate 
asymptotics in z valid for 1 <C z <C \fi or \ft <C £x "C t, 
resp. Above, we extensively made use of the cross-over to 
the interior expansion for z«l. Let us now look into the 
break-down for z > 0(\ft), i.e., for £x > 0(t). This sec- 
ond breakdown immediately follows from inserting into 



(13.42 ) our results q-i (z) — 0(\fz) and go{z) > 0{z) (in 
fact go(z) = O(zlnz) according to Appendix ||]) . This 
new crossover actu ally ne eds to exist in view of our dis- 
cussion in Section II E 1 : The steepness A is conserved 
for i -too for all times t < oo. It will retain the infor- 
mation about the precise initial condition. This region of 
conserved steepness at £x > 0(t) crosses over to the uni- 
versal Gaussian leading edge region for £x < 0(t), which 
determines the universal relaxation behavior as discussed 
above. The region of conserved stepness A at £x > 0(t) 
has no further consequence for the dynamics, if the ini- 
tial steepness is only A > A*. It will disappear towards 
tlx —>■ oo by outrunning the leading edge region with an 
approximately constant speed. This scenario is sketched 
in Fig. 2. 

6) Our result is valid in the pulled regime but it does 
not apply at the bifurcation point fr om t he pulled to 
the pushed regime. For nonlinearity ( 1.1 Cl| ) this means, 
that the analysis applies for e > 3/4 J67[. Only then 
$*(£) oc £ e~£, which is one of the essential ingredients 
of our asymptotic analysis. For e < 3/4, the front is 
pushed, and convergence is exponential, as discussed in 
Sections [II C and II D. For e = 3/4, precisely at the 
pushed/pulled transition, $*(£) oc e~^. In this case, con- 
vergence is still algebraic, but the analysis of this chapter 
does not apply exactly. The convergence analysis, how- 
ever, can be set up along the same lines . As shown in 



Appendix ^ we then get instead of ( 3.66Q 



(3.69) 



Note that the factor 3/2 of the 1/t term is replaced by 
1/2 at the bifurcat ion poi nt. Along the lines of the ar- 
guments of Section [II A 1| this can be understood simply 



from the fact that at the bifurcation point the asymp- 
totic behavior of $* is as $*(£) ~ e~ ? , not as £ e ~ 5 > 
and hence that the simple Gaussian leading edge solution 
e~^~£ ^ 4t /Vt matches to the asymptotic front pro file in 
leading order. However, the velocity relaxation ( 3. 69] ) 
at the pushed/pulled transition does contain a univer- 
sal subleading term of order 1/t 3 / 2 that is absent in the 
relaxation of the linear equation ( [2.46 ) . 

7) Up to now we excluded the particular initial con- 
ditions cj)(x, 0) ~ x~ v e~ x from our discussion, since they 
are neither sufficiently steep nor flat according to our 
definition. It is amusing to see that also such initial 
conditions can be treated with our approach. For suf- 
ficiently steep initial conditions, we discarded the case 
that — ci — 3/2 would be different from an integer or zero 
after Ecl^ 3.5C ) , because it would violate the exponential 
bound (3.37) for large £x- However, for the above par- 
ticular initial conditions, the asymptotic behavior (3.37) 
is replaced by ijj as £ x " for £x > 1. For any v < 2 one 
concludes immediately from Eq. ( 3.50| ) that 



#c,0) 



v(t) = 2- 



2t 



(3.70) 
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a result also derived by Bramson |J77| . In other words, in 
the case in which the initial conditions are intermediate 
between sufficiently steep and flat, the pref actor c\ does 
depend on the initial conditions and may even change 



sign, but the relaxation is still power law like. To get 
the next order term in the expansion for these special 
initial conditions, our expansion will probably have to be 
generalized. We will comment on this in Section (v|. 
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IV. SIMULATIONS OF PULLED FRONTS IN 
THE NONLINEAR DIFFUSION EQUATION 

In this section, we present simulation data for fronts in 
the nonlinear diffusion equation d t (j) = d^cf) + f{4>) (1.1) 
propagating into the unstable state = 0, and com- 
pare these with our analytical predictions. In particular, 
we thoroughly investigate fronts with the nonlinearity 
f(4>) = 4> — cj) 3 , so that the equation becomes 



(4.1) 



This equation forms pulled fronts with v* — 2 and 
A* = 1 = D, if the initial conditions are sufficiently 
steep. As an example of a nonlinearity allowing for 
both pushed and pulled fronts, we also present data for 
/(</>) = ecj) + 3 - 5 for e = 0.56 and 0.96. 
As an initial condition, we here always choose 



#r,0) 



1 



1 -|- gXinit(x—X. ) 



e ->-init{x-X ) x 

1 X 



(4.2) 



According to our analytical results, all initial conditions 
with initial steepness A* < Xinit _ oo exhibit the same 
universal relaxation behavior asymptotically as t — » oo, if 
the front is pulled. We indeed do find this in our simula- 
tions. Below we only present simulations for Xi n u = 10. 

The section is organized into a discussion of t he spe - 
cific numerical features of pulled fronts (Section IV A ), 
the presentation of the raw simulation data for nonlin- 
earities f(<f>) — <fi - <P 3 and f(<f>) — ecj) + <fr 3 - 5 (Section 
IVB ), and a detailed comparison of the sim ulatio ns for 
(1.1) with the analytical predictions (Section IV C|) . 



A. Numerical features specific to pulled fronts 

To integrate a given initial condition <p(x, 0) forward in 
time t for a nonlinear diffusion equation, we use a semi- 
implic it alg orithm which is explained in detail in Section 
VF6j Eq. ( pM14| ). When running the program, we have 
to choose a spatial and temporal discretization Ax and 
At, a system size < x < L, and a position xq of the 
initial condition within the system. Comparing results 
for different parameters Ax, At, L, and x to each other 
and to the analytical predictions in the extreme preci- 
sion of often better than 6 significant figures, we find two 
features specific to the particular dynamic mechanism of 
pulled fronts: 



1. Effect of finite difference code 

The numerical results of the simulation depend of 
course on the step sizes Ax and A t of the finite differ- 
ence code. In fact, in Section VF6 we will have collected 



all analytical tools to calculate the corrections to v* = 2, 
A* = 1 and D = 1, that depend on the numerical integra- 
tion scheme and on the parameters Aa; and At. All data 
presented here are derived for Air = 0.01 = At. For a 
pulled front in a nonlinear diffusion equation solve d with 
a semi- implicit scheme, our analytical prediction ( 5.116| ) 
yields v* = 2.000075, A* = 0.999954, and D = 1.00035. 



2. Effect of finite system size 

In contrast to a pushed front, the final t — > oo relax- 
ation of a pulled front very sensitively depends on system 
size L and front position xq . This effect is closely related 
to the pulled mode of propagation and the breakdown 
of the linear stability analysis. Because the half- infinite 
space x ^> 1 of the leading edge dominates the dynam- 
ics, the very long time dynamics of the front is sensitive 
to the region at i > 1, even if though there ^ < 1. 
More precisely, the d iffus ive s pread ing of the linear per- 
turbation as in Eq. ( |3.3| ) or (3.67) that determines the 
speed, strongly depends on the boundary conditions at 
£x = x - v*t - X(t) = 0(VtDi ). 

For this reason, we shift the front back to its origi- 
nal position £o within the system after every time step 
ti — t\ = 1. This eliminates the a;- interval < x < 
x shift ~ v* on the back side of the front from our 
data, while a new x- interval L — x shift _ % < L 
has to be created. One might assume, that this pro- 
cedure yields good results for integration times T up 
T = O ((L- x ) 2 /(4-D)) because of the diffusive na- 
ture of the spreading. However, the precision notice- 
ably breaks down earlier because of the arbitrariness of 
the newly created x- interval L — x shift < x < L m 
the shift process. Filling this region with the constant 
4>(x) = 4>{L — x shift) creates a flat initial condition, and 
the front accelerates beyond v* for sufficiently long times. 
We therefore use <p{x) = in this region. The observed 
velocities v<j,(t) for L < oo then always will stay below 
those in the infinite system L — > oo. The simulations in 
the finite system are close to those in the infinite system 

up to times T = o((L - x )/v* 



B. Simulation data 

1. f((f>) = — 4> 3 : pulled fronts 

As a n ex ample, we will extensively discuss simulat ion s 
of Eq. (4.1). We present data with initial conditions (4.2) 
and Xinit — 10, where the initial condition is located at 
x n = 100 in a system of size L — 1000. According to our 
estimate above, the simulations then should be reliable 
up to times t of order (L — xq)/2 — 450. We present data 
up to t = 400. The data from this simulation is eval- 
uated in a sequence of figures showing increasing detail 
and precision. 
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Fig. 1 already showed the temporal evolution of a suf- 
ficien tly s teep initial condition under the equation of mo- 
tion (4.1). It shows both, the total displacement of the 
front, and the evolution of the front shape. We now 
choose different presentations that show these two dif- 
ferent aspects of the dynamics separately and in higher 
precision. 

Let us first study the evolution of the front shape: In 
Fig. 4, we present <f>(£x,t) as a function of £x> where 
£x = x — v*t - X(t), Eq. fT?), is adjusted such that 
<f>(0, t) = 1/2 ( fh9| ) for all times t. The remaining dynam- 
ics in this frame is then the pure evolution of the shape 
from its steep initial profile </>(£x , 0) towards its flatter 
asymptotic profile — > $*(£x) as t — > oo. Fig. 4(a) 
shows as a function of £x on the interval — 5 < £x < 5. 
One sees the interior or nonlinear part of the front. Fig. 
4(6) shows \og(p in the range 10 -90 < 4> < 1. This plot 
is appropriate to show the development of the leading 
edge, which here essentially determines the dynamics. 
Accordingly, a very different range of £x has to be plot- 
ted, namely < £x < 190. As is sketched already in Fig. 
2, the leading edge consists of two regions, namely the 
"Gaussian" region, through which the asymptotic steep- 
ness A* spreads in time towards larger , and the region 
of conserved steeness A = Xmit in front of it. In fact, 



Fig. 4(6) shows that the initial Xinit = 10 on the level 
<f> = 1CP 90 is still fully present for times t = 1 and 2, 
while at later times it gradually approaches A* = 1. At 
higher levels, <f> = 10~ 10 say, this process of replacement 
of one steepness by the other is essentially completed at 
time t = 70, while at level 10 -90 , it is not completed even 
at time t = 400, where the simulation stops. 

In Fig. 5, we focus on the second feature, namely the 
displacement of the front. We plot the velocity vj,(t) of 
various amplitudes <f> as a function of t. According to our 
previous definition, we identify Vi/ 2 {t) = v* + X(t). For 
comparison, the predicted asymptotic value v* is plotted 
as a dashed line. In Fig. 5(a), the non-universal initial 
transients up to time t — 20 are shown on the range 
< v < 3. In Fig. 5(6), the velocities are plotted up to 
time t = 400 on the velocity interval 1.97 < v < 2. One 
observes, 

- that for fixed t, the velocity v^,(t) is the smaller, the 
larger (j) is. This is an immediate consequence of the 
fronts becoming flatter in time, cf. Fig. 4. 

— that the v^(t) for large t approach a value largely in- 
dependent of 4>, that is still far from the asymptotic value 
v* . We will see below, that this is the signature of the 
shape relaxing like v ( j >1 (t) — v^, 2 (t) oc 1 /t 2 as t — > oo, while 
the overall relaxation is v^(t) — v* oc 1/t. 



,3,5,7,10,20,40,70,100, 

140,200,250,300,400 




FIG. 4. Simulation of the evolution of the shape of a front under (^l]) at the times denoted in the figure. The initial condition 
is (4.2) with A = A in4t = 10. The comoving frame £x is chosen in such a way that 4>(£,x = 0,t) = 1/2 for all t. a) A plot of (j> 
versus shows mainly the interior of the front, b) A plot of log <f> versus £x for sufficiently large shows mainly the leading 
edge of the front. Note the different scales of £x- 
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FIG. 5. The same simulation as in Fig. 4. Now the velocities v$(t) of amplitudes 4> = 0.99, 0.9, 0.7, 0.5, 0.3, 0.1, 0.05, 0.01, 
0.001, and 0.0001 (solid lines) are shown as a function of time t. The asymptotic velocity v* is marked by the dashed line, a) 
Initial transients for times < t < 20. b) The same data plotted for longer times < t < 400 on an enlarged scale of v. The 
velocities v</>(t) become largely independent of the "height" <f), and together slowly approach v* in agreement with the predicted 
universal algebraic relaxation. 



2. /((/>) =e<j) + (p a - (f> 5 : pushed versus pulled fronts 

A well-known example of a nonlinear diffusion equation 
(1.1) exhibiting both pushed and pulled fronts is given by 
the nonlinearity ( l.lOj ) with n — 2: 



d T ip = dyip + e(f + ip 3 



(4.3) 



This equation for e < is often used as a phenomeno- 
logical (Ginzburg-Landau type) mean field model for a 
first order transition. Likewise, its extension to a com- 
plex field is often used to model a subcritical bifurca- 
tion in pattern formingsystems. Acc ord ing to arguments 
recalled in Appendix fronts of (13) are pushed for 
e < 3/4, and pulled for e > 3/4. 

Th e res caling necessary to b ring (|4.3|) to our standard 
form ( |2.2| ) is discussed in ( |2.3| ) and Q2.4| ), and yields 



di 



when 



i 




_ VTTTe-l 2 

6 = o > ¥>8 = 1 



(4.4) 
(4.5) 



The critical e, where the pushed/pulled transition occurs, 
is e c = 0.5. 

We present data for the pushed front with e = 0.4 
(e = 0.56) and the pulled front with e = 0.6 (e = 0.96). 
The initial condition is the one given before in (|4.2| ) . The 



system size is L = 250 and the front is located at xo = 50. 
The data therefore should be reliable up to time of order 
100, so the data presented extend over < t < 100. 

In Fig. 6, we plot v<j,(t) as a function of t for both val- 
ues of e as solid lines, in the same way as the plot of Fig. 
6 for the other nonlinearity. The dashed lines denote the 
asymptotic pulled velocity v* = 2 predicted for e = 0.6, 
and the asymptotic pushed velocity (cf. Appendix ^|) 



l + 4e 



2.00594 for e = 0.4 



(4.6) 



Fig. 6 shows (i) that the simulated fronts in fact do ap- 
proach the predicted asymptotic velocities, (ii) that up 
to time t < 10 both fronts show quite similar initial tran- 
sients, (Hi) that for time t 3> 10, however, the relaxation 
towards the asymptotic velocity for e = 0.4 is much 
more rapid than that towards v* for e = 0.6. This very 
clearly illustrates the difference between pushed expo- 
nential and pulled algebraic relaxation, despite the tiny 
difference between v* and . 

We do not plot the figures of shape relaxation equiv- 
alent to Fig. 4, the only difference being that in the 
pushed case the region of conserved steepness A = Xmit 
at £x ^> 1 is invaded by a region of steepness rather 
than by the pulled steepness A* . A^ is determined by the 
front interior, while A* is determined by the Gaussian 
region of the leading edge. 



33 



20 



40 



60 



80 



100 



FIG. 6. Plot of v<f,(t) as a function of t as in Fig. 5(6), but now for Eq. (4.3). Simulations for e = 0.4 (e = 0.56) and 
e = 0.6 (e = 0.96) are shown. The dashed lines denote the asymptotic pulled velocity v* = 2 of the front with e = 0.6, and 
the asymptotic pushed velocity v 1 = 2.00594 of the front with e = 0.4. Note the quick exponential relaxation towards tft in 
contrast to the slow algebraic relaxation towards v*. Further away from the transition e = 0.5 (e = 0.75) from pulled to pushed 
front propagation, the relaxation in the pushed regime is even faster and the difference tr — v* is larger. 



C. Comparison of simulations and analytical 
predictions 

We now return to our extensive simulation of the pulled 
front formed by the F-KPP equation dt<p = d 2 4> + 4> — 4> 3 , 
and compare the simulation data to our analytical pre- 
dictions from Table II with v* = 2, A* = 1 = D. 



1. Analysis of the velocity data 

We first concentrate on the analysis of the velocity data 
v<j,(t) from Fig. 5. The prediction for the velocities v$(t) 
of the amplitudes <f> is derived from the expressions in 
Table II through i^(t) = —dt^/dx^QSxed- The result is 



X-X 



Vsh 



O 



> fixed 



(4.7) 



Remember, that X is universal only till order 1/t 3 ! 2 and 
will exhibit contributions in order 1/t 2 , that depend on 
initial conditions. The difference v^. L (t)—v l j >2 (t), however, 
will turn out to be independent of initial conditions up 
to order 1/t 5 / 2 . Let us now test these predictions on the 
simulations in a series of plots with growing precision in 
Figs. 7-9. 

As the velocity correction X is 1/t in leading order, 
we plot Vtf,(t) as a function of 1/t in Fig. 7, for the time 
range 5 < t < 400 in Fig. 7(a), and for 100 < t < 400 in 
Fig. 7(6). The dashed lines present the predicted asymp- 
totes v* + X = 2 + Xi(t) (the lower dashed line), and 
v* + X = 2 + X 3 / 2 (t) (the upper dashed line), where we 
define 



3 

Yt 



x 3/2 (t) 



3 
2t 



(4.8) 



First of all, in comparing Figs. 7(a) and 7(6), we rec- 
ognize the asymptotic nature of the l/\ft expansion: 
whether the X\ or the X 3 / 2 asymptote gives the bet- 
ter prediction, depends on the time scale: If we neglect 
the upper three solid lines with velocities v<f,(t) for the 
very small amplitudes <f> = 0.01, 0.001, and 0.0001, the 
asymptote 2 + X\ clearly fits much better in Fig. 7(a) for 
times 5 < t < 400 — while the asymptote 2 + X 3 / 2 es- 
sentially coincides with «o.ooi(i)> an observation we have 
no analytical explanation for. For times 100 < t < 400 
in Fig. 7(6), however, the coincidence with 2 + X 3 / 2 is 
excellent for all <f), and 2 + X\ very clearly is "far off" on 
this very detailed scale. Hence we will work below with 
the asymptote 2 + X 3 / 2 (t), and we present data for the 
time regime 20 < t < 400 in Figs. 8 and 9. 

Let us no w fu rther zoom in on the 4> dependent velocity 
corrections (L7) to X. Fig. 8(a) shows v ( p(t) — 2 — X 3 



L3/2 as 

a function of X 3/2 = 3 /(2 t 2 ) (l - (3/2)^Jt). Accord 



ing to the prediction (4.7), the plot for small values of 



X 



3/2 



should show essentially straight ^-dependent 



X 



3/2 



as X 



3/2 



0. 



lines, all approaching (t) — 2 
Clearly, that is what they do. 

Fig. 8(6) shows one further step of precision aiming 
now at the precise value of v*: ( [4.7|) predicts 



X 



X 



Vsh 



O 



I <p fixed 



(4.9) 



However, the evaluation of this expression with X 3 / 2 
(4.8) yields (^-independent corrections of order 1/yt: 



X 



3/2 



X 



(4.10) 
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Remember, that the constants C2, C5/ 9 etc . depend on 
the initial conditions. According to ( 4.1Cf) , if we plot 



(v4>{t) — v* — X3/2) I X3/2 as a function of we ex- 

pect these functions to approach a ^-dependent constant 
as -> 0. 

Fig. 8(6) shows, that they in fact do so — but only if we 
choose the correct value of v* ! The dotted lines show the 
function for v* — 2, the fat solid lines for v* — 2.000075. 
The latter value is the analytical prediction of v* taking 
the finite gridsize corrections of t he num eric al cod e into 
account, as explained in Sections IV A 1 and VF6. The 
two values of v* differ in the 6th significant figure. Fig. 
8(6) thus is an extremely precise demonstration of the 
correctness of our analytical arguments from both Sec- 
tions [FFi] and [v|, since it clearly confirms our predictions 
to more than 6 significant figures! 

Our test in Fig. 8(6) is so sensitive, because we divide 
in Fig. 8(6) by the small quantities X 3 / 2 , which are of or- 
der 10 -5 . Without this division the difference of the u*'s 
in Fig. 8(a) is not yet visible. The plot in Fig. 8(6) shows 



that we fully understand the specific numerical features 
of pulled front solutions, both the effect of the finite dif- 



ference code and of the finite system size, cf. Sect. IV A. 

We can eliminate v* and the nonuniversal corrections 
™c 2 /i 2 etc. by plotting (v^t)- v Q . 5 (t))/X 3/2 (t) as a func- 
tion of 1/t. Now (4/7) predicts 



v<t>(t) ~ v . 5 (t) 



Vsh 



X 



3/2 



(4.11) 



Fig. 9 shows this plot with the solid lines for cf> = 0.99, 
0.5, 0.01, and 0.0001. The crosses on the axis are not(!) 
extrapolated from the curves, but they mark the pre- 
dicted asymptotes - r] s h/d^*\ ij) for = 0.99, 0.5, 0.01, 
and 0.0001. The necessary data on r) s h(£) and $*(£) are 
derived from the numerical solution of the appropriate 
o.d.e.'s, and completely independent from the numerical 
integration of the p.d.e. for the initial value problem. The 
coincidence of the extrapolated p.d.e. data with the an- 
alytically predicted o.d.e. asymptote is most convincing. 



a) 5<t<400 b) 1 00<t<400 




FIG. 7. The data v$(t) from Fig. 5, but now plotted over 1/t. The lower straight dashed line is the asymptote v(t) — v*+X\(t) 
(|l.S|), the upper curved dashed line is the asymptote v(t) — v* + X 3 f 2 (t) with v* = 2. a) time regime 5 < t < 400, b) time 
regime 100 < t < 400. Note that due to the i~ 3//2 correction term, the effective slope in this plot is less than 3/2, even at these 
long times. 
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FIG. 8. The data v$(t) from Figs. 5 and 7 for times 20 < t < 400 in different representations, a) Vj,(t) — 2 — X3/2 as a 

X3/2) /X3/2 as a function of l/\/t for = 0.99, 



function of X 3 / 2 - See Eq. (4.8) for the definition of X 3 / 2 . 6) (^(t) 



0.5, 0.01, and 0.0001. Dotted lines: 



2, solid lines: with the corrected value v* 



2.000075 for our numerical scheme and 



gridsize according to Section V F 6 




FIG. 9. The solid lines are again the data from Figs. 5, 7, 8 for times 20 < t < 400, now plotted as (y${i) — wo.s(i)) / ^3/2 
over 1/t for <f> — 0.99, 0.5, 0.01, and 0.0001. As explained in the text, this eliminates a nonuniversal 1/yi term, that depends 
on the initial conditions. The crosses result from solving the o.d.e.'s for $* and r] s h numerically and plotting —T] s h/d^* | for 
4> = 0.99, 0.5, 0.01, and 0.0001. Eq. (4.11) predicts that the lines should extrapolate to the crosses. Since they do, and since 
X(t) is of order 10 -5 at the latest times, these data confirm our predictions with extreme precision. 
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2. Analysis of the shape data 

We now leave the analysis of the velocity data, and 
come back to the shape data from Fig. 4. Table II im- 
mediately yields 



1 



(4.12) 



This gives the clue on how to rewrite the shape data 
<j)(£x,t) at different times as a function of £x- The so- 
lutions of the o.d.e .'s for rj s h and that are needed 
for evaluating (4.12) are derived numerically. They have 
been used for generating the crosses in Fig. 9 and are 
now also used in Fig. 10. 



Plotting the l.h.s. of Eq. ( 4.12 ) allows us to combine 
the information about the interior from Fig. 4(a) and 
the information about the leading edge from Fig. 4(6) 
into one plot. In Fig. 10(a), we do not divide by X, but 



present the data at the small times t = 1, 2, 3, 5, 7, 
10, and 20 as - $*)/ 'rj sh over £ x . In Fig. 10(6), the 
data at the large times t = 20, 40, 70, 100, 140, 200, 300, 
and 400 are shown as ((f) — $>*)/(Xri s h) over £x, where 
we use again the approximation X 



X 3/2 (4 



For 

comparison, both plots also show $*(£x) and £x = as 
dashed lines. Also Fig. 10(6) has the large time predic- 
tion ((f) — 4>*)/ Xrj s h — > 1 as t — > oo as a dotted line. 

Fig. 10(a) shows how the interior of the front rapidly 
relaxes. Fig. 10(6) demonstrates (i) that with X = X 3 / 2 
we indeed have chosen the correct asymptote, and (ii) 
how the predicted asymptotic value (</>—$*)/ (X n] s h) — * 1 
as t — > oo is approached from above in the interior of the 
front and from below in the leading edge. 

Note, that in Fig. 10(6) all lines approximately cross 
one point of height unity far in the leading edge. We 
have no intuitive or analytical understanding of this ob- 
servation. 



a) -(4-4>Vr| sh 
1.5 - 
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FIG. 10. In this figure the shape data from Fig. 4 are represented differently using data from the numerical solution of 
the o.d.e.'s for 3>* and r) a h- a) —(0 — < l > *)/?7sh (solid) as a function of £x for times t = 1, 2, 3, 5, 7, 10, and 20. 6) 
((p - $*)/(X 3/2 r)sh) (solid) as a function of £x for times t = 20, 40, 70, 100, 140, 200, 250, 300, 400. Dotted line: predicted 
asymptote (4>— <&*)/(Xr) s h) — *• 1 as t — *• oo. Dashed lines in (a) and (b) give $*(£x) and £x — for orientation. 
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V. GENERALIZATION OF PULLING TO 
HIGHER ORDER (SETS OF) EQUATIONS 

A. Introduction 

In the last fifteen years, it has become clear that 
many of the observations and intuitive notions concern- 
ing the behavi or of front solutions of the nonlinear diffu- 
sion equation (1.1) generalize to higher order equations 
or systems of coupled p.d.e.'s. First of all, taking the 
spreading velocity v* of a linear perturbation of the un- 
stable state [Eqs. (1.4) and (1.5)] as the generalization 
of v* = 2y/ f'(0) for (1.1 V we observe that there are 
numerous examples P Jl 2| , |l 7| , p2| , |63| , |S6| |68| of fronts whose 
asy mpto tic velocity approaches the pulled value v* given 
by (1.5). So there is no doubt that the mechanism of 



in Table II. The fact that these predictions for various 
examples are fully corroborated numerically in Section 
V F makes us conclude that the velocity selection and re- 



fronts "being pulled along" by the leading edge gener- 
alizes to a large class of equations. Second, there is 
also quite a bit of evidence for the existence of a pushed 
regime in more complicated equations. In a number of 
cases, the pushed regime was again found to be related 
to the existence of a strongly heteroclinic solution with 
velocity i?t > v* . An example of a non-monotonic but 
still uniformly translating pushed front solution in the 
EFK equation is shown in Fig. 7 of j67|. In the quintic 
complex Ginzburg Landau equation, it has turned out 
to be possible to solve for a strongly heteroclinic front 
profile exactly, and in numerical simulations it was em- 
pirically found that this solution does play the same role 
in the front selection process as the pushed front & in 
the nonlinear diffusion equation ^38|. Pushed fronts also 
emerge in coupled amplitude equations for chaotic do- 
main boundary motion pl| . For extensions of the Swift- 
Hohenberg equation there are numerical and perturba- 
tive indications that both pulled and pushed regimes oc- 
cur, and that one can tune the front velocity from one 
regime to the other with one of the nonlinear terms in 
the equation J67| . 

Much of our understanding of the above general find- 
ings has been intuitive and empirical, or based on con- 
jectures. We shall now show that many of our results for 
the second order nonlinear diffusion equation generalize 
to other equations, not only to (sets of) partial differ- 
ential equations of higher order, but also to other types 
of equations like difference-differential equations pi, 109 1 , 
or differential equations with memory kernels [p.!0| . We 
will concentrate here on equations whose relevant front 
solutions are uniformly translating. For p.d.e.'s in this 
class, essentially the whole classification of nonlinearities 
and initial conditions (/>(a;,0) in Table IV applies, pro- 
vided the uniformly translating fronts $„, and in par- 
ticular the fronts <I>* and <^ exist. Many aspects of the 
stability analysis can be generalized, while the relaxation 
of pulled fronts req uires the generalization of the calcula- 
tion in Section [II. This generalization, that we will de- 
velop below, leads to new and explicit predictions for the 
front convergence in the pulled regime, as summarized 



laxation of uniformly translating fronts is now essentially 
understood even for general sets of equations. 

While this paper was nearing completion, it was be- 
coming increasingly clear that even though pattern form- 
ing fronts — both fronts leading to regular periodic pat- 
terns, as in the Swift-Hohenberg equation ]7^ , |6^ , p7t , and 
fronts leading to chaotic patterns as in some parame- 
ter ranges of the complex Ginzburg-Landau equation — 
present additional complications, our most central result 
for the universal algebraic velocity relaxation carries over 
even to these. We will leave this discussion to the fu- 
ture [ [75|j7^1 , and focus here on p.d.e.'s whose asymptotic 
pulled fronts are uniformly translating front solutions of 
the type <j>*(x — v*t), just as in the nonlinear diffusion 
equation. 

In writing this section, we face the following dilemmas: 



(i) The extension of both the stability considerations 
of uniformly translating front solutions of section 
H and of the relaxation analysis of pulled fronts 
of section [II depends quite crucially on two in- 



gredients: First, that the front propagation into 
unstable states is in the pulled regime, and, sec- 
ond, that there is a family of uniformly trans- 
lating front solutions around $*(ir — v*t): Only 
then can the relaxation in the front interior be 
along the manifold of front solutions according to 
<j>{x,t) = $„(t)(6t) + 0(l/t 2 ). However, to our 
knowledge there is no general theory concerning the 
conditions under which fronts are pulled and con- 
cerning the multiplicity of front solutions: For par- 
ticular equations under study or for some restricted 
classes of equations, one can often convince oneself 
that the front should be pulled and that <1>* should 
be a member of a family of front solutions, but a 
general theory is lacking. 

(ii) An immediate jump to the most general (but ab- 
stract) case is pedagogically not justified and more- 
over would assume knowledge of the derivation of 
the pulled velocity v* that most readers probably 
do not have. 

We have chosen to deal with this dilemma by simply 
summarizing our main assumptions and our results con- 
cerning the extensions of Section [n] to more general equa- 
tions below, relegating the details of the analysis to Ap- 
pendices. Then, we proceed with the relaxation analysis 
of p ulled fronts in two steps. We first consider in Section 
V C the analysis of a single p.d.e. which is of first order 
in time but of arbitrary order in space. After that, the 
extension to p.d.e.'s that are of higher order in time is 
discussed in Section VD. The extension to even more 



general classes of equations, including difference equa- 
tions or integro-differential equations, e.g., with me mory 
kernels, is then immediate, as we discuss in Section VE. 
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We there also discuss coupled equations. Section VF 



contains the explicit analytical and numerical results for 
several of the equations listed in Table I. 



B. Basic assumptions underlying the relaxation 
analysis of pulled fronts; generalization of Table IV 

Most of the results discussed in Section || for the non- 
linear diffusion equation can be generalized to higher or- 
der nonlinear partial differential equations, as well as to 
difference or integro-differential equations and to coupled 
equations: 

• The family of solutions can be parametrized as well by 
the steepness A which gives the rate of exponential decay 
of as £ — > oo. 

• If there are one or more strongly heteroclinic solutions, 
then at each velocity where such a solution exists, there 
is a strongly heteroclinic mode of the linear stability op- 
erator which changes stability, i.e., which is such that the 
mode is stabilizing for fronts with velocities larger than 
this value and destabilizing for velocities less than this 
value. This implies in particular that the pushed veloc- 
ity v* is the largest velocity at which there is a strongly 
heteroclinic front solution and that front solutions 
with v < are unstable (see Appendices |H] and |) . 

• Th e line ar spreading velocity v* , given by Eqs. f[l|) 
and ( 5.17 ) below, is the pulled front speed and coincides 
with the minimum of the vel ocities of uniformly trans- 
lating fronts v(X) (see Section V C 2| ). 

• If there are no strongly heteroclinic solutions with 
v > v* , all front solutions with v > v* are stable to 
perturbations which are steeper than A*, while front so- 
lutions with v < v* are unstable: the pulled front solution 
is then the slowest and steepest solution which is stable. 

• The fronts that dynamically emerge from steep ini- 
tial conditions (falling off faster than e~ A x ) converge to 
pulled fronts propagating with speed v* . 

In this Section, we will investigate the front relaxation 
under the assumptions: 



A The front solutions are pulled, i.e., starting from 
a steep initial condition the asymptotic front speed 
Vsei e quals the li near spreading speed v* given by 
Eqs. (5.1 d() and (5.11) below. 



B The asymptotic front is uniformly translating , i.e., 
of the form $*(x — V*t), and it is a member of 
a continuous family of uniformly translating solu- 
tions $> v (x — vt), parametrized by v. 

To put our general assumptions A and B into perspec- 
tive, we note that for a given equation the existence of 
a family of front solutions can often be demonstrated by 
counting arguments. This is shown in Appendix |l| for 
p.d.e.'s of first order in time that are invariant under 
space reflection. Such counting arguments also lead one 
to expect that generically either <I>*(a; — v*t) is a mem- 
ber of a continuous family of front solutions, or there is 



no uniformly translating front solution $* at all. For, 
if there is a discrete set of front solutions (solutions 
exist at isolated values of the velocity), there is no par- 
ticular symmetry reason to have one at v — v* , since the 
existence of an isolated solution depends on the full non- 
linear behavior of the ordinary differential equation, not 
just on the properties near one of the asymptotic fixed 
points. We comment in Section VI on what might hap- 
pen when there is no uniformly translating front solution, 
even though the front dynamics is pulled. 



C. Pulled front relaxation in single p.d.e.'s of first 
order in time 



In the present Section C, we discuss an arbitrary p.d.e. 

d x t)>, d^, d t <t>) =0. (5.1) 

for a single field <fi(x, t). We assume that F is analytic in 
all its arguments, and that the equation admits homoge- 
neous steady state solutions <j> = and 0=1. Moreover, 
we assume <$> = to be linearly unstable and <fi — 1 to be 
linearly stable, and we consid er fr onts connecting these 
two asymptotic states as in (2.14). Also, according to 
our assumption B, Eq. (5.1) admits a continuous family 
of uniformly translating fronts cf>(x,t) — <f> v (x — vt). The 
linearization of some front <j>(x,t) about some <& v gener- 
alizes from Q2.25|) - fl2.27|) to 



(*,t) = *„(£)+ t?(£,t) . d tV = C v (0v + O(r] 2 ) 



(5.2) 



where the linear operator is now 

N 



n=0 



F 



(0 



N+l 



(5.3) 



Here F n (£) denote the functional derivatives of F: 
SF(^°\...,^ N +^) 



(t> {m) =d™<l>v(t), m<N+l 
<£ (JV + 1) =-i;d 5 $„({) 



(5.4) 



In order that the F n have no singularities, i^jv+i should 
be of one sign; for convenience, we take Fjv+i(£) < for 
all £ and rescale time t such that F/v+i(oo) = — 1. We 
also assume, that -F/v(£) neither vanishes nor changes sign 
for any £. 
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1. The pulled velocity v* 

In the pulled regime and with steep initial conditions, 
the asymptotic front velocity equals the linear spread- 
ing velocity v* , i.e., the velocity with which a localized 
perturbation spreads according to the linearized equa- 
tions. Since the calculation of v* forms the basis of our 
subsequent analysis, we summa rize its derivation in the 
context of our fir st ord er p.d.e. (5.1). The general formu- 
lation in Section VE, which is necessary to treat differ- 



ence equations or integro-differential equations, is closest 
to the original "pinch point" analysis [ j58|]60| , from which 
many of these ideas originally emerged. 

In the rest frame (x,t), the equation linearized about 
cb = is 



d t 4> = £o(oo)0 



N 

71=0 



(5.5) 



which is the generalization of ( 2.40| ), and where we intro- 
duced the short hand notation a n = /„(oo). The disper- 
sion relation u)(k) of a Fourier mode e 
by 



ikx—iu(k)t 



is given 



>{k) — 2J a n (ik) r 



(5.6) 



Since we later will again characterize fronts by their ex- 
ponential spatial decay rate A = — ik, we already define 
the growth rate s(A) of the steepness A as 



N 



s(A) = Re = Re > a n (-A) 



(5.7) 



n=0 



for later use. We restrict the analysis to equ atio ns where 
the temporal growth rate Re (-ioj(k)) in (5.6) will be 
negative for short wave length Fourier modes k, i.e., 
where 



Re a N (±i) N < , 



(5.8) 



since otherwise all smooth solutions will be unstable 
against perturbations of arbitrarily short wave lengths. 

A n ar bitrary initial condition <p(y, 0) will develop un- 
der (5.5) as 



<p(x,t) 
G(x,t) 



dy G(x-y,t) </>(y,0) , 



dk 
2^ 



'LKX — IUJ 



(k)t 



(5.9) 
(5.10) 



in generalization of (2.44). 

For sufficiently steep initial conditions (j)(y,0), the 
asymptotic behavior of 4>( x i t) can be obtained fro m the 
large-time asymptotics of the Green's function G ( 5.1C|) 
that can be evaluated by a saddle point integration ]106| 
(also known as "steepest decent approximation"). The 



result will depend on the frame of reference. In an arbi- 
trary coordinate system £ = x — vt with v fixed, a saddle 
point k n is a saddle of —icu(k) + ivk, 



dk 



(—ioj(k) + ivk) 







cLu(k) 



fen 



dk 



(5.11) 



A polynomial of degree N (5.6) generically has N— 1 sad- 
dle points k n , n = 1, . . . , N — 1, ( 5.11 ) in the complex k 
plane. The integral ( 5.10| ) is therefore dominated by the 
saddle point with the largest growth rate through which 
we can lead the fc-contour by continuously deforming it 
off the real axis. If the contour can be deformed to go 
through several, the relevant one is thus that particular 
saddle point k*(v) of the ones we can reach that has the 
maximal growth rate: 



Re (— iu>(k*) + ivk*) — max Re (—iuj(k n ) + ivk n ) 



It will have 



Re D > 



(5.12) 



(5.13) 



k'(v) 



We stress that Eq. (5.11) only expresses the condition 
for the existence of a saddle point. Which saddle point 
the fc-contour can be made to go through by contour 
deformation is a condition that depends on the global 
properties of oj(k) that can only by analyzed for a given 
dispersion relation. It is not a local condition. For a fur- 
ther discussion of this point we refer to Appendix [m| we 
proceed by assuming that the saddl e point indicated by 
a star * is the one that obeys ( 5.12 ) and this condition, 
without distinguishing this underlying condition with our 
notation. 



The expansion of the integral (5.1C) about the saddle 
point k*(v) can be performed in a frame with arbitrary 
velocity v and yields 

G(x,t) = e lk * i+ (- tuj{k ' )+lvk ') t I v (£,t) , i = x-vt. 

(5.14) 

The integral T v (£,t) is expressed after substitution of 
(k — k*) — Kj\ft as 



T — 



dn 

e -e/(4Dt) 



jK,Z/Vt-DK 2 +0(D 3 K 3 /Vt) 



l + O 



DH 



(5.15) 



for large t and arbitrary £. Obviously D plays the r ole of 
a diffusion coefficient. D% is defined below in ( 5.28 ). 

Generically, the growth (or decay) rate of the saddle 
point mode Re (— iu)(k*(v)) + ivk*(v)) will be nonvan- 
ishing. We now define the particular linear spreading or 
pulled velocity v* through Re (— iui(k*) + iv*k*) = 0, or 
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Im w{k*) _ s(-ik*) 



Im k* 



Im k* 



k* = k*(v*). (5.16) 



This means, that in the frame moving wit h velo city v*, 
the absolute value of the Green's function (5.14) neither 
grows nor decays in leading order, v* , k* and oj(k*) are 
determined by ( 5.12| ), ( |5.16| ), and by Eq. (5.11) evaluated 
at v*: 



du(k) 



dk 



(5.17) 



In addition, the solution determines D = D (v*) ( p.l4[). 

Note that the leading order larg e t result (5.14)7 f r >.15| ) 
f or th e Green's function G in (5.9) is diffusive just like in 
( |2.45 ), despite the fact, that we are dealing here with an 
equati on w ith highe r spatial derivatives. We shall see in 
Sects. VD and VE that this even remains true for much 



more general types of equations. 

Note also that in our discussion of p.d.e.'s in this 
and the next Section, we o nly ta ke the spatial Fourier 
transform of G(x,t), as in ( |5.10 ) above. However, the 
most general formulation, which also applies to differ- 
ence equations or integro-differential equations, is most 
conveniently done by taking a Fourier transform in space 
and a Laplace transform in time. In the present context, 
the Green's function G(k,u>) is then defined as 



oo />oo 



1 



G(k,oj) = dt dx e - ikM G(x,t) = 

Jo J-oo b{k,u) 

where S(k,ui) = iuu{k) — iu> , (5.18) 

and the long time asymptotics is determined by the dou- 
ble roots of the characteristic equation S(k,u>(k)) = 
(5.6). We defer this type of formulation, which closer 
follo ws the "pinch point" analysis of |5^,Q, to Section 

In practice, one first will drop conditio n ([5.12 ) and 



generically derive N solutions (k*,v*) from ( 5.16| ), (5.17) 
for a given dispersion relation. But as we already pointed 
out above, not all of these may be appropriate saddle 
points for the dynamics. Typically there are solutions 
with A* = Imfc* > and v* > 0, which describe a pro- 
file spreading to the right and solutions with A* < and 
v* < describing the spreading to the left, and the k- 
contour will have to be deformed through the appropriate 
one for the left- and right-moving front. These solutions 
are related by symmetry, if the ori ginal p.d.e. is sym- 
metric under space reflection: if (5.6) only contains even 



powers of k and if the a n are real, then for every solution 
(k* ,v*) there is a solution (— k* , —v*). Moreover, as men- 
tioned already above, there might be various nontrivial 
saddle point solutions which are not related by symmetry, 
if the degree N of spatial derivatives is sufficiently large. 
The saddle point analysis as well as the arguments of Sec- 
tion II E 1 for the competition between different solutions 
of the linearized equations clearly show that then the dy- 
namically relevant solution is the one with the largest 
velocity v* through which the fc-contour can be led. 

However, choosing the saddle point with the largest v* 
might according to counting arguments (as in Appendix 
H) be inconsistent with assumption B from Sect. VB of 
the existence of a family of uniformly translating fronts, 
since one expects the multiplicity of front solutions to be 
different for every saddle poi nt (v* , k*). The discussion of 
this issue we defer to S ection VI D, as for the applications 



discussed in Sect. VF, this problem does not rise. 



2. Uniformly translating solutions 

In the analysis of the nonlinear diffusion equation in 
Sect. ||, we saw that the uniformly translating solution 
decayed ase"^ with A real for v > v* . Here, 
A = \-(v) Q2.18D is the smallest root of v = s(A)/A, 
where s(A) (|5.7| ) here equals s(A) = A 2 + 1. v > v* 
implied Re k = 0, A = Im k > 0. These front solu- 
tions were found to be stable to perturbations which are 
steeper than the front solution itself provided there 
is no pushed front solution v c = V'. The solutions with 
v < v* had Re k ^ 0, Re u ^ 0, and were unstable. 

We will focus here on the immediate generalization of 
these results, i.e., assume that fronts with v > v* have 
Re k — , so that their asymptotic spatial decay is as 



,-AS 



In particular, this gives for the pulled fronts 

Re uj(k*) = , 
, s(A*) = Im uj(k*) > , (5.19) 



> 0, Im D = 



With this assumption we consider only the generic case, 
that dynamically accessible uniformly translating solu- 
tions of real equations will be characterized by a real spa- 
tial decay rate A and a real growth rate s, and that they 
will leave a homogeneous state (/) — 1 behind. This might 



Re k* 


= , 


A* 


= Im k* 




1 d 2 s 


D 




~ 2dJ? 



Elsewhere |75|, we will discuss an extension of the notion 
of uniformly translating fronts that allows to write pattern 
forming fronts in the Swift-Hohenberg equation as uniformly 
translating solutions of a suitable set of complex amplitude- 
like modes. For these Re k 7^ 0. Similar considerations hold 
for fronts in the complex Ginzburg-Landau equation itself 




41 



exclude some pathological cases of uniformly translating 
front solutions, that are not c harac terized by a real A 8 . 

If the sadd le point obeys (5.19), the expression for 
t ^> 1 ( 5.14 ) and v — v* for the Green's function G 
reduces to 



„-£7(4£>t) 



V4ttL>< 



1 + 



/ flag 



(5.20) 



Except for a rescaling of time and length scales with the 
real cons tants A* and D, this is precisely the functional 
form of (p~45| ). 

If we consider the velocity v(X) of the family of front 
solutions whose asymptotic spatial decay is as e~ x ^ with 
real A, then it is straight forward to see that $* is the 
slowest of all these uniformly translating fronts: Accord- 
ing to the linearized equation fl5.E| ) the solution in the 
leading edge is as e - Xx - w { lX ) 1 , The resulting velocity v 



v(X) 



—iu){i\) s(A) 



for all A . 



A A 
The minimum of this curve is given by 







dv(X) 



d 2 v{X) 



dX 2 



1 

A 

2D 
A*" 



ds(X) s(A) 



A 



> 



(5.22) 
(5.23) 



Ta king into account that w(k) is anal ytic, Eqs. ( 5.21) 
and ( |5.22 ) are equivalent to Eqs. ( [5.16 ) and ( 5.17 ), be- 
cause at a saddle of an analytic function, the maximum 
as a function of real k coincides with a minimum as a 
function of imaginary k. 

The analysis of the stability of the uniformly trans- 
lating solutions proceeds largely as in Section |II C| : the 
existence of a family of front solutions implies, accord- 
ing to counting arguments as given in Appendix [h], that 
there is at least a continuous spectrum of eigenmodes of 
the stability operator. Indeed, if we again write the tem- 
poral behavior of the stability eigenmodes as e~ <Tt and 
the steepness of the modes as A, and if we first focus on 
the spectrum of perturbations that is also continuous in 
A, then we have for the front solutions with v{X) > v*: 
a = — (s(A) — w(A)A). Expanding the A of the perturba- 
tion about the A of the front, we then get 



a(A) 



-A 



ds(X) 

dX 
dv{X) 
dX 



-v(X) (A -A) 



(A -A) 



(5.24) 



using (5.21) in the second line. Since we showed above 
that dvfdX < for A < A* (v > v*), a{A) > for 
A > A. This generalizes the result ( p!4| ) for the non- 
linear diffusion equation that the front solutions are 



stable to modes from the continuous spectrum which are 
steeper than the front itself. In addition to the continu- 
ous A spectrum there again may be discrete perturbation 
modes associated with the existence of pushed front so- 
lutions. 

We show in Appendix || that the existence of a strongly 
heteroclinic front solution $^ implies the existence of un- 
stable strongly heteroclinic stability modes for v < v\ 
again in parallel to the results for the nonlinear diffusion 
equation. The central assumption of our further analysis 
is of course that we are in the pulled regime, and hence 
that such solutions are absent. 

We finally note that the fact that v(X) has a min- 
imum for A = X*,v = v* , implies that for v < v* 
front solutions decay to zero in an oscillatory manner 
for £ — > oo as they have Refc ^ 0. By expanding the 
function v(X) about the bifurcation point at v* , A* , it is 
easy to show that for small \v — v*\, this branch of solu- 
tions has Im(fc-fc*) = A- A* w (X*) 2 v"' /(l2D)\v - v*\, 

Rek w y/X*\v-v*\/D, where v'" = ^j^U*. One usu- 
ally has v'" < and then such solutions are unstable 



(5.21) according to a slight generalization of (5.24) 



3. The leading edge representation 

As in our analysis of the pulled dynamics of the nonlin- 
ear diffusion equation, we will find it expedient to study 
the large time asymptotics in the leading edge by using 
the leading edge representation ip. For uniformly trans- 
lating fronts, the immediate generalization of the trans- 
formation ( [2.5C| ) from Section || is 



W,t) = <f>(x,t) e 



i = x -v*t 



(5.25) 



The linearized dynamical evolution equat ion fo r the lead- 
ing edge representation now generalizes (2.49) to 



d t rb = Vrb 



2 e -A*€ 



where 



V = e A * ? £„»(oo) e 



n=2 



D n d n e 



(5.26) 



(5.27) 



A short calculation (Appendix Q) reveals that the con- 
stants D n can be expressed in terms of the dispersion 
relation u(k) (|5.6|) as 



D,, 



1 d n 
n\ d(-X) 
1 d n 



d(-X) n 



n ( - iuJ ^ ~ 

s(X)-v*X 



v*X 



A=A* 



(5.28) 



Note, that in this generalized leading edge representation 



Q5.26 ) the coefficients of ij) and d^rb again are vanishing. 
This is an immediate consequence of the proper cho ice o f 
v* and A*. In fact, for uniformly translating fronts ( 5.19| ) 
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Dq = i s equivalent to the proper choice of the velocity 
v* ( 5.16 ) and Di = is equivalent to the saddle point 
equation (5.17) fixi ng A* for given v* . D2 is obviously 
identical to D from ( 5.13 ). We will see below, that in the 
leading edge, the contribution proportional to D2 = D 
gives the dominant contribution, while D3 appears only 
in the sub domi nant term, similar to what we already ob- 
served in ( 5.15|) . We therefore will essentially rec over the 
results of the nonlinear diffusion equation (1.1), which 
had the particular property of D n = for n > 2. 



4- The relaxation analysis 



this as a matching procedure.) In this analysis, the 
fact that the parameter a ^ in the asymptotics 
$*(£y ) = (a^x +/3)e~ A ^ x is nonzero (see Section 
[I E 2 ) plays a central role. 



The generalization of these steps to our equation (5.1) 
which is of higher order in space, is actually quite 
straightforward. We again use the general coordinate 
£x ( 5.30| ) with X(t) to be determined. Th e inte rior ex- 
pansion r){(;xii) = 4>— &*(£,x) from Section IIIB applies 
literally, except that we no w ne ed to use the linear op- 
erator C* = Cy (£ x) from fl5.3|). Accordingly, also the 
resummation ( |3.31 ) again is valid, and we again have 



We have now laid the ground work for the extension of 
the analysis of the rela xation of pulled fronts for our more 
general equation ( |5.lD in the case of sufficiently steep ini- 
tial conditions which as before are characterized by the 
requirement that 



lim (f>(x,0) e Xx = 



(5.29) 



for some A > A* . The an alysi s in Section [II for the non- 
linear diffusion equation (1.1) was based on the following 
steps: 

Step 1. 

The proper choice of the comoving coordinate sys- 
tem 

U=x-v*t~X{t) , X = ^ + ^ + ... , 

(5.30) 

allowing for a logarithmic shift X(t) oc lni in com- 
parison with to the asymptotic coordinate system 
f = x- v*t. 

Step 2. 

An expansion of <f> in the nonlinear interior part of 
the front about the asymptotic front profile $*(£x)) 
taken, however, not in the frame moving with ve- 
locity v* , but in the frame £x = x — J dt'v(t') with 
velocity v(t) = v* + X(t). 

Step 3. 

A resummation of this expansion of <f) in the cross- 
over region towards the leading edge, where the new 
variable z — £x/(4i) is introduced for the region 
with ix > 0(y/i). 

Step 4- 

An analysis of the leading edge in variables z and t, 
where <f> now is linearized about the unstable state 
= 0, and not about <&*. The two boundary con- 
ditions that (f> crosses over to the functional form 
of Step 3 for z <C 1, and that <j> is steeper than $* 
for z> 1, now determine both the functional form 
of <f> and the constants c„/ 2 m X. (We can think of 



(5.31) 



with <f> t , a uniformly translating solution of (5.1) with ve- 
locity v. The correction 0(l/t 2 ) is again non- vanishing 
and non-universal, as it depends on the precise initial 
conditions. 

The expansion of the interior shape towards the lead- 
ing edge ( 3 .36) ) depends on both the differential operator 
C* for £ — > 00 and on the shape of the asymptotic front 
$* ( |3~32| ). Eq. ( |H2| ) is generalized to 



« + /?) e" 



(5.32) 



since the saddle point expansion in Section V C 2 implies 
that for a pulled front $*, two roots of the dispersion 
relation coincide. Generally, 



(5.33) 



since a calculation resulting in a generalization of (2.54) 
can be set up along similar lines: If there is a bounded 
uniformly translating solution < I > *(£), then upon going to 
the leading edge representation and integrating the equa- 
tion for ty* (£) once over £, we find that a can be expressed 
in terms of the spatial integral over the nonlinear terms. 

How does the leading edge develop under inclusion of 
the higher spatial deri vativ es? Fir st of all we observe that 
the large- ^-so lutions (5.2C) and (2.45) of the linearized 
equation ( |5.5| ) are in leading order identical up to rescal- 
ing. In other words, the saddle point approximation ren- 
ders the spreading around the asymptotic exponential 
solution diffusive. This suggests that the leading edge 
can be analyzed by the same type of similarity variables 
(z, t ) as i n (3.40). In fact, in our shifted coordinate frame 
£x (5.30) the leading edge representation is 



(x,t) = e- x *tx^ x ,t) , (5.34) 
d t ip = V^ + X{d i -\*)i) + o{e~ x " £ - x ) , (5.35) 



with the differential operator T> from Eq. (5.27). After a 
rescaling with 
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Cy = A*£x i r = D 2 X* 2 t , d n = ^"^ 

£> 2 A* 



y = 



AA* _ Ci C| 



7- 7-3/2 

n-l 



A>A*< 
C„ = c„A* ^D 2 A 

this equation takes the form 

/ N \ 



(5.36) 



<9 T ?A = U? + ^ <z„d£ U + y (a c - 1) </> , (5.37) 



n— 3 



that is the same as Eq. (3.38) except that there are now 
the higher derivatives <9". As we show explicitly in Ap- 
pendix [k|, the lead ing ed ge ca n be analyzed with the same 
ansatz as in ( 3.40 ) and ( 3.42 ), 



(5.38) 



G(z, t) = y/r g-i (z) + g (z) + 



and in rescaled variables, one gets 



Ct = 



-3 



Cs = 



3^ 



9=±{z) = 2a-Jz , (5-39) 



3a 



go(z) =/3(l-2z) + 3a(l + d 3 )z - 2ad 3 z* - — F 2 (z) 
+ &a s/¥z (l - M ( -y, |,z 



In these variables the result is identical with that for the 



nonlinear diffusion equation in Section III, except for the 
additional terms proportional to d 3 in go{z). In particu- 
lar, the velocity parameters C\ and C 3 / 2 and the leading 
order contribution g=i (z) are independent of the value of 



^3, just like the subdominant term (3 from (5.32) enters 
go(z) but not the other quantities. That for the prob- 
lem written in variables z and r, d 3 can only contribute 
in subleading order, is in fact immediately obvious after 
the transformation. It is surprising, however, that the 
subleading velocity coefficient C 3 / 2 is independent of the 
value of g?3. We will find it to be unchanged even for 
much more general equations. 

In terms of the unsealed variables, the universal alge- 
braic convergence of the velocity is given by 



v(t) 



2X*t 



(X*) 2 Dt 



(5.40) 



where v* and A* are determined by the saddle point equa- 
tions ( |5.16 ) and (5.17) tog ether with ( 5.12Q , and wher e 
the diffusion coefficient D (5.13) equals D 2 from (5.28). 
The central results of this analysis are summarized in 
Table II. 



D. Generalization to single p.d.e.'s of higher order in 

time 



We now proceed in two further steps of generalization. 
In the present Section we first discuss partial differential 
equations for a single field cf>(x,t), which include higher 
order temporal derivatives as well as mixed temporal and 
spatial derivatives. These are of the form 



F^, d x cj>, 
d t d x c 



, 



(5.41) 



generalizing (pTj) to M > 1. In Section VE, we then also 



deal with difference or integro-differential equations and 
coupled equations. 



The extension to equations of type (5.41) presents no 
conceptual difficulty — we will follow here a route that is 
the immediate generalization of the discussion in the pre- 
vious Section. The new elements in the discussion will 
be the fact that higher order temporal derivatives and 
mixed spatial and temporal derivatives are generated in 
the dynamical equation for the leading edge representa- 
tion rp, but as we shall see these turn out not to affect 
the expression for the velocity relaxation and for the re- 
laxation of t he sh ap e in t he interior front region. The 
notation in ( 5.48 ) - flS.54 ), which may strike the reader 
at first sight as unnecessarily heavy, prepares for the dis- 
cussion of ev en more general equations and sets of equa- 
tions in VE, where finding a proper scalar leading edge 



representation is less straight forward than here. 

If we linearize (5.41) about <f> — 0, we get an equation 
of the form 



M N 



£ ]T a mn dr fl£ c/>(x, t) + o{4> 2 ) = . (5.42) 



m— n— 



For solving the initial value problem in time, it is conve- 
nient to Fourier-transform in space 



(*,*)=/ — e ikx cf>(k,t) . 
-00 27r 



(5.43) 



Below we will use the superscript ~ to denote a quantity 
Fourier transformed in space. 



The Fourier transformation of (5.42) results in an 
o.d.e. of order M for every Fourier mode <p(k,t): 



M 



N 



]T A m {k)dF 4>{k,t) = , A m (k) = Y,a mn (iky 



m=0 



n=0 



(5.44) 



Obviously, we need M functions to specify the initial 
conditions. We write these as an M-dimensional vector: 



.0. 



M-l ' 



(5.45) 
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The equation of motion (5.44) can now be written in 
Fourier space as 



d t £(k, t) = -fjk) ■ ±{k, t) , 
with the M x M matrix 



(5.46) 



£(*) 



( -1 

0-1 




\ As. Ax. _m_ 

\ A M A m A m 



\ 





-1 

A M -i 



(5.47) 



For later use, we here already define the matrix 
J(fc,w) = A M (k) (z(k)-iul 



(5.48) 

which later will result from a Fourier-Laplace transforma- 



tion as in (5.1q). Here and below, we use the superscript ' 



to denote a Fourier-Laplace transformed quantity to dis- 
tinguish it from spatially Fourier transformed quantities, 
which arc indicated with a tilde. 

The M eigenvalues u> m (k) (m = 1, . . . , M) of the ma- 
trix T_(k) are determined by the characteristic equation 
S(k, cj m (k)) = 0, where S(k,u>) is the characteristic poly- 
nomial 

M 

S(k, uj) — det §_(k, uS) = ^ A m (k) {~iuj) m 

M N 

= E E a ™« (-*")"* ( ?fc )" • ( 5 - 49 ) 



m— n— 



Defining the eigenvectors U_ m (k) of the matrix T_{k) 
through 



Xik)-U m (k)=icj m (k)U m (k) , 
and their adjoints through 

U i m (k)-U n (k) = S mn , 
the matrix T(k) can be written as 

M 

±(k)=Y,^ k ^rn(k)xU\ n (k) , 



(5.50) 



(5.52) 



where x den otes the outer product. 

Now (5.46) is easily integrated in time and the Fourier 
t ran sform ation inverted. We find in generalization of 
(p|) and (|57T0t): 



#M)= / Q(x-y,t)-£(y,0) 
Jy 

M 



(5.53) 



m=l J 



{k)t U m {k)*U: m {k) ■ 

(5.54) 



Obviously, the quickest growing mode U_ m {k) — charac- 
terized now by Fourier mode k and branch of solutions 
m — ■ again will be determined by a saddle point (5.16), 
( |5.17|) . Even more so than in the case of a first order 
equation, we can in general have more than one saddle 
point, as each branch of the disperion relation can in prin- 
ciple have one or more saddle points (a trivial example 
for two coupled equations is discussed in Appendix [I]). 
Again, the relevant saddle point is the one through which 
the fc-contour can be deformed and which has the largest 
velocity v* in the comoving frame. The associated saddle 
point values are denoted as k* = D etc. As before, 



we assume uniform translation as in (5.19), so that k* 
and u>(k*) are purely imaginary. Suppose, that v* lies on 
the branch uji(k). We then find in the comoving frame 
£ = x — v*t for long times t: 



m,t) = e 



-e/(iot) _ 



U_i(k*) xU\(k*) + 



(5.55) 



in generalization of (5.2C). 

This result shows that in the long time limit the 
Green's function G projects onto the eigendirection 

U_i(k*). T he re sult (5.55) is not restricted to the ex- 



plicit form ( 5.47 ) of the matrix T, so it applies to sets of 
coupled p.rf.e .'s ju st as well as they also can be written 
in the form (5.46). Projection onto the eigendirection 
U_i(k*) then defines the scalar leading edge equation re- 
sulting from coupled p.d.e.'s. We will further exploit this 
property in the following section. 

In the present section we just use ( |5.55| ) to calculate v* 
and A*, and to demonstrate why the leading edge trans- 
formation catches the relevant dynamic s. Pr oceeding as 
in earlier Sections, the scalar equation ( |5.42 ) now trans- 
forms under the leading edge transformation with v* and 
A* to 



( 5 - 51 ) <f>(x,t) = e- A *« ip(Z,t) , £ = X-V*t, 

M N 



(5.56) 



0= J2 E fl ™« (dt-v'd t +v t \ t ) m (d i -\') n i> 

m— n— 
M M+N 

= E E b mn d?d^(^t) . (5.57) 



m— n— 



Just as the a mn from Eq. (5.42) can be written in terms 
of derivatives of the characteristic polynomial S(k, u)) 
( |5T9l ) as 



S(k,u) 



(5.58) 



(fe=o)=0) 



so the b mn can be written as derivatives as well, simi- 



lar to (5.28). In showing this, it simplifies the notation 
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to use coordinates expanded about the saddle point by 
introducing the variables 



SI = lu — v*k 
and by defining 



k — iX* 



S* (q, fi) = S(k* + q, uj* + v*q + SI) , 
k* = iX* , uj* = v*k* 



(5.59) 



(5.60) 



When we will later consider the Fourier-Laplace trans- 
form of t) in the frame £, the frequency in this frame 
will turn out to be SI and the wavenumer will turn out 



to be g, since e 



-iujt+ikx 



). Accordingly 



the long time-small gradient expansion ofip(£, t) will cor- 
respond to a small 17-small q expansio n. In deed, in line 
with this interpretation, inspection of ( [5.57 ) shows that 
the b mn are simply 



(id^r (-i) n (d k + v*d u y 



ml n\ 

(id n ) m {-id q ) 



S(k,w) 



S*(q,Sl) 



(k* ,v*k* ) 

(5.61) 



( q =n=o) 



We will discuss the precise correspondence between the 
formulation in terms of S and the dispersion relation 
u>i(k) below, and just note here that the saddle point 
equations that determine A* and v* are expressed by 



■JQO 



S* (0,0) = , &oi = -id q S*(q,n)\ q=sl =o = 



(5.62) 



After dividing the whole equation ( 5.57 ) by 610 and in- 
troducing the notations 



D,. 



bon 
ho 



w 



hi 
ho 



Tl 



ho 
ho 



etc. , (5.63) 



the terms with the lowest derivatives are 



(d t + n d 2 



D 2 SS - D 3 & 3 



(5.64) 



+ w d t d £ + . . . ) ^ + o(^ 2 e" A * ? ) - 



= . 



This is the leading edge equation in its most general form. 
Note that after the leading edge transformation, the coef- 
ficient w may be nonzero even if the coeffi cient an = of 
dtd x <j) in the original equation of motion ( |5.44 ) vanishes. 

To show the connection with our discussion of first 
order equations in earlier Sections, it is instructive to an- 
alyze the relation between S and the dispersion relation. 
The various branches uj m (k) or Sl rn (q) of the dispersion 
relation are defined implicitly through the roots of 



S(k,uj m {k)) = 



S*(q,Sl m (q)) = , (5.65) 



As before, let u>i(k) (Sli(q)) be the branch on which the 
saddle point determining v* lies. Upon differentiating 



( |5.65 [ ) onc e with respect to k or q and using Eqs. (5.61) 
and (|5.62 ), we get our familiar result 



duji(k) 



dk 



dfii(g) 



dq 



(5.66) 



9 =0 



Likewise, by differentiating (5.65) twice we get 



d 2 fti(9) 



dq 2 



d 2 u 1 {k) 



dk 2 



d 2 Q S(q,Sl) 



dnS(q,Sl) 



g=fii(0)=0 

(5.67) 



If we combine this with the expression D = —b^/bio, we 
recover our familiar expression 



D = 



d 2 S(q,Sl) 



2id n S(q, SI) 
id 2 coi(k) 



id 2 Q-i(q) 



g=Oi(0)=0 



2dq 2 



q=Q 



2dk 2 



(5.68) 



For the case of an equation which is of first order in 
time, one can easily check that our gener al exp ression 
for D n reduces to the one given before in ( ^28^ , D n = 
{-i/n\)d n uj/d(ik) n \ k ,. 



Before we discuss the consequences of ( |5.64| ), we note 
in passing that formally we could have p rocee ded directly 
from the linearized equati on of motion ( |5.42| ) to the lead- 
ing edge representation ( 5.57 ) and hence to ( 5.64 ), by 
choosing the two parameters v* and A* such that the 
two conditions £>oo = = feoi are obeyed. The detour 
from this straightforward transformation via the saddle 
point analysis was taken to bring out the physical origin 
of the transformation in this context and to show why 
one has to use the saddle point (v* , A*) with the largest 
v* through which the contour can be deformed. In ad- 
dition, it explicitly shows how a particular "direction" 
U_i(k*) of the vector field <j) corresponds to the slow lead- 
ing edge dynamics. We will see in the next section that 
for coupled equations there is some freedom in choosing 
the projection onto a scalar leading edge variable. 

Let us now a nalyz e the implications of the leading edge 
representation (5.64). First of all, we observe, that a uni- 
formly translating pulle d front $*(£) = e" A *^^*(^) still 



will have the form ( |5.32j ) = (a$ + /?), and that the 



argument for a ^ from Section II E 2 still does apply. 

Can the extra terms t% dfi/j, w dtd(j ) etc. change our 
relaxation prediction from Section VC? A short inspec- 
tion shows, t hat a fter rewr iting t he eq uation in variables 
z and t, cf. d5~36l ) - (|3|) and (^39j) , w d t d^ will be 
of the same subleading order in 1/yi as D3 while 
both the terms Tid 2 ip and D4 d^ip will be one order lower. 
Also, when rewriting the equation in the variable £x = 
x — v*t — X(t), higher temporal derivatives will create 
terms like X and X 2 from the exponential factor in the 
leading edge transformation <fi(£,x,t) — e~ A ^ x *li(£,x,t). 
Since these are of order 1/t 2 , they do not influence the 
leading and subleading terms. 
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We do not repeat the detailed calculation here, be- 
cause it completely follows the lines of th e ear lier one. 
One finds that the result again is given by ( [5.3S| ), except 
that the subleading go(z) picks up another polynomial 
contribution from w besides the one from D3, namely 



g (z) = g (z) ((5.39)) +2q wX* 



z 2 - 



(5.69) 



The uniform velocity relaxation is invariably 



v{t) = v* 



2X*t 



(X*) 2 Dt 



, (5.70) 



and the interior part of the front is again slaved to the 
tip like 



(M) = *«(*)(£*) + 



(5.71) 



So the predictions from Table II also apply to p.d.e.'s 



with higher temporal derivatives like (5.40), if the front 
is pulled. 

Thus we reach the important conclusion that the uni- 
versal power law convergence is not an artefact of the 
diffusion-type character of the nonlinear diffusion equa- 
tion: it holds generally in the pulled regime of uniformly 
translating fronts, because the expansion about the saddle 
point, which governs the dynamics of the leading edge 
representation ip, is essentially diffusive. 



E. Further generalizations 

We now complete the last step in our discussion, and 
show that our results hold much more generally: even if 
the original dynamical equation is not a p.d.e., the dy- 
namical equation for the appropriate leading edge vari- 
able tp is still the same diffusion type equation (5.64), 
and consequently, our results for the velocity and shape 
relaxation from Table II do apply. 

When we have a set of coupled equations, we can view 
the m as components of a vector field, using a notation as 
in (5.46) with a different matrix T_(k). The main com- 
plication we are facing in this case is that the leading 
edge dynamics then not only "selects" a velocity v* in 
the pulled regime, but also an associated eigendirection 
U. m (k) in this vector space — this eigendirection deter- 
mines the relative values of the various fields in the lead- 
ing edge of the front. The long time dynamics in the 
frame moving with the pulled velocity v* is then asso- 
ciated with a slow dynamics along this eigendirection, 
while the dynamics along the other eigendirections is ex- 
ponentially damped. The appropriate scalar leading edge 
variable -0 will then turn out to be nothing but the pro- 
jection of the dynamics along this slow direction. 

The second complication is that we now consider equa- 
tions whose temporal dependence is not necessarily of dif- 
ferential type d^: they may just as well be of difference 



type or contain memory kernels. To treat such equations, 
we also perform a Laplace transformation in time besides 
the Fourier transformation in space just as in (5.18) by 
defining 

poo 

<fi m (k, us) = / dt e iut 4> m (k,t) . (5.72) 
Jo 

We thus consider dynamical systems that after the 
Fourier-Laplace-transformation of the equations, lin- 
earized about the unstable state, are of the form 



M 



^ S nm {k,uj) (j> m (k,uj) = 2J H nrn (k) 4> m (k,t = 0) , 



n= 1,...,M 



(5.73) 



The terms on the right hand side generally arise upon 
partial integration of temporal derivative terms, when 
we take the Laplace transform. They contain the initial 
conditions. Before exploring the implications of (5.73), 



we first discuss in more detail the type of systems whose 
linear dynamical equations can be written in the above 
form. 

Sets of p.d.e.'s: Single or coupled p.d.e.'s can gen- 
erally be written in the matrix notation (j)t + 7Z(k) 



4>(k,f) — 0, Eq. (5.44), and after Laplace transformation 
immediately yield ( |5.73| ), with the matrices S_(k,uj) — 
A M {k) (Z(k) - as before in ( pl8l ), and f£(k) = 

AM(k) 1. The leading edge behavior of single p.d.e.'s, 



where the matrix T_{k) has the explicit form (5.47), was 
discussed in the previous section. For coupled p.d.e.'s, 
the derivation of a scalar leading edge equation is not 
as straight forward, and also leaves some freedom, as we 
discuss below and for an example in Appendix [l| Never- 
theless, we will see that the results summarized in Table 

II are robust, in that they do not depend on the particu- 
lar choice made. We discuss examples of single p.d.e.'s in 
Se cts. V F 1 and VF4, and an example of sets of p.d.e.'s 

III |VF 2| . Of course, if one has a p.d.e. for a single scalar 



field 4>, one can directly take the Fourier-Laplace trans- 
form without writing as a vector field. This yields a 
slight generalization of (5.73), the most important differ- 
ence being that H then also depends on u>. Our results 
can ob viously also be obtained via this route - see Section 
V E 2 for further details. 

Difference- differential equations: When we have dif- 
ference equations in space, the equations can also be re- 
duced to the above form — the only difference is that 
upon Fourier transformation in space, the /c-values can be 
restricted to lie in a finite interval (the "Brillouin zone" , 
in physics terminology). An example will be discussed 
in Section VF3. Likewise, when we analyze a dynami- 
cal equation with finite time difference, the Laplace in- 
tegral can be replaced by a sum over integer times but 
the "frequency" remains a continuous variable. The only 
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difference is that upon Laplace inversion, the integral is 
over a finite interval of u values. Examples of difference 
equations in both space a nd tim e, arising from numerical 
schemes, can be found in 



VF6 



Equations with memory or spatial kernels: If the equa- 
tion has memory and/or spatial kernels of the type 
J dx' f* dt'K(x -x',t- t')<j)(x' ,t') JllO| , pTll , then upon 
Fourier-Laplace transformation these just give rise to 
terms of the form K(k,uj)(j)(k,ui) in ( 5.73| ), as will be 

1 The 



VF5 



illustrated with a simple example in Section 
only difference with the case of p.d.e.'s from this point 
of view then is that the elements S mn are not polyno- 
mials in oj and k, but more general functions of these 
arguments. 



1. Long time asymptotics of the Green's function via a 
Fourier- Laplace transformation 

We now return to the problem of extracting the 
long time behavior of the dynamical equation (5.73) in 
Laplace-Fourier representation. In analogy with our ear- 
lier analysis of p.rf.e.'s, and following |38|,[30|, we intro- 
duce the Green's function 9 G(k,u>) of the linear equa- 
tions, defined by 



G(k,u) =S(k, ojy 1 



(5.74) 



S_ is the inverse of the matrix S_. Eq. ( |5.73 ) now im- 
mediately can be solved as 



<j>(k, w) = G(k, uj) ■ H(k) ■ 0(fc, t = 0) 



(5.75) 



We write th e eige nvectors and eigenvalues of S_ in analogy 
to (|5~50l) - (|5~52l) as 



g(k,uj) - U m (k,u) = u m (k,u) U m (k,uj) 
The determinant of S can now be written as 

M 

S(k,uj) — det S_(k,oj) = J^J u m (k,oj) . 

m— 1 

and the characteristic equation 

u m (k,u m (k)) = 



(5.76) 



(5.77) 



(5.78) 



determines the dispersion relation Lo m (k) of the mode 
with eigendirection t£ m (k,oj). Note that each eigenvalue 
u m (k, uj) may be a nonlinear function of k and uj. There- 
fore it can happen that the equation u m (k,uj) — spec- 
ifies more than one branch uj(k) of the dispersion rela- 
tion. For simplicity, we will not distinguish this possi- 
bility with our notation, but we stress, that our results 
are generally valid. For equations of the form (5.46), 
we can identify u m (k,uj) — A]\.[{k) {iuj m {k) — iuj) and 
U m (k,uj) = U m (k). 

Upon inverting the Fourier and the Laplace trans- 
formation, where the Laplace inversion requires a suffi- 
ciently large real 7, we now find for the Green's function 
in the comoving frame £ = x — vt: 

t) = J dy g{£ -v,t)- 1 dy' K(y - y ) ■ Hy, o) , 

- l7+oc duj [°° dk 



D ik^—i(uj—vk)t 



, ■ ■ \.i 27T J _ oc 2"7T 



u m (k,w) 



G(k,u) , 



(5.79) 



The exp ressio n for G(£, t) is the immediate generaliza- 
tion of (5.54). When we evaluate the Fourier-Laplace 
inversion ofG(£, t) in the long time limit, each term in 
the sum ( |5.79| ) can be evaluated by a so-called "pinch 
point" analysis J58|,|30| making use of expansions about 
zeroes of u m (k,Lu). We then need to deform not only the 
contour of ^-integration, as in the saddle point analysis 
in the previous sections, but also the contour of uj inte- 
gration. The pinch point analysis is based on first eval- 
uating the /c-integral, and then the resulting w-integral. 
Alternatively, we can extract the long time dynamics by 
first closing the w-contour, and then performing the k- 
integral. This last route is closer to the one of Sect. 
For a further discussion of both approaches and of 
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the global conditions that determine which of the saddle 
points or pinch points is dynamically relevant, we refer 
to Appendix [M|. As before, we use the * to denote the 
appropriate solution that satisfies these conditions. 

As always, there can in principle be several saddle 
point or pinch point solutions through which the inte- 
gration contour can be deformed, and if this happens, 
the relevant one is the one corresponding to the largest 
velocity v* . If we again write iti(fc,o>) for the eigenvalue 
on which this solution lies and as before use a superscript 



A different choice for the definition of the Green's func- 
tion is G(k,dj) = < jJ )~ 1 ' iLik), which avoi ds th e con- 
volution of the initial condition with H_(z) in (5.79), and 
also for equations of the form ( |5.44| ) leads to the easier ex- 
pression Gj k, uj) = (f(k) . The advantage of the 
choice (5.74) is that we consistently work with derivatives of 
S = det S. 



48 



* for functions which are wri tten in ter ms o f the trans- 
formed variable Q and q as in (^59|) and (^60|) , the saddle 
or pinch point equations assume their familiar form 



u\{k* ,uj*) = 
(dk + v*d UJ )ui(k,uj)\k*,u* 



ul(0,0) = 0, (5.80) 
= <J=^> d q u[(q, fi)| ,o = . 



Note that since S" is the product of all eigenvalues, cf. 
Eq. ( 5.77 ), these equations are eq uivalent to those given 
before in terms of S, Eq. ( |5.62| ). Likewise, we get for 
the long time asymptotics o f the Green's function the 
immediate generalization of ( |5.55j ), 



(5.81) 



e -e/(iDt) O^k^u*) x U.\(k*,uj*) 



V4vrL>i 



ia w ui(fc,w)|( fe » )n «) 



which is our usual Gaussian ex press ion again, with D 
given by its familiar expression (5.68). 

Our strategy in deriving the long-time front dynamics 
is always to use the long-time evaluation of the Green's 
function just to show how the pulled velocity v* and 
the dominant exponential behavior e~ x * emerge, and 
to motivate why the leading edge variables ip(£,t) have 
essentially slow diffusive dynamics. The analysis of the 
slow i/j dynamics and the matching to the front interior is 
most properly done by going back to the p.d.e.(s) for the 
spatio-temporal evolution of ip. Switching back to the 
space-time formulation for %p comes out most directly 
from Fourier-Laplace inversion of the small-q and small- 
n expansion of the '(/'-equation. Indeed, for ip(£,t) the 

approp riate Green's function is e x ^G(£, t) and accord- 
ing to Q5.79 ) we have 



(5.82) 



which confirms that Q and q are the proper Fourier- 
Laplace variables of the leading edge variables ip. 



2. The case of a single field 

In contrast to our earlier matrix notation, a single 
equation for a single field <f>(x, t) after Fourier-Laplace 
transformation can also be written in a scalar form: 



S(k, to) 4>{k, oj) 



(5.83) 



function jfc, 0(fe,O), dt<t>(k,t)\t=o, ■••} 



The most common and direct way to arrive at the above 
equation is by performing a Fourier-Laplace transforma- 
tion on the original dynamical equation. In this case 
one immediately gets the characteristic function S(k,u>) 
on the left hand side, while the partial integrations (or 
partial summations in the case of difference equations, 
where also the derivatives in the initial condition terms 



are replaced by finite difference versions) of higher order 
temporal derivatives yield w-dependent initial condition 
terms on the right in (5.83). Of course, we can also arrive 
at this equation via the route of Section VP , where we 
introduced a vector notation for a scalar p.d.e. of higher 
order in time, so that the dynamical equation is of matrix 
form ( |5.73| ). Indeed, when we then calculate det S_(k,u>) 

with S_(k,uj) = A\i(k) (x_{k) — iwlj by developing the 

determinant along the last row of the matrix, one easily 
sees that one just retrieves the above result. 

Of course, the asymptotic analy sis o f t) paral- 
lels the earlier discussion of Section VD| , irrespective of 
whether or not the equation is written in vector form. 
Again, the asymptotic spreading speed is given by a sad- 
dle point of S(k,u/). However, as we have seen, for ana- 
lyzing the proper front dynamics we want to return to the 
dynamical equation for the leading edge variable ip. For 
the case of a p.d.e., this can be done simply by transform- 
ing the original equation for <fi to the leading edge rep- 
resentation ip(£,t) = e A i), but for difference equa- 
tions or equations with memory terms, additional steps 
are clearly necessary. The general analysis is based on 
the observation that in the leading edge representation, 
the dynamical equation is of the form 

S*(q, f2) ip(q, ft) = initial condition terms. (5.84) 

If we expand S* in q and Q m and perform an inverse 
Fourier -Lap lace transform, we immediately arrive at the 
p.d.e. (5.57) for ip(£,t) with coefficients b mn given in 
terms of the derivatives of S* according to ( |5.6l| ) ! From 
there on, the analy sis c ompletely follows the one in the 
last part of Section V D , and we recover again all our fa- 



miliar expressions for the relaxation of the front velocity 
and the profile. 

We stress that for a given equation, the transformation 
to the leading edge variable can be done exactly. If this is 
done for a p.d.e., we again get a p.d.e. of finite order. As 
no approximations are made, the resulting equation still 
allows one to study the fast or small scale dynamics in 
the linear region as well. For finite difference equations 
or for integro-differential equations, the transformation 
to the leading edge variable ip still results in a finite dif- 
ference equation or an integro-differential equation: the 
usual p.d.e. for ip then only emerges if in addition a gra- 
dient expansion is made for ip. Such an expansion will 
obviously contain an infinite number o f term s . (We will 
see ex plicit examples of this in Sections VF5 , VF 5 , and 
V F q ) . Normally, such an expansion is not of much use. 
However, when we turn to the long time relaxation to- 
wards pulled fronts in the leading edge, ip becomes arbi- 
trarily smooth and slow and hence the derivatives become 
nicely ordered. Moreover, the long-time large-scale re- 
laxation of ip corresponds precisely to the low-frequency 
small- wavenumber behavior of the Fourier-Laplace trans- 
form and this is why the expansion of S* gives the proper 
evolu tion equation to analyze the front relaxation: As 
([5.61 ) shows, the coefficients b rnn in this equation are 
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then nothing but the expansion coefficients of the char- 
acteristic equation S*(q, f2) for small q and Q. In other 
words, independently of whether we started from a differ- 
ential, a difference or an integro- differential equation, we 
find at this point always the same p.d.e. for the leading 
edge variable tp, and hence the same expression for the 
velocity relaxation! 

Let us finally remark, that instead of the leading edge 
transformation, we could also have performed a leading 
edge projection onto the slow dynamics, as discussed in 
the following section. We will show, that the results of 
Table II do not depend on this choice. 



3. The case of a set of fields and possible projections 

For dynamical equations which inherently consist of 
sets of equations for more than one field, one obviously 
can only arrive at an equation for a scalar variable iji by 
some kind of projection onto the slow direction. The way 
in which one projects out the slow dynamics clearly en- 
tails a certain freedom of choice. For a given equation 
the "best" choice may be obvious, but in general there is 
some ambiguity. We illustrate this explicitly in Appendix 

B 

We note first, that a vector field i/j(q, can be decom- 
posed into its dynamical components 7r m (q, f2) as 



M 



(5.85) 



rn — 1 



.* t 



n m (q,n) = U m (q,n)-l(q,n) , (5.86) 

where the superscript * on the eigenvectors U_ m and 
eigenvalues u m is to remind us that these are written 
in terms of the variables q and Q. 



Each n m (q, f2) has its own dynamics, cf. (5.84), 



u* m (q, Q) TT m (q,Q) = initial condition terms m . (5.87) 

The natural projection onto a scalar leading edge vari- 
able is thus onto the eigendirection with the largest v*, 
which we denote with U_i(q,£i). We then identify the 
scalar leading edge variable with 7ri(q, f2). Inverting now 
the Fourier-Laplace transformation, we find a p.d.e. for 
7Ti(£, t) of the form (5.57) with the coefficients 



(5.88) 



( g =n=o) 



Defining the saddle point parameters just as in ( 5.62| ) 
and (5.63) for Eq. (5.57), they in general will depend 
on whether we derived the coefficients from S or from 
u\. However, we will argue below, that the saddle point 
parameters v* , X* and D do not depend on this choice. 

Though the projection onto U_i(q, 0) is formally the 
simplest one, the direction of projection is actually not 
very practical, as it depends on q and f2. In practice, one 



will want to project along a fix ed di rection. Our previous 
analysis, summarized by Eq. ( 5.81 ), indeed suggested to 
project the long time dynamics of the Green's function 
onto U_ x {k* ,u)*) = C/i(0,0). Projection of ip(q,£l) onto 
this eigendirection yields 



r(q,n) = u* 1 \o,o)-±(q,n) 



(5.89) 



X^mM^! (o,o)-u m (q,n) . 



Now only for q w w O, we have ip p (q,Q) w 7tx(q,£l), 
while for finite q and fi, also Tr m (q,Cl) with m > 1 will 
contribute. Inverting the Fourier-Laplace transform and 
working in the frame £ = x — v*t, we find the contribu- 
tions from TT m >i to decay exponentially in time. Such 
contributions we encountered already a number of times 
before, for the first time in Sect. HE. The more impor- 
tant contribution comes from the coefficient of 7Ti , which 

is IZiV) ) ' IZifofi) = 1 - 0(q,Q). These algebraic 
corrections in q and O actually modify the b mn for the 
projection ij> p (q, f2) in comparison to Q5.88 ), except for 
the diffusion coefficient D, as we will see below. 

Still other projections might be physically useful as 
is illustrated on the explicit example of Appendix |]. 
We now turn to the consequences of all these different 
choices. 



4- The freedom of projection and the universality of Tables 
II and V 



At first sight, the leading edge transformation or the 
different leading edge projections each determine their 
own saddle or pinch point e quati ons or expa nsion pa- 



rameters b mn , compare, e.g., (5.61) with (5.88) 



Nevertheless, the definition of the saddle or pinch point 
parameters v* , A* and D in Table V does not depend on 
the choice of the leading edge transformation or projec- 
tion, and hence the universal relaxation results for the 
velocity v(t) and the shape & v (t) W Table II are indepen- 
dent of these as well. 

For the saddle/pinch point equations of Table V this 
conclusion is based on two observations: (?) S(k, ui) con- 
tains Ui(k,u>) as a factor (5.77). The saddle point is 



determined by a double root in k of ui(k,u), which can 
be written as 



u\{q, 0) = 6$ (-in + Dq 2 + ...) 



(5.90) 



= ui(k, lu) = bH' (-i(u - v*k) + D{k - k*f + ...) 



^02 Mo witn 



D here obviously is defined as D 

from ( |5.8§| ). The root (5.90) fully determines the lowest 
derivatives of S = Y[ m u m at the saddle point q = = Q, 
- up to a constant prefactor, resulting from the other 
factors in S. (ii) The saddle point parameters are defined 
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by ho mogeneous equations ( |5.62 ) or ratios of derivatives 
( |5.63 ). So the prefactors depending on differentiation of 
either u\ or S will cancel in the equations that determine 
v* , X* and D. In particular, D defined by D = — 602/^10 

in ( |5.63 ) is identical with D = —602 /^10' nere an< ^ with 
other D's resulting from different projections. 

The subleading ter ms D 3 and w for the scalar lead- 
ing edge variable in (5.64), in contrast, do depend on 
the choice of projection. Henc e, as there always will be a 
leading edge equation of form (5.64), and as the universal 
results summarized in Table II do not depend on the val- 
ues of D3 or w, Table II is a universal result, independent 
of the particular projection chosen. The subleading con- 
tribution go(z) in the leading edge will always be solved 



as in (5.69), so it will not depend on initial conditions, 
but it will depend on the direction of projection through 
the parameters D3 and w. 

In conclusion, we reiterate that the relaxation results 
also apply to dynamical equations other than p.d.e.'s, 
because the dynamics of the leading edge representation 
tj) becomes arbitrarily slow and diffusive for long times. 
This allows one to do a gradient expansion in time and 
space for ip, even if the original equations are not p.d.e. 's! 
In this case the path of analysis via the Fourier-Laplace 
transformation and pinch point analysis is necessary. For 
equations, that are of differential form in time, Fourier 
transformation in space and saddle point analysis is suf- 
ficient. 



Definition of ui m (k) : S(k,oj m (k)) = 



Saddle point equations: 
S(k*,uj*) =0 

(d k + v*d w ) =0 



Comoving frame: 
Im (uj* — v*k*) =0 

Diffusion constant: 



D 



-i(d k + v*d UJ ) S 



(fc*,w*) 



(definition: A* = ik*) 

lu* = UJ m (k*) 

„ _ du> m (k) 



dk 



(fc* 



Im u) m (k*) 
Im k* 



D 



id 2 uj m (k) 

2 dk 2 



In general, only saddle points with Re D > are relevant 
In this paper only saddle points with D real are considered. 



Table V: The saddle or pinch point equations, determining v* , k* — iX* and D for a given characteristic 
function S(k,u>) = det S_(k,oj). If there are several saddle point solutions that satisfy the global conditions 
that determine which saddle point solution is dynamically relevant (see Appendix 0) , take the one with the 
largest v* . 



F. Applications 



In this subsection we support the above arguments by 
summarizing the results of numerical simulations of three 
equations — a spatially fourth order p.d.e., a set of two 
coupled p.d.eJs and a difference-differential equation — 
which are all in complete agreement with our predicted 
universal relaxation trajectory as i n Table II, consist- 
ing of the velocity convergence ( 5.70| ), the slaved interior 
Q5.71D , and the cross-over to a diffusive type of dynamics 
in the leading ege for £ > y/i. We also briefly consider a 



p.d.e. with second order temporal derivatives, an exten- 
sion of the nonlinear diffusion equation with a memory 
kernel, and the discretization corrections in the Euler and 
in the semi-implicit numerical integration method for a 
nonlinear diffusion equation. The last results were used 
already in Section IV in our numerical study of the non- 
linear diffusion equation. 
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1. The EFK equation 

The EFK ( "extended Fisher KolmogorofT" ) equation is 
an extension of the nonlinear diffusion equation f37| , |66f , 
which has been inv estig ated quite intensely in the math- 
ematical literature 108 1 . ft reads 



dr 



d 2 x cf> - + - ^ 



(5.91) 



A straightforwar d cal cu lation |65fl shows that the saddle 
point equations ( pM6 ), ( 5.17 ) and ( 5.13 ) yield 



2A* (1-2 7 A* 2 ) 



1 - VI - 12 7 
67 



D = y/1 - 12 7 



1/2 



for 7 < 



12 



(5.92) 



For 7 > 1/12, the saddle point solution has Re k* 7^ 0, 
and in agreement with this, the pulled fronts in this equa- 
tion are then found to be non-uniformly translating and 
to generate periodic patterns |6(| . We will therefore focus 
here on the regime 7 < 1/12. The arguments of the Ap- 
pendix of |67| for the multiplicity of front solutions (sum- 
marized in Appendix |H|) give evidence that this equation 
indeed admits a family of uniformly translating fronts 
in this regime. One also can prove that the front cannot 
propagate with a larger velocity than v* if the initial con- 
ditions are sufficiently steep j74|,(75) . The convergence to- 
wards the pulled front solution should therefore be given 



by Eq. j5lo|) for v(t) and Eq. (|5~3ll ) or Eqs. (|3|)-(|3|) 
for the interior of the leading edge of the front profile. 
Fig. 11 shows some of the results of our numerical sim- 
ulations for v(t) at 7 = 0.08. This value of 7 is closely 
below the bifurcation value j c = 1/12 = 0.083. The 
plot is of the same type as in Fig. 8(a) for the nonlinear 
diffusion equation. 

The numerical grid sizes of the simulation are Ax — 
0.01 = At. The system size is L = 200, the initial con- 
dition is characterized by \ n u = 20 and x$ = 25. The 
analytical prediction for 7 = 0.08 is according to ( |5.92 ) in 
the limit Ax -> 0, At -> 0: D = 0.2, A* = = 1.29, 

and v* = 4.4 ■ A*/3 = 1.89. The ratio between the l /t- 
and the l/t 3 / 2 -contribution in v(t) according to ( |5.40| ) is 
measured on the time scale 



T= 1/(A* 2 L>) 



(5.93) 



as in the dimensional analysis (5.36). For 7 = 0.08, we 
have T = 3. The plot of Fig. 8(a) gave good results from 
time t = 20 on, where T = 1. It is therefore consis- 
tent, that the plot of Fig. 11 with T = 3 is good from 
times t = 60 on. We thus plot here the time interval 
60 < t < 200. One can already anticipate from the plot 
that again a correction of v* for the numerical finite dif- 
ference code will be required if we proceed to even higher 
precision. In conclusion, we find the results to be in full 
accord with our analytical predictions. 
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FIG. 11. Velocity relaxation in the EFK equation (5.91) for 7 = 0.08: Plot of v^,(t) — v* — X 3 / 2 as a functi on o f A3/2 as in 
Fig. 8(a) for times 60 < t < 200. System size L — 200, front position xo — 25, initial steepness Xinit = 20 in (4.2). Grid sizes 
Ax = 0.01 = At. 



2. The streamer equations 



Streamers are discharge patterns which result from the 
competition between an electron avalanche formation due 



to impact ionization, and the screening of the electric 
field by charges. For planar streamer fronts, the equa- 
tion for the electron density a and electric field E are 

11 
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d t cj = D a d 2 x a + 3 x {oE) + af str (E) , 

d t E = -D a d x a - aE , (5.94) 

where we have assumed that in the region x 3> 1, where 
the electron density vanishes a + = a(x — > oo,i) = 0, 
the electric field E + = E(x — > oo,t) does not change 
in time: dtE + — 0. The field dependent ionization 
rate has a functional form like f s tr{E) = \E\ e~ x l\ E V 
This is the functional form we use in our simulations. 
The state (<7,E) = (0,E + ) is unstable, and also for 
these equations it is known fl5|| , that they admit a one 
parameter family of uniformly translating front solu- 
tions. The dispersion relation for linear perturbations 
about the unstable state a — 0, E = E + < reads 
—iuj(k) = ikE + + f s tr(E + ) — D a k 2 , where we choose to 
analyze the leading edge in a projection onto the tr-axis. 



The saddle point equations ( 5.16| ), ( |5.17| ) and ( 5.15 ) then 
yield 



= — _E H 



2^D a f str {E+) 



A* = Vfstr(E+)/D c 

D = D a . 



(5.95) 
(5.96) 



Again, the simulations of these equations show that the 
velocity convergence follows our analytical prediction 
( |5.40 ). An example of our results is shown in Fig. 12 
in a plot as in Figs. 8(a) and 11, where we track various 
level curves of the electron density a. The dimensionless 
time is T = l/f str (E+) = e 1 = 2.718 for E+ = -1. We 
plot our data for times 40 < t < 200, and again find our 
predictions to hold. 
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FIG. 12. Velocity relaxation in the streamer equations (5.94) for E = —1 and D = 0.1, plotted as in Figs. 8(a) and 11 for 
times 40 < t < 200. Initial condition: Gaussian electron density a(x,0) = 0.9 e~ x (thus \i n it = oo), E(x,0) = —1. System 
size L = 400, front position shifted back to xo = 100, after it is reached. Grid sizes: Ax — 0.01, At = 0.0025. 



3. A difference-differential equation 

We now summarize some key elements of our analysis 
[109 1 of the difference-differential equation 



d t Cj(t) = -Cj + Cf_! , C (t) = , C J>x (t) = 1 , 

(5.97) 

with j integer. This equation originates from kinetic the- 
ory [17j . If we transform with (t) = 1 — Cj (t) to 



0O(*)=1 . 03>l(*)=O, 



(5.98) 



we have our usual notation with the state <pj = be- 
ing unstable and the state (pj = 1 stable. As usual we 



consider fronts between these states starting from suffi- 
ciently steep initial conditions. It is easy to see t hat such 
initial conditions will create a pulled front [109]. 

Equation ( 5. 98] ) p rovid es the first illustration of our 
argument from Sect. VE that our analysis applies to dif- 
ference equations as well — with the only difference that 
the spatial Fourier modes k now extend over a finite in- 
terval or "Brillouin zone" < k < 2w only. Substitution 
of the Fourier ansatz ~ e -^t+ikj - m ^ ^ ne e q ua tion of 
motion linearized about the unstable state <f>j = 



d t <t>j = -<j)j + 2<t>j-i 
yields the dispersion relation 



iuj(k) = 2e 



- ik 



1 



(5.99) 

s(A) = 2e A - 1 . (5.100) 



As discussed before, the long time asymptote of the lead- 
ing edge is again determined by the saddle point which 
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obeys (5.12), ( 5.16| ) and (5.17). This results in 



9tV = Dd U ! 



(5.104) 



2e A 



2e 



A* 



A* 



(5.101) 



When we choose the solution with v* > 0, the saddle 
point equations are solved by 

— ik* = A* > real 
v* = 2e A * = 4.31107 , 

A* = (2e A * - l)/(2e A *) = 0.768039 , (5.102) 
D = D 2 = v*/2 , D n = {-) n v*/n\ . 



The D n are determined from ( HH ). We now perform 
the leading edge transformation 



/(*) 



-A*? 



(5.103) 



The large-time, small-gradient expansion in the leading 
edge now results in the p.d.e. 



The velocity convergence is again given by (5.40), with 
v* , A* and D given by ( 5.102[ ). We do find indeed that the 
fronts in this equation arc pulled, and that the velocity 
convergence follows Q5.40| ). This is illustrated in Fig. 13, 
where we plot (v(t) — v* + 3/(2A*i))/i~ 3 / 2 as a function 
of 1/y/t. v(t) = x{t) is the velocity of the front defined 
as x(t) = Yl°jLo4>j{k)- The curve in Fig. 13 should ex- 
trapolate to 3/(2A*) • \Jtt/(\* 2 D) = 3.0699 as 1/Vt -> 0. 
This predicted asymptote is marked by the cross on the 
axis. Indeed, the data of (v(t) - v* + 3/ '(2A*i)) /i~ 3 / 2 
for 40 < t < 4000 extrapolate very well to the pre- 
dicted asymptote, especially in view of the fact, that 
i 3/2 w 2 • 10 5 at the latest times. The slight offset at 
the end might be due either to finite system size L or to 
finite numerical discretization At. 
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FIG. 13. Velocity relaxation for the difference-differential equation (5.97), where v(t) = x(t), and x(t) = ^j'Wi see 



( p.98|) . Plotted is here (y(t) - v* + 3/(2A*t)) /t~ 3/2 as a function of 1/v^ for times 40 < t < 4000. The curve is predicted to 

r~ — 2 

extrapolate to C3/2 as 1/yt — > 0. The predicted value of C3/2 is marked by the cross on the axis. Initial condition 4>j(0) = e 3 . 



System size TV = 4000 grid points. Front shifted back to no = 75, after it has been reached. Temporal grid size At = 0.0005. 



4- Diffusion equation with second order time derivative 

Quite recently, it was shown ]8(J that, not surprisingly, 
fronts in a second order extension of the F-KPP equation, 



T-2 



d 2 (f> dcf) _ d 2 4 



(5.105) 



are also pulled. One interesting aspect of this equation 
is that while the diffusive spreading in a first order dif- 
fusion equation is, in a sense, infinitely fast, the second 
order term gives a finite speed of propagation of the dis- 
turbances. 



As discussed in Section |VD| , our results immediately 
apply to this equation, so t he ve locity and front relax- 
ation is then given by Eqs. (5.70) and (5.71), with 



vT + 4t 2 ' 
A* = VI + 4r 2 , 

D 



(5.106) 



(l + 4r 2 ) 2 



The expression for D nicely illustrates the efnctive renor- 
malization of the diffusion coefficient due to the second 
order time derivative. 



5. An extension of the F-KPP equation with a memory 
kernel 

As an example of an equation with a memory kernel, 
consider the extension of the F-KPP equation 
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d t (j>(x, t) = d 2 (f>{x, t) + / dt'K(t - t')4>(x, t') 
Jq 

~cj) k {x,t), (k>l). (5.107) 



Upon Fourier-Laplace transform ation as in ( 5.18| ), this 
equation is a scalar version of ( 5.75 ) with S(k,u>) — 
iu/ — k 2 + K (lu), and so according to our discussion of 
Section |VE| , our analysis directly applies. If we take for 
instance 



K{t-t') 



1 



-(t-t') 2 /4r| 



(5.108) 



the equation reduces to the F-KPP equation in the limit 
T3 — ► 0, and the characteristic equation becomes 



A 2 - a- 



2 2 



erfc(r 3 s) = , 



(5.109) 



where we follow the notation of Section V C 2 in writing 
s = Im u>, A = Im k, and where erfc is the complementary 
error fu nction. The results for v*, A* and D, obtained by 
solving ( |5.109 ) together with the saddle point condition 
ds/dX = s/A\\* numerically, are shown in Fig. 14. 

Other example s of equa tions with memory kernels can, 
e.g., be found in [110,111]. 




FIG . 14. Plot of v* , A*, and D as a function of T3 for the extension (5.107) of the F-KPP equation with a memory kernel 

(fjoi). 



6. Exact results for numerical finite difference schemes 

The fact that our results also apply to finite differ- 
ence equations has the important implication that if wc 
study a p.d.e. with pulled fronts numerically using a finite 
difference approximation with gridsize Ace and timestep 
At, we can calculate v*(Ax, At) as well as v(t; Ax, At) 
exactly. This allows us to estimate analytically the in- 
trinsic discretization error in these quantities, and hence 
to decide beforehand which grid and step size are needed 
to obtain a given accuracy. 

As a first illust ratio n, suppose, that one integrates the 
F-KPP equation (1.1) numerically with an explicit Euler 
scheme. This amounts to approximating the p.d.e. by 



Uj(t + At) - Uj(t) 



11 



At 



2uj{t) 



Uj (t)-u^(t). (5.110) 



si — Aa 



(Ax) 2 

Upon substitution of Uj(t) ~ e 
linea rized equation (we again follow the notation of Sec- 
tion V C 2 by writing s = Im uj, A = Im k), we obtain 



x = j Ax into the 



- 1 



At 



/sinhAAx/2 
Ax/2 



(5.111) 



which is straight forward to solve for s(A; Ax, At). As we 
emphasized above, by solving the saddle point condition 
ds/dX — s/X\\* = v* , we can obtain the exact values 
of v* , X* and D for any step and grid size, and in this 
way determine the accuracy of the numerical scheme. In 
general, these equations have to be solved by a simple nu- 
merical iteration routine, but for small Ax and At, the 
result can easily be calculated analytically: Expanding 
in Aa; and At, we find the dispersion relation 



s{X;Ax,At) = 1 + A 2 



For At -> 0, Ax -> 
result s(A) = 1 + A 2 , 
parameters, we find 



A 4 (A.x) 2 _ (l + A 2 ) 2 At 
12 2 



(5.112) 



0, this reduces to the continuum 
as it should. For the saddle point 
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v* = 2 - 2At + —(Ax) 2 

Euler: X* = 1 + At - l(Ax) 2 + 
8 

D = 1 - AAt + ^(Ax) 2 4 



(5.113) 



In practice, the Euler scheme is not used very often, 
because it is numerically very unstable and not very ac- 
curate. We have done all our simulations in Section IV 
and in the present section with a more stable and ac- 
curate semi-implicit method [112], which for the F-KPP 
equation amounts to the discretization 

Uj (t + At) - uj (t) 



At 



Uj + i(t) — 2uj(t) + Uj-i(t) 



(Ax) 2 

u j+1 (t + At) - 2uj(t + At) + Uj-i(t + At) 
(Ax) 2 



+ 



+ -[ Uj (t)+u j (t + At)}- 

- \ [2u)(t) + ku)-\t){ Uj {t + At) - Uj (t))] . (5.114) 

The term on the last line is obtained by expanding 
Uj (t + At) about Uj (t) to first order in uj (t + At) — uj (t) , 
so that one obtains a linear equation for the Uj(t + At). 
This expansion makes what would otherwise have been 
an implicit method, into a semi-implicit method. This 
difference, however, does not matter for the leading edge 
dynamics evaluated below. 

The dispersion relation is now given by 



tanhsAi/2 
At/2 



/sinhAAa;/2 
V Ax/2 



(5.115) 



which immediately yields s(A; Ace, At). For small Aa; and 
At the result is 



s(\;Ax,At) = 1 + A 



_> X 4 (Ax) 2 (l + X 2 ) 3 (At) 2 



12 



12 



(5.116) 



For this integration scheme, it is now straight forward to 
find 

v*=2+ 2 -(At) 2 + ^{Ax) 2 + --- 

2 1 

Semi-implicit: A* = 1 - -(At) 2 - -(Ax) 2 + ■■■ (5.117) 

3 8 

D = 1 + 3(At) 2 + ^(Ax) 2 + ■ ■ ■ 

We stress that these are the exact expressions for the ap- 
plication of this numerical scheme to the nonlinear dif- 
fusion equation, scaled to the normal form (2.1), ([2~^). 
They are therefore the "ideal" finite difference correction 
terms in the absence of numerical instabilities, round-off 
errors et ce tera. The correctness and accuracy of the 
prediction ( |5.117 ) for v* is demonstrated in Section IV 
in Fig. 8(6). 

We finally note that an early example of pulled front 
relaxation observed in a finite difference equation in space 
and time was seen in a mean-field model of ballistic 
growth J2^]. In this paper, the prefactor of the 1/t term, 
obtained by plotting v versus 1/t, was found to be about 
9 percent too small. Presumably, this discrepancy is due 
to the corrections from the 1/t 3 / 2 term: According to 

( |5.70D , the term (1 — J it/ (X* 2 Dt)) generally gives rise to 

a lowering of the effective slope in a v versus 1/t plot, as 
Fig. 7(6) clearly demonstrates. 
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VI. SUMMARY AND OUTLOOK 

A. Summary of the main results 

The essential result of this paper is that for front prop- 
agation into unstable states, starting from steep initial 
conditions, the convergence of front velocity and shape is 
given in the pulled regime by the universal expressions 



v(t) = v* + X(t) , 
3 



X(t) 



1 



2\*t V V {X*) 2 Dt 
1 



= *„ct)(ex)+oi ¥ 

£x=x~v*t- X(t) , 



-or 2 ) 

U < Vi , 



(6.1) 
(6.2) 

(6.3) 
(6.4) 



provided the asymptotic front profile is uniformly trans- 
lating. All terms in the expression for v(t), A*, v* and 
D, are given explicitly in terms of the dispersion rela- 
tion of dynamical equation, linearized about the unsta- 



ble state [see Eqs. (5.16), (5.17) and (5.19) or Table V] 



These results are also summarized in Table II. The depen- 
dence on pushing or pulling and on the initial conditions 
is sketched in Table IV. 

With universal we mean that not only the asymp- 
totic profile is unique, but also the relaxation towards 
it, provided we start with sufficiently steep initial con- 
ditions which decay exponentially faster than e~ A x for 
x — ► oo. Moreover, the relaxation is universal in that it 
is independent of the precise nonlinearities in the equa- 
tion, and of the precise form of the equation: It holds for 
p.d.e.'s, sets of p.d.e.'s, difference-differential equations, 
equations with memory kernels, etc., provided fronts are 
pulled and that the asymptotic front solution is uniformly 
translating, and provided that we are not at the bifurca- 
tion point from the pulled to the pushed regime, or at the 
bifurcation point D = towards pattern forming fronts 
(e.g., at 7 = 1/12 in the EFK-equation). The fact that 
the results also apply to finite difference equations has 
a nice practical consequence: If a p.d.e. is studied nu- 
merically using a finite difference approximation scheme, 
both v* and the prefactors of the algebraic relaxation 
terms can also be calculated exactly for the numerical 
scheme. This allows one to estimate in advance how big 
step and grid sizes need to be, in order to ach ieve a par- 
ticular numerical accuracy (see Section VF6 ). 

The remarkable relaxation properties are reminiscent 
of the universal corrections to scaling in critical phenom- 
ena, if we think of the relaxation as the approach to a 
unique fixed point in function space along a unique tra- 
jectory. An alternatively way to express this in more 
mathematical terms is to say that we have constructed 
the center manifold for front relaxation in the pulled 
regime. 

The above expressions contain all universal terms: 
those of order t~ 2 depend on the precise initial condi- 
tions and on the nonlinearities in the equations. The 



order of the limits is important here: Our results are the 
exact expressions in a \ jt expansion, i.e., when we take 
the large time limit while tracking the velocity of a par- 
ticular fixed value of <j>. To order 1/t 2 , this is equivalent 
to keeping ^x fixed. When we interchange the limits by 
taking t;x large at fixed time, there is a cross-over to a dif- 
ferent intermediate asymptotic regime for £,x > Vi- The 
different dynamical regions of a pulled front are sketched 
in Fig. 2. 

The slow algebraic convergence of pulled fronts to the 
asymptotic velocity has important consequences, as it 
prohibits the derivation of a standard moving boundary 
approximation for patterns in more than one dimension 
that consist of propagating pulled fronts whose width is 
much smaller than their radius of curvature |S7[ . 

While we have limited the analysis in this paper to 
equations that admit uniformly translating front solu- 
tions, it turns out that most elements of our analysis 
can be extended to pattern forming fronts for which 
Re k * and Re uj* ^ 0. In this case, the expression 
(|J) with 1/VD replaced by Re (X/s/D) applies 

In addition to our derivation of the above expressions 
for the convergence of pulled fronts, we have reformulated 
and extended the connection between front selection and 
the stability properties of fronts. This leads to an es- 
sentially complete picture also of front relaxation in the 
pushed regime and in the case of leading edge dominated 
dynamics resulting from flat initial conditions. For an in- 
terpretation of these results, again a consideration of the 
different dynamical regions of a front as in Fig. 2 is help- 
ful. The relaxation behavior in the pulled regime with 
sufficiently steep initial conditions can not be obtained 
simply from the properties of the stability operator of the 
pulled front solution, and therefore had to be obtained 
along a different route, which is summarized below. 



B. Summary of the main conceptual steps of the 
analysis 

The derivation of our central result on pulled front re- 
laxation is based on the following steps: 

1. From the dispersion relation uj{k) or from the char- 
acteristic function S(k,uj), we obtain v* , A* and D (see 
Table V). 

2. The double root condition which determines v* 
and A* implies that the asymptotic large £ behavior 
of uniformly translating front solutions is as <&*(£) = 
(a£ + P)e~ x where generically a^O. 

3. The double root condition which determines v* 
and A* also implies that the lowest order spatial deriva- 
tive term in the dynamical equation for the leading edge 
representation ip — e A ^ (/>(£, t) is o f the diff usion type, 
Dd 2 ip/d£ 2 (see Sections |VC 3| , |VD| and (Vj|). 

4- The diffusion type dynamics implied by 3. shows 
that in the co- moving frame £ = x — v*t, the front pro- 
file shifts back with the collective coordinate X(i) which 
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grows logarithmically in time. Linearization about the 
asymptotic front solutio n 3>*(£) in th e £ fram e is therefore 
impossible (see Sections [II A 3 and III A 4). Instead, we 
introduce the frame £x — x — v*t — X{t) with the expan- 
sion X{t) = ci/t + c 3 / 2 /t 3 / 2 + ■ ■ ■ and the corresponding 
leading edge transformation ip(£,x,t) — e x ^ x cf>(x,t). 

5. In the front interior, the long time expansion for 
X generates an expansion for the corrections to the front 
profile in inverse powers of t. To order t~ 2 temporal 
derivatives of the front corrections do not come in, so 
that to this order the equations for the profile shape re- 



duce to those for 3>„(t). This immediately leads to (6.3) 



for the time dependence of the front profile. 

6. In the leading edge, where nonlinearities can be 
neglected, we use an asymptotic expansion for ip(£x,t), 
linearized about ip = 0, in terms of functions of the sim- 
ilarity variable z = ^ x /(ADt) of the diffusion equation. 
Now for small values of z the expansion has to match the 
boundary condition rp w e~ A ^ x $*(£,x) ~ ct£,x + P (im- 
plied by observation 2.), and for large £x the terms in 
the (intermediate) asymptotic expansion have to decay 
as a Gaussian e~ z = e^^ x '^ iDt ' times a polynomial in 
the similarity variable z. These two requirements fix the 
constants c\, ■ ■ ■ in the expansion of X, and hence 

(ph. 



C. Open problems 

What one considers as remaining open problems con- 
cerning pulled front propagation, will depend largely on 
one's background and standards regarding the desired 
mathematical rigour. While our results are exact and 
yield an almost complete understanding of the general 
mechanism of pulled front propagation, they have, of 
course, not been derived rigorously. In physics, such 
a situation is often not just quite acceptable but even 
quite gratifying, but more mathematically inclined read- 
ers may wish to take up the challenge to provide a more 
rigorous justification. More work could also be done on 
enlarging the classes of equations for which the assump- 
tions underlying our approach can be shown to hold, i.e., 
for which one can show that fronts are pulled and that 
there exists a family of uniformly translating front solu- 
tions. 

Within the realm of our approach, one can consider 
slight extensions of our method to two nongeneric spe- 
cial cases. First of all, we have focussed on the case of 
sufficiently steep initial conditions such that the steep- 
ness A = — lim^—too In 4>(x, 0) is larger than A*. As we 
discussed at the end of Section [II , the intermediate case 
in which for large x <f>(x, 0) ~ x~ h ' e~ x x with v < 2, 
does give a ^-dependent result for the coefficient of the 
l/t term. According to Bramson ]77j], the next order 
correction is of order 1/ (t In t) . This suggests that for 
this special case logarithmic terms will have to be in- 
cluded in the expansion. Second, at the bifurcation point 



from uniformly translating solutions t o pat tern forming 
fronts, which in the EFK equation ( 5.91 ) happens at 
1/12, the diffusion coefficient D vanishes [see Eq. 
■92j) ]. At this bifurcation point, the equation for the 
leading edge representation ip is not of the diffusion type, 
so our asymptotic expansion breaks down right at this 
point. We have not investigated what happens then. 

As mentioned before, we will elsewhere address what 
we consider the most interesting remaining challenges, 
the extension of (part of) these results to pattern form- 
ing and chaotic fronts |75|,|76| and the question whether 
weakly curved fronts can be analyzed with a moving 
boundary approximation fl37| , an issue which is of cen- 
tral importance for understanding fronts in two and three 
dimensions like streamers 1L5|. 



D. The multiplicity of front solutions and of 
solutions of the saddle point equations 



As we discussed in Section VB , our general discus- 
sion of the convergence of pulled fronts to their asymp- 
totic velocity and shape is based on the assumption that 
a uni formly translating front solution < I ) *(£) exists (see 
(|5.19 ) for a definition), and that it is a member of a 
one-parameter family of front solutions. What happens 
if this family of front solutions does not exist has, to 
our knowledge, not been investigated systematically for 
real equations. However, experience with various pat- 
tern generating fronts — especially with a similar case 
in which no generalized uniformly translating solutions 
exist in the quintic complex Ginzburg-Landau equation 
p8[ , even though the dynamics is pulled — yields the sce- 
nario that the leading edge just spreads according to the 
linearized equations, and that the front interior "just fol- 
lows" , in the sense that if there are uniformly translating 
fronts solutions, the front interior and the region behind 
it relax smoothly, while if there are none, it is forced to 
follow the spreading in some other way. This leads one to 
conjecture that if there is no family of uniformly trans- 
lating front solution s, th e velo city relaxation will still be 
described by Eqs. fl6.l| ) and (6.2) in the leading edge, 
but that in the interior front region the dynamics will 
inherently time-dependent, e.g., incoherent J76| . 

This can occur in particular in the following situation: 
As mentioned in Section V C 1 , it can happen that the 



dispersion relation is such that there is more than one 
allowed nontrivial solution for the equations for v* and 
A*. According to the linearized equation, arbitrary suffi- 
ciently steep initial conditions will spread out asymptoti- 
cally with the largest of the speeds v* . Hence the asymp- 
totic spreading speed of pulled fronts emerging from steep 
initial conditions is simply the largest velocity v* . Now, 
according to a counting argument for the multiplicity of 
uniformly translating front solutions, the multiplicity of 
front solutions associated with different solutions of the 
saddle point equations for v* will differ: If there are two 
solutions v\ and with X\ < A|, the multiplicity of 
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front solutions with velocity near and an asymptotic 
spatial decay rate near A2 will be smaller than that of 
those with velocity near v\ and a spatial decay rate near 
\\ . Investigations of the issue of the competition between 
various solutions v* will therefore also bear on the issue 
raised in the beginning of this section, the question what 
happens when there is no uniformly translating solution 
<f>*. In particular, the dynamics in an equation that has 
a family of uniformly translating fronts associated with 
the solution t>J, should show a transition from smoothly 
relaxing interior dynamics for v*[ > v% to incoherent in- 
terior dynamics for < v^. 



E. A step by step guideline for applying these results 

If one just wants to apply our results to a given dy- 
namical equation with a given initial condition without 
worrying about the derivation and justification, one can 
simply follow the following guidelines: 

(i) Linearize the dynamical equation about the un- 
stable state, and determine the characteristic equation 
S(k,oj) — for modes e~ wt+ m the linearized equa- 
tion. 

(ii) Solve the double root or saddle point conditions 
from Table V to determine v* , k* and D. 

(Hi) Check whether the leading edge of the initial con- 
ditions is steeper than e~ A x with Im k* = A* . Only then 
the front is a candidate for pulling with an asymptotic 
velocity v* . 



F. The subtle role of the nonlinearities: an 
alternative intuitive explanation 



(iv) Check whether the conditions (5.19) under which 
fronts are expected to be uniformly translating, Re k* = 
Re uj* = 0, Im D = are satisfied. If not, the fronts will 
be pattern g enera ting rather than uniformly translating 
(see Section VI D| above) . 

(v) Assuming the conditions under (iv) are obeyed, so 
that the asymptotic front is expected to be uniformly 
translating, investigate by a counting argument or other- 
wise whether there is a one-parameter family of uniformly 
translating front profiles ^(C) that includes $*(£). 

(vi) Determine, by using bounds, comparison theorems 
or physical arguments, whether the fronts will be pushed 
or pulled. This determines, which particular regime from 
Table IV applies. 

(vii) If according to points (iv)-(vi) there is a family 
of front solutions that includes <&*, and if the dynamics 
is pulled, then our predictions (|6~l|)-(|6~3|) or Table II ap- 
ply. If the conditions under (iv) are satisfied but there 
is no family of uniformly tra nslating solutions according 
to (v), then our formula (|6.2[ ) should apply but one then 
expects intrinsic nontrivial dynamics in the front interior 
to remain, so that (6.3) does not app ly. If (iv) is not sat- 
isfied (as for the EFK equation (5.91) for 7 > 1/12), one 
expects pattern generating fronts with a similar algebraic 
convergence J75|,[76[. 



As we have seen in ( |2.46| ) and ( 5.20| ), the convergence 
of the linear spreading velocity to the asymptotic value 
v* is as v(t) = v* — 1/(2 X*t) + ■ ■ ■, while the convergence 
of nonlinear fronts is as v(t) — v* — 3/(2A*t) + The 
prefactor of the 1/t in the latter case is just three times 
larger than for the linear spreading velocity. What is this 
subtle difference due to? 

In this paper, we have attributed the difference to 
the presence of the term at; in the large-£ asymptotics 
(at; + P)e~ x 5 of $*(£)■ We used an argument closely 
relat ed to the one presented below, to prove in Sec- 
tion |IIE2| , that a ^ 0. The functional form of <f>* 
leads to the requirement that the leading term in the 
expansion in similarity solutions in the leading edge 
is (g/t 3 / 2 )e-£ 2 /( 4Df >, not (l/t 1 /2) e -? 2 /(4Ct) (gee SectiQn 
III A 1 ). Nevertheless, one may want to have a better 
intuitive understanding of why the asymptotics of the 
linear spreading velocity is not correct for the nonlinear 
front relaxation — after all, one might at first sight think 
that the linear spreading results should be correct suffi- 
ciently far into the leading edge where the nonlinearities 
can be neglected. The following picture allows us to un- 
derstand why this is wrong, and why the same type of 
algebraic convergence also applies to pattern forming and 
chaotic fronts [ ]75| , |76| |. 

Consider for simplicity the F-KPP equatio n (|lTl ). As 
discussed in the introduction and Section |II E 2 , the 
dynamical equation for the leading edge representation 
t) of is 



-2A*£ 



(6.5) 



We can think of the nonlinear ?/; 3 e~ 2A ^ term as a lo- 
calized sink term in the diffusion equation for ip: the 
term vanishes for positive £ due to the exponential term, 
and for large negative £ since tjj vanishes expone ntiall y 
in the region to the left where <fi saturates, see ( [2.54 ). 
Thus, if we think of ip as representing the density of dif- 
fusing particles, then in the region where this term is 
nonvanishing it describes the annihilation of particles. 
For the half space to the right of it, where the particles 
freely diffuse, this term therefore acts like an absorbing 
boundary on the left. This is actually all that remains 
of the nonlinearities in the equation! Whenever the inte- 
grated sink strength a [the s patial integral of the nonlin- 
ear term, in agreement with (2.54)] is nonzero, the prob- 
lem in the leading edge reduces to that of the buildup of 
a diffusion field in the presence o f an absorbing bound- 
ary (and at the same time, as ( 2.54| ) shows, a 7^ 0). 
In this language, the pulled to pushed transition occurs 
precisely when the absorption strength a vanishes, and 
indeed precisely at this point the veloc ity c onvergence is 
as v(t) = v* - l/(2A*t) + • • • [see Eq. (|j36|)]. 

There is one complication: unlike the usual problems 
of diffusion in the presence of a given absorbing bound- 
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ary, the "sink" in ( |6.5| ) depends on the relaxing field ip 
itself. In fact, as we discussed extensively in the paper, 
the diffusive dynamics of ip leads to a logarithmic shift 
of the sink in time, in the frame £. That is why in this 
interpretation we have to go, for selfconsistency, to the 
frame £x = £ — X(t). In this frame, the "sink" or "ab- 
sorbing wall" remains essentially fixed in time, and so 
the dynamics of ip is, in leading order, that of a diffu- 
sion field in the presence of a fixed absorbing wall. As 
is well known, in such linear gradient ip ex £x 

will build up in front of the wall, to balance the constant 
annihilation of particles in the wall region. 

Clearly, even if the "sink" strength is not stationary 
in time, the buildup of the linear diffusion gradient far 
ahead of it will not be affected. The present interpreta- 
tion therefore yields a natural starting point for analyz- 
ing the velocity relaxation of non-unifo rmly translating 
fronts. This will be explored elsewhere [109 75 7(|. 

We end this paper by stressing that while we have 
shown that nonlinear fronts relax according to the "3/2 
law" v(t) = v* — 3/(2A*t) + • • •, one can not apply this 
result completely with closed eyes. An amusing illustra- 
tion of this warning is the following. It has been noted, 
that the spreading velocity in the equation 

2 

(6.6) 



follows the "1/2 law" v(t) = v* - l/(2A*i)--< = 
2 — 1/ (2t) + - • • [ 1 13 1 . At first sight, this equation therefore 
might appear to yield a counterexample to our assertions. 
In fact, it does not. Our results only hold for equations 
where the growth of the dynamical field saturates be- 
hind the front, not in the case in which the growth is 
unbounded. If the growth is unbounded, our arguments 
for why a/0, and hence for the "3/2 law", break down. 
The above equation is precisely an example in which the 
growth does not saturate: For e > and positive <j>, the 
nonlinear term only increases the growth. Hence there 
is no saturation and the spreading velocity v n i(t) in the 
presence of the nonlinearities is larger than the one of 
the linear equation: v n i(t) > v* — 1 / (2t) H . Appar- 
ently, in practice the equality is obeyed asymptotically. 
Of course, if we add a saturation te rm of the type — <p k 
with k > 1 to the right hand side of (5.6), we obtain reg- 
ular fronts and our usual expression for v{t) is recovered. 
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APPENDIX A: AN UPPER BOUND FOR Vc IN 
THE NONLINEAR DIFFUSION EQUATION 



APPENDIX B: THE GENERALIZED 
NONLINEAR DIFFUSION EQUATION 



With a generalization of th e leading edge transforma- 
tion introduced in Sect. [IF, it is straight forward to 



prove the well-known upper bound v c 



< v sup , where 



2 sup 

0<<p<l 



(Al) 



for the selected front velocity in the nonlinear diffu- 
sion equation, if the initial conditions have steepness 
A > v sup /2. The steepness A of a fro nt i s defined in 



(2.6). To prove this bound, transform ( |2.l| ) to a frame 
£ = x — vt, and write 



(x,t) . 



The equation of motion is now 

2 

dtij) = dl<if> - V — i> + e v( - j2 f (V e-^ /: 

If the initial steepness is A > v/2, then 

lim ^(f , t) = for all < t < oo 



(A2) 



(A3) 



(A4) 



HE 



1 



since the steepness of the leading edge (£ — > oo) is con 
served for all finite times, cf. the discussion in Sect 
and since convergence at £ — > — oo is garantueed by c 
behind the front together with the transformation (|A2j). 
Thus the decay of ip at £ — > ±oo is exponential in £ for 
t < oo. Hence, the whole equation can be multiplied by 
ij) and integrated over £. This yields 



IT 

T 



-</2 



^ e -</2 



(A5) 



where all integrals are finite. The r.h.s. of this equation is 



strictly negative, if v > v sup (Al). Therefore, in a frame 



moving with velocity v > v sup , the integral ip 2 decays 
in time. This means, that the frame is propagating too 
rapidly, so that the front shrinks away in the leading edge 
representation ip ( |A2j ) . Only a frame moving with veloc- 
ity v < v sup can propagate along with the speed of the 



front. 



vsup 



is therefore an upper bound for the asymp- 



2. Hence, these fronts are pulled with 



totic velocity of any initial condition with A > v sup /2. 

For nonlinearity Jkpp = <P> — 4> k i we have v sup = 2. 
But on the other hand, we know (see Sect. |J), that 

Vsup > V c > V* 

v c = v* — 2. For nonlinearity f e = e 
we have v sup = y/1 + 4e > 2\fe = v* . 

This version of the argument for v < v sup J74[ can 
be generalized to equations with higher spatial deriva- 
tives, forming both uniformly translating fronts or pat- 
tern forming fronts Pq| . 



Analyse a general equation with first temporal and sec- 
ond spatial derivative: 

F (0, d x cj>, 3% d t cf>) = . (Bl) 

A front translating uniformly with velocity v solves 

F ($„, d^ v , 4$„, -vd^ v ) =0 , £ = z-vt. (B2) 

The stability analysis of such a solution and the further 
treatment of convergence is identical with what we did 
for the nonlinear diffusion equation (1.1) in Sections [0] 
and III. We only need to transform the linear operators 
as discussed below. Our ana lysis is directly relevant for 
the equation studied in [114]. 

We us e the definition of functional deriv ative s as in 
(|5.32j) - ( fc.36| ). A linear perturbation rj(g, t) ( fc.25| ) about 
a uniformly translati ng sta te then solves the linear 
equation d t r\ = L v r\ ( 2.26 ) resp. ( |5.30 ) with the linear 
operator being now 

Cv = MO d\ + f hv (0 d( + MO , fi, v = v + h. 

(B3) 

For transforming to a Schrodinger problem dfip = H v ip + 
o(ip 2 e~ a ), H v = —dy + V v (y), we now have to make the 
coefficient of the first order derivative vanish, and the 
coefficient of the second order derivative d 2 constant. 
Thi s ca n be achiev ed through a transformation similar 
to ( [j.2§| ) and (2.29), combined with a nonlinear transfor- 
mation y(£) of the length scale £: 



n v ( y ) 

dy(0 
V v (y(0) 



da{0 = 2 Jl*LZ M d {: 



4/ 2 



V v (y) 



«*■ d y = Vm)dt 



VMO 

flv - 4 /o/2 h djji^y - fi, v rfg/2 
4/2 2/ 2 
3 (d e f 2 ) 2 - 4/ 2 d\h 

+ I6/2 



(B4) 
(B5) 
(B6) 



(B7) 



We use again the convention lim^oo & v {0 = 0. By 
construction the pulled velocity v* is the velocity, where 
V v * (00) = 0. Accordingly now 



v* = 2^/0(00) Mpo) - /i(oo) 
The steepness of the leading edge is 



A* =d f 



£ — M30, V—V * 



I Moo) 
Moo) 



(B8) 



(B9) 



61 



(In the convention of Section VB: / n (oo) = c n .) In the 
leading edge region, the relation between y and £ is lin- 
ear: y = £/y/f 2 (oo). 

If T4« (y) > for all y, there are no d establizing linear 
modes within the Hilbert space of (B5). Then the front 
propagating with v* is stable. The remaining analysis 



translates from Sections |IJ and [II step by step with only 
the explicit form of the linear operators C v and TL V and 
the transformation operator e a being more involved. 

If there is a range of y such that V v » (y) becomes neg- 
ative, there might be a destabilizing mode in the spec- 
trum of linear perturbations. In this case, there must be 
a pushed front solution with some velocity > v* with 
steepness A = X + (v^) > Ao(v^) = v< /2. Such a pushed 
front might even be integrated analytically, if one can 
find an analytic solution of 



F 



vip = , 



(BIO) 



equivalent to (p3|). £(0) can then be integrated as in 
(C5). (Again, a closed form for ip(<fi) cannot be found 
for pulled fronts, except possibly for equations at the 
pushed/pulled transition.) 



APPENDIX C: ANALYTICAL SOLUTIONS FOR 
PUSHED NONLINEAR DIFFUSION FRONTS 
AND TRANSITION TO PULLING 

We here discuss, how to find analytical solutions for 
uniformly translating fronts 0(£) in the equation 



dl<j> + vd^ + f(<f>) = o 



(CI) 



• We rephras e and straighten the method from (see 
also [115, 1 16 1 ) how to find analytical front solutions. 

• We recall, that analytical solutions can be found only 
for pushed fronts (propagating either into a meta- or into 
an unstable state, Cases I and II from Table IV) , but not 
for pulled fronts (Case IV). 

• We recall, that only a strongly heteroclinic orbit, i.e., 
a front approaching = with A > Xq(v), is a candidate 
for a pushed front. This allows us to calculate the criti- 
cal e for the push ed/pulled transition in the case of the 
nonlinearity ( 1.1 0| ) 



Write the equation as a flow in phase space as in (2.24) 



( 



-vip - f{4>) 
4 



(C2) 



where £ parametrizes the flow. If is monotonic in £, ijj 
can be parametrized by instead of by £. This substitu- 
tion yields for ip(<fr) 



V>|| + ^ + /(0) = o 



(C3) 



This is the differential equation for the trajectory in 
phase space, where now the translational degree of free- 
dom is removed together with the parametrization £ of 



the flow. The resulting differential equation is one order 
lower than the original differential equation (CI). Ac- 
cording to ( 2.16] ) , the initial condition for the integration 
at « 1 is 



-0(0 = 1-5) = -\-5 + o(S 2 ) , 



(C4) 



A_ = v/2 - vV/4-/'(l) , 



so the front trajectory is unique and can be integrated. 
In some cases, the integration can be done analytically, 
if one i s luc ky enough to find an analytical solution ip((j>) 
of Eq. (C2) for a given /(0). If we have a solution ip (</>), 
then the function £(</>) can be integrated as 



0(0) V#) 



(C5) 



The final step consists in finding the inverse function 
= <H£) i ^ this is possible. 

Note now, that solutions ip((j>) can be found analyti- 
cally only, if approaches = with a single exponential 
oc e~ A ^, since only then ?/>(</>) has the simple analytic 
form tp((f) = — Xcj) + o{4> 2 ) - Any other form of the ap- 
proach to (j> = 0, cf. ( 2.17 ), would not be expressible in 
a simple analytic expression for ip(<f>). In particular, a 
generic $* front with $* oc (a£ + f3)e~ x *t in the lead- 
ing edge does not have a simple analytical expression for 
■0(0) since ^(0) = ~ A*0+a0/(a£+/3)+o(0 2 ), so a pulled 
front generically cannot be integrated analytically. Again 
ffi^O spoils the conventional tools of analysis! 

Given an analytical front solution with velocity v and 
decay rate A, one has to check the nature of the front. A 
pushed front is a strongly heteroclinic front, i.e., it has 
leading edge steepness A = X+( v) > X (v). (For the no- 
tation of A's, compare Eq. ( 2. 20| ) .) If A = Xq(v) = X±(v), 
we have found a front at the transition point from pushed 
to pulled with leading edge behavior oc e~ A ^. This is 
the only pulled front, we can integrate. If A = A_(u) < 
Xq(v), we have a particular flat front, that has evolved 
from an initial condition with the same flatness in the 
leading edge. 

Finding analytical solutions for pushed fronts can even 
be turned into a machinery, if we don't fix / and look for 
a ip, but if we define V>(0) and then calculate /(0). For 



i, = -A0(1 - 0") 



we calculate, e.g., 

/(0) = X(v - A)0 + A(A(n + 2)-v) 



(C6) 



in+l 



-X 2 (n+l) 



k n+l 



(1 



i2n+l 



t2n+l 



(C7) 



where we have to identify v = (n + 2) A — 1/A, and 
e = A(i? — A). Th e an alytic front solution for ( |C6| ) can be 
calculated from (C5) and inverted to yield 
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0(0 = [i + (0(0)- 



i 



-1/n 



(C8) 



This solution is a pushed front, if A > Ao(v) = v/2, 
which implies e < 1/n. For such e, we find pushed 
fronts with decay rate A = y/(e + l)/(n+ 1), velocity 
V = (1 + e(n + 2))/y(n + 1)(1 + e), and analytical form 
(C8). For e = 1/n, the solution is a front on the transi- 



tion point from pushed to pulled fronts with asymptotic 
decay (p oc e~ A *^ + o((f> 2 ). For e > the solution 

(C8) is a flat front evolving from flat initial conditions. 
Fronts evolving from sufficiently steep initial conditions 
are then pulled, propagate with velocity 2\/e, have decay 
rate A* = y/i and no analytic form of the front solution 
can be found. 

Eq. ( |l.l| ) with nonlinearity (C7) can also be rescaled 
to bring the equation to the more familiar form 



(C9) 



e = e(l + e) , (1 + e) = 1 = (^J = (jj . (CIO) 

This re prod uces precisely the form of the nonlin- 
earity ( 1.10| ) with the stable state now at tp s = 



[(l + y/l + 4e) /2] '". Accordingly the critical e is now 
e c = (n + l)/n 2 . Fronts propagate for e < e c with the 
pushed velocity w f = [(n + 2)\/l + 4e - n]/[2\/n + l] 
and decay rate A + (i; t ) = [l + y/l + 4e\/[2y/n + lj . For 
e > e c , they propagate with the pulled velocity v* = 2^/e 
and decay rate A* = yfe. 



APPENDIX D: LINEAR STABILITY ANALYSIS 
OF MOVING FRONT SOLUTIONS 

In this appendix we study the linear perturbations r\ 
of a uniformly translating front <&•«(£) in the nonlinear 
dif fusio n equation. The problem is defined in Eqs. ( 2.25| ) 
- (2.27) and can be transformed to a Schrodinger prob- 
lem ( p. 28 ) - (2.31) with linear operator TL V . Since Tt v is 
self-adjoint, we can decompose functions that lie in the 
Hilbert space of 7i v , into the orthonormal set of eigen- 
functions of TL V . Eigenfunctions in this Hilbert space 
form a complete set. However, it is obvious that not all 
linear perturbations with |n| <C 1 are in this space: Only 
perturbations with 

lim \n\ e Xo{vH < oo with X (v) = - (Dl) 

>oo 2 

can lie in the Hilbert space (which consists of square in- 
tegrable functions and of solutions proportional to plane 
waves e lk( - as £ — > ±oo). 



1. Schrodinger stability analysis 

The general properties of the spectrum and eigenfunc- 
tions of Tl v within the Hilbert space can be immediately 



obtained from a few well known results which to physi- 
cists are known from quantum mechanics (see, e.g., |p9fl ) ; 
since Ti. v is the Hamiltonian operator for a (quantum) 
wave in a potential in one dimension. The potential is 
asymptotically lower on the right than on the left, since 



V(oo) 



/'(0) < 



f'(l) = V(-oo), (D2) 



according to (|2.2| ). If we write the temporal behavior of 
an eigenfunction as ip a (0 e _<Tt , one finds that the spec- 
trum of 



(D3) 



is continuous for a > V(oo), and that the eigenfunctions 
are distributions, i.e., essentially plane waves e lk ^ with 
k = ±yer — V(oo) as ^ — > oo. One immediately con- 
cludes that a front <F„ with velocity v < v* = 2^//' (0) 
will be unstable against the continuous spectrum of linear 
perturbations with "energies" V^(oo) < a < 0. 

For a front with velocity v > v* — 2\J f (0), there 
still might be a point spectrum of bound and square in- 
tegrable states with a < 0. Bound states have a finite 
number of nodes, and there is a one-to-one correspon- 
dence between the number of nodes and the eigenvalue 
of the bound state "wavefunction" ?/v : the eigenfunction 
with the lowest eigenvalue a is nodeless (if it exists) , the 
eigenfunction corresponding to the next largest bound 
state eigenvalue has one node, etc. Therefore, the point 
spectrum is bounded from below by the "energy" a of 
the nodeless eigenfunction, if it exists. Now, one eigen- 
function is known: the translation mode f/'o clearly has 
cr = 0. It can be generated by an infinitesimal translation 
of$„: 



(D4) 



If is monotonic, ipo will be nodeless. If is non- 
monotonic, tpo will have nodes. 

From this one might be tempted to immediately 
draw conclusions on the stability of monotonic or non- 
monotonic front solutions. However, this is only po ssible 
if ipo is in the Hilbert space! Comparison with ( 2.17 ) 
shows, that this is the case, if either v — v* and a = 0, 
or if v > v* and A v = 0, i.e., for one of the strongly 
heteroclinic orbits. 

If a front obeys one of these conditions and if it is 
monotonic, then i/jq is the eigenfunction in the Hilbert 
space with the lowest "energy" a = 0. Therefore all 
other eigenfunctions will have a > and will decay in 
time as e~ at . An arbitrary linear perturbation in the 
Hilbert space can be decomposed into the complete set 
of eigenfunctions, and therefore it will decay too (apart 
of course from the nondecaying translation mode ipo)- 

If such a front is non-monotonic it will have n ex- 
trema, with n > some integer. The translation mode 
then has n nodes, and hence there are then n bound 
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eigenfunctions ip a with negative a. The front profile 
is then linearly unstable with respect to these modes. 
Since any generic initial condition will have a nonvanish- 
ing contribution from these destabilizing modes, a non- 
monotonic will generically not be approached for long 
times. Such a is called dynamically unstable. 

The analysis of the spectrum and eigenfunctions of the 
Schrodinger operator in the Hilbert space therefore yields 
the following results: 

1) A front <E>„ with velocity v < v* is intrinsically un- 
stable against a continuous band of linear perturbations 
from the Hilbert space. Such a front generically will not 
be approached under the dynamics. 

2) A front <&* with velocity v = v* and a — is un- 
stable against perturbations from the Hilbert space, if 
it is non-monotonic, and it is stable, if it is monotonic. 
There is a continuous band of linear perturbations with 
a > 0, that continuously extends down to a — 0. Accord- 
ingly, there is no gap in the excitation spectrum, which 
already hints at the non-exponential convergence towards 
a monotonic $*. 

3) A strongly heteroclinic orbit with v > v* and 
A v = 0, if it exists, is unstable against perturbations 
from the Hilbert space, if it is non-monotonic, and it is 
stable, if it is monotonic. If strongly heteroclinic orbits 
exist, by construction (see Section II B ) only the one with 
the largest velocity v — is monotonic, and the front 
$* with velocity v — v* < is non-monotonic and thus 
unstable. So only for <I>^, the spectrum of linear pertur- 
bations is purely positive: a > 0. For $^ there is at best 
a discrete spectrum of linear perturbations in the Hilbert 



space in the range < a < V(oo) = (v^ 



')/4, and 



the continuous spectrum begins at a > V(oo). Conver- 
gence of all perturbations in the Hilbert space will thus 
be exponential in time like e _<T *, with a the smallest pos- 
itive eigenvalue. 

Note the restrictions of this analysis: 

(a) Up to now, we have no predictions for fronts with 



velocity v > v*, whose translation mode ipo (D4) is out- 
side the Hilbert space. We will see that the equivalence 
of stability and monotonicity extends beyond the Hilbert 
space analysis. 

(6) The analysis of general initial conditions might re- 
quire linear perturbations, that lie outside the Hilbert 
space, even if ip is in the Hilbert space. 



2. Linear perturbations outside the Hilbert space 

The mapping to the Schrodinger problem is a power- 
ful method for perturbations 77 about a front $ v , that lie 
within the Hilbert space, because we then can work with 
a complete set of orthogonal functions. In general, how- 
ever, this space of perturbations needs to be completed 
by functions from outside the Hilbert space. 

To see this, consider for simplicity an initial condition, 
that is close to some $> v with v > v* , but steeper than 



this asymptotic front: lim x _>oo <fr(x, Q)/<& v (x) — 0. Then 
the steepness in the leading edge of r\ = <j> — will be 
dominated by and 



ip = r] e 



at; + (3 for v = v* 
e -A<w£ f or v = v i or generally 
for v > v* and A v = 



for v > v* and A v ^ 
(D5) 



with fi(v) — y/v 2 /4 — 1 > from ( [2.171 ). Accordingly, 



only for a pushed front propagating with velocity v 



(or more generally for a strongly heteroclinic orbit with 
v > v* and A v — 0) or for a pulled front with veloc- 
ity v* and a — 0, the linear perturbation 77 e^ '"^ is in 
the Hilbert space of Tt v . The decay of the zero mod e t/>o 
(D4) is asymptotically the same as that of ip in (D5). So 
a treatment of linear perturbations outside the Hilbert 
space is clearly called for. 

In general, we want to decompose perturbations r\ that 
obey 



lim |r?(C,t)|<l 

{— >±oo 



(D6) 



This is required for the linearization of <fi about in 
( |2. 26 ). We aim at a decomposition of r)(£,t) into eigen- 
functions ??a-(£) e~ at . We therefore return to the eigen- 
val ue eq uatio n for such an eigenmode, which according 



val ue eq uatio n tor sucn an eige 
to (|2.26|) and ( |2.27| ) is given by 



Vc 



= 



(D7) 



Our previous analysis in the Hilbert space already has 
identified m any of these eigenmodes, in fact all those, 
which obey (pl|). This criterium on 77^ is too strict at 
£ — > 00, so we now need to additionally analyze pertur- 
bations with e' Xo{v ^ < \7]a(£,)\ < 1 

as £ — > 00, which 

lie outside the Hilbert space. On the other hand, for 
£ — > — 00, Eq. (Dl) is less restrictive than (D6). This 
gives us the freedom to impose only |f? CT (£)| ^ e*°^m as 
£ — ► — 00, since such a divergence can be compensated for 
by perturbations from inside the Hilbert space, where we 
make use of its completeness. We therefore now impose 
the boundary conditions 



lim \r] a (£)\ < 00 

£ — >oo 



lim e A °^|7? CT (£)| < 00 , (D8) 

£ — > — 00 



where perturbations that additionally obey 
e _A °^^|77o-(£)| < 00 as £ — ► co, are in the Hilbert space 
oiH v . 

First o f all , we note that the translation mode rjo (£) = 
djj&viO (D4) now is always included in the larger space 
(Dq) of perturba tion s. 

Sec ond, solve ( p7|) for £ — > 00 and find in analogy to 
(§T7|) that 



-A_i 



-A+£ 



(D9) 
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with 



A ± (a,v) 



± 



/'(O) 



(DIO) 



For brevity of notation, we here allowed A±(a,v) to be 
complex. In Fig. 15 we plot A± versus a, both for 
the case of a front propagating into an unstable state 
(/'(0) > 0), and for the case of a front between a stable 
and a metastable state (/'(0) < 0), and for /'(0) > 0, 
we furthermore distinguish bet ween v > v* and v = v* . 
The leading edge solution (|D9|), of course, precisely co- 
incides with the leading edge behavior of the Hilbcrt 
space functions, except that one case was excluded from 
the Hilbert space: A leading edge with A a ^ and 
a < V(oo) = v 2 /4 — /'(0) does not obey the bound- 
ary condition (Dl). It does obey the boundary condition 
(|M), if a > -TlO). Let us therefore now focus on the 
additional perturbations with 



/'(0)<a<F(oc) = T -/'(0) 



(Dll) 



li A a ^ 0, such perturbations are outside the Hilbert 
space, but they do obey (D8). 

Are there such perturbations for a given a, and how 
many? For answering this question we need to analyze 
rja globally, in close analogy to th e glo bal analysi s of the 
as a function of v in Section |ll H , Solving (D7) at 
£ — > — oo yields two exponents 



A±(<r,v) = - ± 



/'(I) 



(D12) 



= \ {v)±y/V(-oo)-<r 



in analogy with (2.15) and (2.16). Since V (-co) > V(oo) 
(|D3]), for a < V(oo) we certainly have V(— oo) — a > 0. 

The coefficient of e~ A+ ( a ' v ^ therefore needs to vanish for 
rj a to obey (DS). Behind the front for £ — » — oo, we 



therefore find that 



Va(0 = ±e 



-A-(S-So) 



+ o e 



-2A_{ 



(D13) 



for an rj a obeying (D8) and (Dll). Eq. (D13) deter- 
mines r\ a uniquely because the arbitrary constant coeffi- 
cient ±e A -£° can be scaled out of a linear equation like 
(D7). Such a linear equation can always be integrated 
towards £ — > oo, where it u niquely determines the co- 
efficients A a and B a in (D9). Accordingly, A,, and B G 
generically are non- vanishing , in c omplete analogy to the 
argument for A v and B v in (2.17) to be generically non- 
vanishing in 

What do we gain with these extra solutions? The 
eigenfunctions in the Hilbert space had a continuous 



spectrum for a > V(oo) 



")/4 > and at best 



a discrete spectrum defined by A a = for a < V(oo). 
Adding the solutions, that obey ( |D8| ) , we extend the con- 
tinuous spectrum down to a > —F{0) = V(oo) — w 2 /4 < 



and find at best a discrete spectrum defined by A a = 
for a < — /'(0). These discrete solutions for a < — /'(0) 
all lie in the Hilbert space. 

Let us now look at the steepness in the leading edge 
of the solutions outsi de th e Hilbert space. They have a 
a from the interval ( Dll ), and A a ^ 0. For these we 
observe (cf. Fig. 15) that 

A_(cr, v) > A_(u) for a > (decaying) , 

A_(cr, v) < \-(v) for a < (destabilizing), (D14) 

A_ (0, v) = X- (u) for cr = (marginal), 



with \_(v) from fl2.16| ). This means, that these linear 
eigenmodes rja- of will decay (a > 0), if they are 
steeper than e _A -^^, and that they will destabilize a 
front if they are flatter. Note that the spectrum of 
decaying modes is continuous down to zero, as a j as 
A_ 4 A_. 

It is tempting to conclude here immediately, that a 
front $„(£) with velocity v > v* will be stable against all 
perturbations, which are steeper in the leading edge than 
e - A -(")£. However, the possible existence of the discrete 
set of solutions with A a = and a < requires spe- 
cial attention, since these perturbations are steeper than 
but destabilizing (er < 0). Now, if <E>„ is stro ngly 



-Ao(w)£ 



heteroclinic (A v — 0), we already found in Section D 1 
that such destabilizing perturbations exist, if and only if 
$ v is non-monotonic. We now need to show that this ar- 
gument also holds for fronts with v > v* and A v ^ 
or for fronts with velocity v* and a ^ 0. The fol- 



lowing five steps (i)-(v) prove this: (i) Impose (D13) at 
£ — ► -co. This defines a unique solution of equation (D7) 
for rj a for every a < V{— oo). In fact, we only need to an- 
alyze a < V(oo), sinc e we know the spectrum for larger 
a. (ii) Integrate (D7) forward towards £ — > oo for a very 
large negative a. The variation of f'(&v(0) m space then 
can be a lmost neglected. Therefore at £ — > oo, we will 
find (D9) with \A a /B a \ 3> 1. For our further construc- 
tion it is crucial to observe, that such a perturbation for 
sufficiently large negative a will be nodeless. It does not 
matter, on the other han d, t hat this solution typically 
will not obey our bound ( |D8| ), since we only use it as a 
means for construct ing the solutions with A a — 0, which 
will not only obey ( |D8[ ), but even lie inside the Hilbert 
space. (Hi) Upon increasing a continuously, at discrete 
values of a < V(oo), r\ a will gain an extra node. Since 
the generation of every new node is associated with a 
change of sign of the perturbation at £ — * oo, if the sign 
at £ — ► — oo is kept fixed, the appearance of an additional 
mode can only occur at a a, where the sign of A a changes. 
(iv) We know the number of nodes of the zero mode 770 • 
It is identical to the number of extrema of Q v . We there- 
fore know the number of particular perturbations with 
A a = and a < 0. (v) From this it follows that if 
is monotonic, there are no particular perturbations with 
Aa = and a < 0. If is non-monotonic, there are 
such perturbations. 
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In Section II B, we have counted the multiplicity of 
front solutions $ v as a function of v. Here we have 
counted the multiplicity of perturbations r\ a of a front 
<fr v as a function of a. This counting was based on the 
proper asymptotics of the solutions at £ — * ±00, which is 
of the same structure for both and r/a , so the counting 
argument follows exactly the same lines in both cases. 

The c oncl usions from this appendix are summarized in 
Section fl 



APPENDIX E: STABILITY ANALYSIS, 
SELECTION AND RATE OF CONVERGENCE 



In this appendix, we analyze the implications of the 
results of the stability analysis of Appendix [S] for under- 
standing the selection of fronts and for the rate of conver- 
gence towards the asymptotic front solution. For pushed 
fronts, the stability analysis implies that the relaxation 
towards the pushed front solution is exponentially fast, 
while for pulled fronts the spectrum is gapless, and the 
convergence can not be obtained from the stability spec- 
trum. 



1. Pushed regime: v c = 

We first consider equations with the nonlinearity f(cf>) 
such that the slowest stable front is a stron gly h etero- 
clinic orbit in phase space with A v \ = in (2.17). We 
have denoted this asymptotic front with <J>^ and its ve- 
loci ty wi th v'. Its steepness is X + (v^) = Xo(v') + fi(v'), 
cf. ( [2.23 ). There is a continuous family of stable front so- 
lutions <!>„ with velocity v > which are all flatter than 
A_(?/f) = Ao(w^) — ft(v^). Their steepness A = X-(v') is 
related to their velocity v through 



v(\) = \+j 



(El) 



as can be obtained by inverting ( 2.1§| ). 

CASE I: Consider an initial condition with steepness 
Knit > Ao(i^). We let <j> evolve so me t ime, and then 
linearize it about $t. According to (2.37) the perturba- 
tion 77 will have steepness A^ > Ao(f T ). It then is in the 



Hilbert space analyzed in Section D 1 . We can decompose 



the perturbation into the known eigenperturbations. The 
spectrum of decay rates has no negative eigenvalues, one 
eigenvalue zero and then a gap above zero. A contri- 
bution from the zero mode can be made vanishing by 
adjusting the position of the subtracted asymptotic front 
by making use of the translational freedom of <Ev 
The perturbation then can be decomposed into Hilbert 
space functions r\ a with a all positive and bounded away 
from zero. Thus, for large times the perturbation will de- 
cay exponentially. This means that an initial condition 
with Xinit > Ao(w^) will converge to <&t exponentially 



in time, generically with e _<Tlt , where oi is the smallest 
positive eigenvalue. 

CASE II: If the initial steepness is A_(t?t) < Xinit < 
Xo(v'), the perturbation of cf> about <&t will not be in the 
Hilbert space. However, we do know from the results il- 
lustrated in Fig. 15 that there is an eigenmode r\ a of the 
linear stability operator of $t with the proper steepness 
\nit = A,,, that will decay exponentially in time, see Sec- 
tion p 2j . The remaining linear perturbation r\ — r\ a might 
lie in the Hilbert space, in which case we are back to Case 
I. If it does not, we have to identify the subleading A, its 
corresponding eigenmode r\ a etc. The iteration of this 
construction leads us to conclude that the perturbation 
indeed will decay exponentially in time. (Examples of 
exponential convergence towards pushed fronts which is 
dominated by such modes can be found in Fig. 19 of p7|.) 
Another way of putting the argument is that only per- 
turbations with A < A_(t^) can grow in time, but these 
cannot be involved in the decomposition of a perturba- 
tion with A r) > X-(v'). A more elegant way of analyzing 
this case and the following ones is discussed in Section 
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CASE III: If the initial steepness is Xinit < A_(«t), 
and we linearize <j> about $t, there is a perturbation 
rj a with steepness A r) = Xi„it that is growing in time 
(er < 0). So such an initial condition cannot approach 
<Ff or any other asymptotic front with steepness 
^asympt > \nit ■ If we linearize (f> about the asymptotic 
front <!>„ with the same steepness Xinit = ^asympt, the 
remaining perturbation will be steeper, so contributions 
from the zero mode are excluded by construction, and the 
perturbation can be decomposed into eigenperturbations 
of <!>„, which all decay in time. 

In summary: All initial conditions with Xi n u > \—(v<) 
converge exponentially in time to the "selected" front 
with velocity v se i — tr and steepness A se ; = X+(v'). Ini- 
tial conditions with A, n jt < A_(i>^) will converge to a 
quicker asymptotic front with steepness X asymp t — A, n ,t 
and velocity v(X in i t) given by ([El] ). 

In Section [II A| , we have termed an initial condition 
sufficiently steep (Ai„# > A steep ), if it approached the 
"selected" front for large times. We have denoted the 
steepness of the selected front with A se ; . In the pushed 
regime, one can thus identify these parameters with 



it 



A steep = A„(w f ) = —- /j(7J f ) , 



(E2) 



Asei 

Vsel = V 



A+(«t) 
t 



+ M(w t ) , /i(«t) 



7jt2 _4 



2. Fronts into metastable states 

The only difference between a pushed front propagat- 
ing into an unstable state, i.e., with a nonlinearity / such 
that /'(0) > and v c — > v* , and a front propagating 
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into a metastable state, i.e., with /'(0) < 0, is the sign 
of A_(u): For a front into a metastable state, we have 



given by (El) and with steepness A a 



syrapt 



A; 



Suf- 



/i(v) 



liffl > | for /'(0) < , (E3) 



so \-(v) < and X + (v) > for all v > (the sign of 
Ao(u) is the same as the sign of v). Suppose, that the 
selected front still travels with positive speed v. 
(otherwise reverse x). Because now A_(u) < 0, 



scl 



= «t 



^steep 



(E4) 



so all initial conditions are sufficiently steep and converge 
to xhe continuous spectrum of asymptotic solutions 
with Xasympt < ^steep ceases to exist, and the asymp- 
totic front <E>t therefore now is unique. 

For the convergence of an initial condition (f> towards 
we still need to distinguish whether A,- n ,t is larger 
or smaller than Ao(i^) = /2. If Xinit > Xq(v>) the 
perturbation about $t lies in the Hilbert space, while 
for Xinit < \o(v<), it does not. This corresponds to the 
Cases I and II for v c = v< above, which apply literally. 
In both cases the initial conditions converge to $t expo- 
nentially in time. Case III does not occur for fronts into 
metastable states. 



3. Pulled regime: v c — v* 

At the transition from fronts propagating into 
metastable towards fronts into unstable states, /'(0) 
changes sign, and so does X—(v). At this point a contin- 
uum of possible attractors & v of the dynamics comes into 
existence, but the convergence behavior of sufficiently 
steep initial conditions is completely unchanged. In other 
words: Cases I and II are completely unchanged and only 
Case III needs to be considered additionally for initial 
conditions with X mit < X steep . 

A qualitative change in the convergence behavior of 
sufficiently steep initial conditions Ami* > X s teep only 
takes place at the transition from the pushed to the 
pulled regime. This happens for / changing such that 
ft approaches v*. Then 



A 



steep 



X (v*) = A se ; 



(E5) 



This transition leaves the multiplicity of possible attrac- 
tors unchanged, but the resulting changes in the spec- 
trum have deep consequences for the convergence behav- 
ior of sufficiently steep initial conditions. 

We now need to distinguish but two Cases for the ini- 
tial condition, namely A > A* and A < A* , where we use 
the short hand notation A* = X a (v*) = X±(v*) = v*/2. 

For flat initial conditions Xinit < A*, the arguments 
from Case III above apply literally. Such an initial condi- 
tion will approach a front $ t , with velocity v (Xinit) > v * 



ficiently steep initial conditions, however, exhibit a new 
behavior: 

CASE IV: Consider a sufficiently steep initial condi- 
tion with A > A* . As before we linearize the profile <f>(x, t) 
after a sufficient evolution time about the selected front 
$*. The corre spond ing perturbation rj = <f> — then 
decays like $* ( 2.21 ), because the steepness of <fi(x,t) re- 
mains larger than that of <&* at any finite time t, cf. Eq. 
( |2.37j ). As a result, r\ is ju st out side the Hilbert space in 
the generic case of a ^ fl2.17| ), just like the zero mode 
(|D4|). The Hilbert space has a continuous spectrum for 
all decay rates a > 0, and there are no growing pertur- 
bations with a < 0. The perturbation 77 can be written 
as a multiple of the zero mode rjo plus a remainder inside 
the Hilbert space. From this we might be tempted to 
argue that the perturbation will decay, and that we only 
can not tell how quickly — probably non-exponential, 
because the spectrum is gapless. However, in contrast 
to Cases I - III, there is no way to get rid of the zero 
mode, because no matter at which position £0 one places 
the subtracted $*(£ — £0), $* will always dominate the 
large £ behavior, and therefore the coefficient of the zero 
mode in the decomposition of the perturbation will al- 
ways be non-vanishing. A convergence argument based 
on simply neglecting the contribution from the zero mode 
is bound to be wrong: In the very same way we could ar- 
gue, that a steep initial condition converges to with 
just any v > v* . Strictly speaking, the linear stability 
analysis does not even allow us to conclude, that suffi- 
ciently steep initial conditions approach $* at all. Wc 
only can reason that there is no steeper attractor than 
$*, and that one therefore expects that the pulled front 
solution (f>* is selected from steep initial conditions. The 
diffe rent analytical tools that are developed in Section 
pTE] to analyze the convergence behavior, confirm this. 



APPENDIX F: GENERAL INTEGRATION OF 

Gn P /2(Z) 

We here show how to find special solutions 3*y 2 ( z ) °^ 
inhomogeneous equations like (3.44) or ( 3.45| ) in general. 
The general form of such an equation is 

T n [z,d z ] g(z) = i n (z) , (Fl) 

with i n (z) the inhomogeneity and T n [z, d z ] the operator 

* n , - .. 



T n [z,d z ] 



dz 2 



z 

2 



d 
dz 



n 
2 



We search for a particular solution g(z) of Eq. (|Fl|). A 
particular solution of the homogeneous equation (i n (z) = 
0) can be expressed by Hermite polynomials: 

f n [z,d z ]h n (z)=0 (F3) 
ho(z) = 1 , hi(z) = v2z , h 2 (z) = 1 — 2z etc. 



G7 



The ansatz g{z) = h n {z) u n (z) reduces (Fl) to an equa- 
tion for d z u n {z) of first order: 



, d d \nh n {z) 
T n g = z h n (z) [ — + 



1 — 2z\ du n {z) 



z h„(z) 



dz 



d z (M n (z) d z u n {z)) 



2z 



dz 



(F4) 



M n (z) 

where in the last line we introduced the integrating factor 



M n {z) = h 2 n {z) V-z e 



(F5) 



Identify now T n g n = i n , integrate twice, and substitute 
M n by the full expression. A special solution of (Fl) then 
reads 



g{z) = h n (z) 



dx h2t \ r -x > ( F6 ) 



where the integration constants a and b are free. If we 
in particular choose b = oo, the integrated exponential 
e x ~ y cannot exceed unity, and g^{z) can at most diverge 
algebraically, if the integrated inhomogeneity i n {z) is al- 
gebraic. 



Integrating Eq. (3.55 ) for go(z) as in (F6) with b — oo, 
we find for the algebraic divergence of go(z) for large z: 



9o{ z ) ~ 3a z lnz as z 



(F7) 



while the solution of the homogeneous equation diverges 
only as h^z) ~ z. For determining the small z expan- 
sion of ([Fq), it must be noted that the factor h n {x)~ 2 is 
singular at the zeroes of h n {x). Hence, ( |F^ ) needs to be 
evaluated separately in each interval between the zeroes 
of h n (x). This can be done by a proper choice of a. It can 
be shown, that the results in each interval join smoothly. 



APPENDIX G: ALGEBRAIC CONVERGENCE 
AT THE PUSHED/PULLED TRANSITION 



In Section III we have analyzed equations, that are 
within the pulled regime. We here analyze equations, 
that are at the pushed/pulled transition. Leading edges 
of fronts within the pulled regime have the form $* = 



« + 0) e"« oc £ e-« (£ > 1), cf. (p7l7b . Leading 



edges of fronts within the pushed regime are given by 



$t (x g-A+^'K, cf. ( p.23| ). Leading edges of fronts at 
the pushed/pulled transition accordingly behave as 



$* = /? e - e for£»l , \+{v*) = \* = l 



(Gl) 



For our example nonlinearity (1.10), fronts are within the 
pulled regime for e > {n+l)/n 2 and at the pushed/pulled 
transition for e = {n + l)/n 2 . The analysis below can 
again be extended to more general equations along the 
lines of Section |V|. We will come back to this at the end 
of this Appendix. 



At the pushed/pulled transition, the spectrum of linear 
perturbations is still gapless, and convergence therefore 
is algebraic. On the other hand, the form of the lead- 
ing edge played a crucial role in determining the veloc- 
ity corr ections X. Compare our qualitative disc ussion in 
Section HI A 1 . The leading edge behavior ( |Gl[ ) immedi- 
ately lets us expect, th at no w v(t) = 2 — l/(2i) + . . ., in 
contrast to (3J5) and ( |3.66 ) for fron ts within the pulled 
regime, and in agreement with ( 2.46| ) for the spreading of 
perturbations under the linearized equation. Intuitively, 
we can argue, that the slower convergence of fronts within 
the pulled regime is due to the leading edge having to pull 
the interior part of the front along. This also makes the 
leading edge flatter. The quicker convergence of fronts at 
the pushed/pulled transition and in the linearized equa- 
tion then resembles the fact, that the leading edge and 
the interior part of the front "impose the same speed" . 

Let us now do the explicit convergence analysis for 
fronts developing from initial conditions steeper than e~ x 
for x ^ 1 and aproaching (Gl) for large time s. Th e anal- 
ysis of the interior is identical with Section III B , where 
constants c.j are yet undetermined. When expanding 
the interior shape towards the leading edge as in Section 
III C , the inhomogeneities created by $* flGl| ) are differ- 
ent, because now a = 0. The differential equations for 
the V'f result from ( 3.34 ) with a = 0, 7 = j3 and start 
with 



= ci/3 



9|^3 = C3/3 



(G2) 



9^2 = [-l + ci(l-0 € )]Vi +c 2 /3 + o(e^ x ) etc. 
Integrating and resumming, we now find for £ 3> 1 
i/f - 13+ (G3) 



2! t 



t \ t 



c^e x . „(tx 



2! t3/2 



O 



t 3/2 



ci( Cl -l)/3££ , c 1 {c 1 8-5-c 1 0)8 L {& 



4! t 2 



+ 



3! t 2 



t 2 



4! < 5 / 2 



Us/2 



Here 6 is an unknown integration constant fixed by con- 
dition (|3.9|). We will see below, that it is not involved 
in fixing the velo city, just as also the sublcading j3 for 
the leading edge (U|) within the pulled regime is n ot 
involved in fixing the velocity, cf. calculation till (3.65). 

Again, for £x 3> \ft we have to reorder the expansion 
in powers of ^fz = \J /{At) and 1 / y/i, and find 



2! 



4! 
2! 
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d(d<J - S - ci/3) (4z) 3 / 2 C|(2 Cl - (4z) : 



3! 



4! 



(G4) 



The structure of the expansion is the same as in ( 3.42 ) 
except that now the leading order term is of order t°: 



G(z,t) = e z V =go(z) + 



(G5) 



The equations of motio n for the l eadin g and subleading 
term are derived from (3.43) - (3.45) through putting 
q -i = 0. For go we find now the homogenous equation 



1 \n 1 

--z d z --- Cl 



9o = 



(G6) 



Just like (l3~43|) was solved by (|3.52f) , we now solve flGq ) 
with 



ci = -y j £to(*0 = P ■ 



The equation for gi is now, cf. (3.45) and (G7) 



5i = /8 



(G7) 



(G8) 



Again a special solution of the inhomogeneous equation 
can be found, and the general solution contains the con- 
stants of integration ki and l±: 



9± 







2c3 



>z (l)„_i z n 



E 



-1 1 



+ k h M [T'r z ) +l i Vz 

2 = 1 2pC3+ki+h^+0(z) 



4z 



2.2 



(G9) 
(G10) 
(Gil) 



Comparing ( G10|) to the o rder l /y/i in (p4|) and impos- 
ing proper convergence of (Gil) for z — > oo, we find 



2/3c|+/ci=0 , Zi=-<5 , /Sx/tt + 2fei = . (G12) 
With these constants, the velocity correction ci is 



In summary, we find for the convergence to a front at 
the pushed/pulled tr ansit ion, whose leading edge accord- 
ingly takes the form (Gl), that the velocity correction is 
given by 



(G15) 



In the interior, i.e., for £x V~t, the front is given by 
( |3.31 ) just like a front within the pulled regime. In the 
leading edge, where £x 3> y/i, the front is given by 

=e-**-&/(">G^|,i) , (G16) 



where 



G(z,t) = /3 + 



Vt 



o 



(G17) 



The extension along the lines of Section [v| to more gen- 
eral equations is straightforward. The general expression 
for X(t) is 



X 



2X*t 



1 

2X* 



O 



t 2 



(G18) 



but the subleading function g i (z) will depend on the ad- 
ditional terms i n the expansion, just like the subleading 
go{z) in Section | 



APPENDIX H: MULTIPLICITY OF FRONTS 

AND LINEAR EIGENMODES FOR 
REFLECTION SYMMETRIC EQUATIONS OF 
FIRST ORDER IN TIME 

The generical multiplicity of uniformly translating 
fronts can be determined by counting arguments anal- 
ogous to those performed in Section II B. Uniformly 



translating solutions of ( |5.l| ) can be understood 

as a heteroclinic orbit in TV-dimensional phase space be- 
tween fixed points characterized by = 1 at £ — > — oo 
and $„ = at £ — > oo. For a linear pert urb ation 
5 = 1 — $> v about the fixed point <j> = 1 from (5.2), we 
get the equation 



£„(-oo) 5 + 0(S 2 ) = , 



(HI) 



C3 



4 ' 



and for qi we find 

(3^ 



9k 



2 

- Sy/z 



I 1/ | — ,-,z 
2 2' 



(G13) 



2 ^ (l)n-l 2 r 



E 



7T ' if) 7l! 

n— 1 \2/n 



(G14) 



which is a linear ordinary differential equation with con- 
stant co effic ients with the linear operator C being de- 
fined in ( [5.3] ). The same is true for a linear perturbation 
= + S of the fixed point <j> = 0, which solves 



C v {oo) 6 + 0(S 2 ) = 



(H2) 



In linear order of 5, each of these equations has N solu- 
tions e - A "^ £,n=l,...,N. 
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Let us restrict the analysis to real equations which are 
isotropic in space, i.e., where ( |5.41 ) is invariant under 
x — > —x. Such equations are even in d x , so N needs to 
be even. Accor ding to ar gum ents presented in Appendix 
A of §7|, Eqs. @ and @ for v > will have N/2 + 1 
eigenvalues A„ with positive real part and N/2 — 1 ones 
with a negative real part, if the state, about which we 
linearize, is linearly unstable against a range of Fourier 
modes. If it is stable, we will have N/2 eigenvalues with 
positive real part and N/2 ones with a negative real part. 
We assume = 1 to be stable, so at £ — > — oo there are 
N/2 directions in phase space with negative real part 
of A, that need to be excluded. If <j> = is unstable, 
we have only N/2 — 1 bad eigendirections at £ — * oo. 
We then generically have a front connecting these fixed 
points for arbitrary values of v. If, however, the state 
<f> = is metastable, there are N/2 bad eigendirections 
at £ — > oo. Then also v needs to be tuned to find a solu- 
tion. So for fronts propagating into unstable states, we 
generically have a front solution for a continuum of 
velocities, while for fronts into metastable states, there 
are solutions only for discrete valu es of v, in general- 
ization of the arguments from Section II B . 

The multiplicity of linear perturbations is determined 
along the same lines. We again decompose the linear 



perturbations n ( 5.3CQ into t) = ?7 CT (0 e at by sepa- 



ration of variables. The rj a then solve the o.d.e. 



(H3) 



For counting the generic multiplicity of solutions, we need 
to linearize the equations about £ — > ±oo, which amounts 



to a problem equivalent to (HI) and (H2), except for a 
shift of the constant contribution of £■«(£) by a. For 
fronts propagating into unstable states, we in general ex- 
pect a continuous spectrum a of linear perturbations at 
least in s ome finite interval of a, in generalization of Sec- 
tion Inc. 



APPENDIX I: STRONGLY HETEROCLINIC 
ORBITS AND CHANGE OF STABILITY AT V t 

According to the counting argument from Appendix 
||, the front $*(£) propagating uniformly with velocity 
v* does exist. The question is now, whether it is stable 
and whether it will be approached by steep initial condi- 
tions. In particular, we want to analyze initial conditions 
4>(x,0), that are steeper than e~ x x in the leading edge. 

This amounts to the question, whether in the spectral 
decomposition r\ a (H3) of a generic 4>{x, 0) — $*(#), there 
are destablizing modes with a < 0. As in Section HC, 



the contributing modes in general will all decay at least 
as quick as <!>* in the leading edge. The leading edge 
properties of the r\ a in general will depend smoothly on 
a, just as in (D14), so generically $* will still be stable 
against all perturbations, that in the leading edge decay 
quicker than $*. 



An exemption is again the generalization of A a = 
from (|DS| ). For an equation of order N with (p — un- 
stable, there are N/2 + 1 exponents A n (a,v) > 0. The 
leading edge will be a superposition of all the exponen- 
tials 



N/2+1 



-An* 



as £ > 1 



(II) 
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The condition = 0, where Ai is the smallest one of 
the positive A„, fixes a discrete set (which can be either 
empty or not empty) of negative cr's whose eigenfunctions 
rja have a steepness in the leading edge larger than $*. 
The stability of the pulled front $* thus again depends 
on the "strongly heteroclinic" perturbations. 

If there are strongly heteroclinic perturbations, that 
destabilize the pulled front propagating with velocity v* , 
then there will be a steeper and quicker front & , which 
can be constructed as a strongly heteroclinic orbit of 
zero mode then again is a strongly 

heteroclinic perturbation, and as discussed in Appendix 
[l5]on the stability of front solutions, we can conclude that 
the quickest of all strongly heteroclinic orbits cannot be 
destabilized, so it will attract all sufficiently steep initial 
conditions. 

We conclude, that Table IV generalizes to higher or- 
der equations, which form uniformly translating fronts, 
if we only appropriately adjust the explicit definitions of 
the velocities v and steepnesses A. 



APPENDIX J: RELATION BETWEEN THE 
GENERALIZED DIFFUSION CONSTANTS D N 
AND THE DISPERSION RELATION 



If we use the expansion fl5.q ) for the dispersion relation 
u>(k), we get 



N 



V = e x *Z ^a ro cf -v*dt\ e" A *« , 
\m=0 / 

N 

= ^a m (%-A*r-u*(d s -A*) , 



m=0 
N m 

EE' 



' n\{m — n)\ 



^ ^*^(m— n) Qr, 



-v*{d t - A*) 



N ( F) n N \ i 

n=0 \ m=n / 



V(%-A*) 



(Jl) 



This immedia tely y ields the expansion (5.27) with the 
identification (|5.28 ). 
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APPENDIX K: EDGE ANALYSIS OF 
UNIFORMLY TRANSLATING PULLED FRONTS 
WITH M = 1 



We analyze the leading edge representation ( |5.37 ) for 
a uni formly translating front whose equation of motion 
( 5.41 ) is of arbitrary order N in space and of first order 
in time M = 1: 

d T rl> = ( $ + d " d <^j *l> + Y(d c -l)il>- (Kl) 

We g ener alize the leading edge analysis from Sections 
ilTc| and pID| , 
With the notions and ansatz 



N 



Cn 



71—3 



i ! ((Y,T) = a( Y +(3+^ + ^ + ^ + ... , (K3) 



n=2 

t/>i V's 

T T'' 



(K2) 



the expansion of the interior in the region of (y 3> 1 at 
the crossover towards the leading edge reads 



2?Vi = Ci(" +7) . J = l3-a 



(K4) 



in generalization of ( |3.34 ). These equations can be inte- 
grated explicitely. The result can be written in leading 
edge variables z = £y/(4rr) as 



,1/2 , ^(48)8/2 



if> = V^a (4z) i/z + 



3! 



Ci(^-q(l + d 3 ))(4z) 
2! 



0(1/V7) 



(K5) 



This generalizes the results of Section [II C and supplies 
us with the small z expansion of the leading edge function 



V>(C,r) 



e~ z G{z,t) 



4r 



(K6) 



G solves [compare Eq. (3.41)] 



zd 2 



— — z ] d z — — — rd T — C\ 
2/2 



G 



G 



O 



1 



^ dz -l) 2 + z(d z -lf 



G 



(K7) 



where we wrote all operators of order t° on the l.h.s. 
of the equation and the operators of order t -1 / 2 on the 
r.h.s. 

With the ansatz 



gi (z) 

G{z,t) = Vr g=±{z)+g (z) + -S— + 



(K8) 



as in (3.42), we find that q-i (z) solves again ( [3.43 ), so 
we copy from Section III D, that 



_3 

G\ = — , g^.(z) = 2ay/z 



(K9) 



For go(z) we then find instead of (3.55) 
1 



--') 2 + [ |- : ] ih - 1 



9o 



(K10) 



2a 



3(l + rf 3 ) , r- 3 a ' 
^ \- c^z — - z + d 3 (z — dz) 



A special s olutio n of the inhomogeneous equation is now 
instead of ( |3.5S| ): 



9o P ( z ) = 2a 



3 (1 + ds) 



r 2c3 — - F 2 (z) — d 3 z 2 

4 2 4 



(Kll) 



with F 2 (z) from ( 3.56 ). The general solution is 

g {z) = g s p (z) + k (l-2z) + loV^M(^-,^,z^ 

.^1 ^ (1 + + ^ + +1q ^ + 0{z) 

- ^_ a ^+ l A z -m e ^ (K12) 



Note, that ^3 7^ does not cause any divergences at 
z — > 00. It only shifts the constant contribution at z — > 0. 



Suppressing the divergence at z — > in (K12), and 
comparing its small z expansion to ( |K5| ) yields again 



C3 = 



30? 



(K13) 



and 



g (z) = P (1 - 2z) + 3a(l + d 3 )z - 2ad 3 z 2 



71 



APPENDIX L: LEADING EDGE PROJECTIONS 
FOR COUPLED EQUATIONS: AN EXAMPLE 

As a simple illustration of the vario us ques tions related 
to the projection discussed in Section VE3, we consider 
two coupled F-KPP equations, 



d t <t>i = d 2 x <t>i + <f>i - <t>{ , 

d t fo = Ddlfc + 4> 2 -4> 3 2 



(LI) 



K<h. 



The dynamics of this set of equations for fronts propa- 
gating into the state 4>\ — <\>i — with steep initial con- 
ditions, is of course immediately obvious: when K = 0, 
the two equations are uncoupled, and fronts in the first 
equation propagate with speed v\ — 2, while those in the 
second equation propagate with speed v% = 2\f~D. The 
dynamics of 4>i is always independent of that of 02 , even 
for K 7^ 0, so for K > and D < 1, the dynamics of the 
coupled equations amounts to a normal F-KPP </>i front, 
with relaxation given by our usual expressions. This front 
entrains a front with speed v = v\ = 2 in (f>2 . For D > 1, 
the 4>\ and <p2 fronts keep on propagating with different 
speeds. We consider the case D < 1 and make a lead- 
ing edge transformation <j>i — e~^ipi, 02 = e~^ip2 (where 
A* = 1) to the frame moving with velocity £ = x — v*t. 
The linearized equations then become 



d t ipi = d^ip! , 



(L2) 



dti>2 = Ddlfa + 2{D - 1)5^2 + {D- 1)^2 + KiPt . 

The matrix S_* (g, ft) of the linearized equations is in this 
case 



s*( q ,n) 



itt-q 2 
K 





ift - q 2 + J(q) 



(L3) 



where J(q) = (D — 1)(1 + 2iq — q 2 ). Since the element 
Si2(<7, ft) = 0, the eigenvalues u{ and u\ are simply the 
diagonal element of S_*(q, ft), u*(q, ft) = ifl — q 2 and 
u^{q, ft) = ift — q 2 + J{q)- However, the eigenvectors are 
not both along the ipi and tp2 axis. Indeed, we have in 
the notation of |V1 



HI (?) = ( -K/J(q) ) ^ = (!. °) . ( L4 ) 

U$(q) = (K/J(q) ) l) , (L5) 



£S = , ! 



The appropriate saddle point is ft 
J(0) = (D- 1), we have 



£3(0) = 



-K/(D-1) 



0, and since 



(L6) 



The fact that the second component is nonzero just ex- 
presses the fact that the variable tp 2 is entrained by the 
leading edge in ipi. We can now illustrate our asser- 
tion that different choices of projection lead to differ- 
ent dynamical equations for the projected leading edge 



variable i\j v , but that the universal results from Table 
II are independent of the particular choice of projection. 
Clearly, one obvious intuitively appealing choice is to take 
if} p = ipi , since the ipi dynamics is independent of that 
of ip2- In this case, the dynamical equation for i\) v is 
nothing but the single F-KPP equation, and all the re- 
sults for this equatio n carr y over in detail. Likewise, the 
choice ip p = 7ri(g, ft) ( |5.86|) leads to the linearized F-KPP 
equation for ip p since u\(q, ft(<?)) = gives the dispersion 
relation of the F-KPP equation. However, this choice 
is more formal than practical, since the direction in the 
vector space (^1,^2) is not fixed, but depends on the 
variable q which influences the dynamics. A more practi- 
cal choice for the coupled variables would be to take ip p as 
the component along U_l (0) , as this corresponds to a fixed 
ratio of ip-L and ip 2 - Since U$ ■ U_l = K(J(q) - J(0) = 
—2Kiq/{D — 1) + 0(q 2 ), the projected equation for in 
this case picks up a third order derivative term D^,d^ip v , 
amoungh other ones. 

Thus, we observe in this particular example, that in- 
deed the universal results from Table II on velocity and 
shape relaxation are independent of the choice of projec- 
tion, while the subleading contribution go(z) in the lead- 
ing edge is universal in the sense, that it is independent 
of the precise initial conditions, but it does depend on the 
direction of projection. 



APPENDIX M: PINCH POINT VERSUS SADDLE 
POINT ANALYSIS 

In this appendix, we briefly discuss the major differ- 
ences and similarities between the saddle point and the 
pinch point approach for evaluating the integral 



T 



17+00 dw_ 
■v_^ 2tt 



dk 
2^ 



D ik^ — i{uj—vk)t = 



M (k,to) 



U m (k,LO) 



(Ml) 



from Eq. (5.79) on a given branch m. Here 7 > 
needs to be large enough, that the integrand is an- 
alytic along and above the path of u> integration in 
the complex uj plane. We introduced the abbreviation 

M (fc, uj) = U_ m (k, uj) x U_ m (k, uj). In the moving frame 
£ , it obviously is convenient to transform to the variable 
Q = uj — vk, and to introduce 

u^(fc, ft) = u m (k, ft + vk) = u m (k, uj) , ft = uj — vk . 

(M2) 

The characteristic equation 

Um (k,u m (k)) = ^ u v m (k,Q m (k))=0 (M3) 

defines the dispersion relation uj m (k) or ft m (fc). The in- 
tegrals are now of the form 
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r 7+o ° an 

'17-00 ^7T 



dk iki _ mt KJk,n + vk) 

2n 



(M4) 



The "sadd l e po int" type approach, that we follow in 
Sects. V C - VE of this paper, is based on first evaluating 
the f2 integral by closing the fl contour in the lower half 
plane for t > around the simple pole oc (f2 — f2 m (fc)). 
The integral then yields 



dk 
2^ 



o4 *-«M* )t KJk,n m (k) + vk) 



, (M5) 



where 7 needs to be larger than maxfc rca i (Im O m (fe)). 
From here on, t he sa ddle point analysis proceeds essen- 
tially as in Sect. VD: the k contour is deformed continu- 



ously such that it passes through a saddle point of fi TO (fc) 
that allows for a steepest descent evaluation of the k in- 
tegral. A saddle point is a double root in k of u^ n (k 1 O), 
so that 



u v m (k,n)\ sp = o 

and 

d k u v m (k,n)\ sp = o ■ 



^m{k S p) 



vk s 



(d k + vduj) u m (k,u>)\ sp 

d k u m (k,Lj)\ 

v = — — . 

d u u m (k,u})\ 



(M6) 

} 

= 
(M7) 



By expanding about such a saddle point, we get for large 
t to leading order 



7 — 



K(k,w) 



idujU m (k,uj) 



sp 



with the diffusion constant 



-i (d k + vd^) 



D sp = 



sp 



-i dlul 



k "m 



2 a 



l sp 



2 cL<J 



(M8) 



(M9) 



The remaining integral over real q = k — k sp is a sim- 
ple G aussia n integral of the form discussed previously in 
Sect. VC 1. As before, we are in the comoving frame, if 



Im tt sp = 



Im Lo m (k 



sp) 



Im fc. 



(M10) 



Differentiating the dispersion relation u(k,uj m (k)) = 
with respect to k: dku(k,uj m (k)) — 0, and comparing to 
(M7), we can immediately identify 



duj m (k) 



Ok 



(Mil) 



From dlu(k,uj m (k)) = and (M9), we get 



D 



id 2 u> m (k) 



2dk 2 



(M12) 



sp 



■7S|) t he branch m with the largest velocity 
(|5.8l|) immediately results. 



Choosing in ( 
v sp = v* , Eq. 

If the denominator of an integral like (M4) contains a 
product of characteristic functions f| m=1 u^ n (k, fi), then 
each factor wi ll co ntribute with its pole and 

yield an integral as in (M5), so that the t otal integral 
amounts to a sum of M integrals of the form (M5). Again 
the dominating contribution for £ fixed and t 1 will be 
the one with the largest velocity v sp through which the 
contour of k integration can be deformed. 

Th e pinch point analysis ]53] is based on evaluating 
(Ml) by a different order of the integrations, i.e., by first 
closing the k contour to get k = fc(O) and then evaluat- 
ing the remaining f2 integral. (For £ > 0, the k contour 
must be closed in the upper half plane.) As discussed 
most clearly by Bers j58j, this is done as follows. 7 in 
(M4) has to be large enough to lie above the maxima of 
the dispersion relation Q m (k) for real k. When O varies 
along the integration path, the poles in the k plane move. 
Now when 7 is lowered sufficiently, that it approaches the 
maximum of the line Q m (k) traced out by the real k val- 
ues, a pole in the k plane will approach the real k axis. 
When that happens, the k contour can be continuously 
deformed to avoid this pole. This in turn allows one to 
lower the value of 7. This process can continue until two 
poles in the k plane approach the k contour from oppo- 
site sides, and "pinch off" the k contour at a particular 
value of f2*. Clearly, that point corresponds to a dou- 
ble root, since for that given value of fl the two k roots 
coincide. When the k contour is closed, this point gener- 
ates a branch-cut in the Q plane, since near f2* we have 
k-k* = ±yj{n - Sl*)/D. When the fi contour is sub- 
sequently closed in the lower half f2 plane, these branch 
poi nts th en gene rate the usual leading asymptotic behav- 
ior (|t|), ftjjj ). 

In both approaches, there are conditions for a saddle or 
pinch point to be dynamically relevant; these arise from 
the global properties of the dispersion relation cu(k). In 
the saddle point approach, only saddle points that will 
dominate the k integral along the deformed contour of 
integration, are relevant for the dynamics. Pictorially, 
a saddle point that obeys this conditions is located be- 
tween "valleys" of Im u>(k) in the direction of real k that 
are not completely separated by "ridges" from the real k 
axis. In this formulation, the condition Re D > nat- 
urally comes out. If there is more than one such saddle 
point, the one with the highest velocity v* determines the 
asymptotic spreading velocity. In the pinch point formu- 
lation, the condition usually mentioned is that the poles 
in the k plane "pinch off" the k contour, while the condi- 
tion Re D > is usually not mentioned, but it is actually 
hidden in the formulation as well: it just expresses that 
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the pinch point is associated with a point of the disper- 
sion relation, where the growth rate is maximal. In fact, 
the examples discussed on pages 466, 467 in Q for so- 
lutions of the saddle point equations which are no pinch 
points, are just cases where Re D < 0, i.e., solutions 



a) f(0)<0 




b) f(0)>0, v>v 
a 



which are excluded by a saddle point formulation as well. 
In the pinch point formulation, the improper solutions of 
the saddle point equations correspond to solutions where 
two poles in the fc-plane do not "pinch off" the deformed 
£;-contour, but instead just merge by themselves. 

c) f'(0)>0, v=v* 
a 



v74-f'(0) 




A + (a,v) 



FIG. 15. Steepness A(ct, v) (D1C) versus decay rate a of linear perturbations r\ a (D7), (D8) of a given front $„ with velocity 
v > v* . The solid curve denotes real A, the dotted curve the real part of complex A. X±(v) and A* are the steepnesses of and 
of the zero mode 770 = d^ v - They are marked by circles on the A axis. The generic steepness of a front & v with v > v* is X—(v), 
while in the particular case of A v = 0, it is X + (v). The continuous spectrum of a is denoted by a fat solid line on the a-axis, the 
interval in which there may be discrete eigenvalues a by the fat dotted line. The continuous spectrum within the Hilbert space 
of 7i v exists only at v 2 /4 — /'(0) < a. The continuous spectrum for — /'(0) < a < u 2 /4 — /'(0) is on the A_-branch. There 
might be discrete solutions characterized by A„ = 0. They lie on the A+-branch, might exist for all a < n 2 /4 — /'(0), and need 
to be constructed, a) The front propagates into a metastable state (/'(0) < 0). Its steepness is X+(v). It is stable against all 
linear perturbations with A < X+(v). The discrete spectrum of steep perturbations with A > X+(v) needs to be investigated. 

b) The front propagates into an unstable state (/'(0) > 0) with velocity v > v* . It is stable against all linear perturbations with 
X—(v) < A < X+(v), it is unstable against the continuous spectrum of very flat perturbations with < A < X-(v), which might 
be excluded by the initial conditions. The discrete spectrum of steep perturbations with A > X+(v) needs to be investigated. 

c) The front propagates into an unstable state (/'(0) > 0) with velocity v = v*. The discussion is as for (6) after identifying 
A±«) = A*. 
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